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FOREWORD 


The  study  reported  herein  was  conducted  in  support  of  DA  Project 
4A663712D860,  "Military  Geographic  Systems,"  Task  04,  "MGI  Data  Base," 

Work  Unit  04,  "Automation  of  Performance  Prediction  Models,"  sponsored 
by  the  U.  S.  Army  Engineer  Topographic  Laboratories.  The  work  was  per¬ 
formed  by  personnel  of  the  Terrain  Analysis  Branch  (TAB),  Mobility  and 
Environmental  Systems  Laboratory  (MESL),  U.  S.  Army  Engineer  Waterways 
Experiment  Station  (WES),  under  the  general  supervision  of  Messrs.  W.  G. 
Shockley,  Chief,  and  S.  J.  Knight,  formerly  Assistant  Chief,  of  the  MESL, 
and  under  the  direct  supervision  of  Mr.  W.  E.  Grabau,  Chief,  TAB;  Mr,  M.  P. 
Meyer  was  Project  Manager.  Computer  programs  were  developed  by  Miss  M.  H. 
Smith;  Mr.  C.  A.  Carlson  assisted  in  developing  the  model.  This  report  was 
prepared  by  Miss  Smith  and  Mr.  Meyer.  Mr.  A.  Vazquez  and  Mrs.  M.  P.  Terry 
assisted  in  computer  programming  and  compiling  information. 

COL  Ernest  D.  Peixotto,  CE,  was  Director  of  the  WES  during  this  study 
and  preparation  of  the  report.  Mr.  F.  R.  Brown  was  Technical  Director. 
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AUTOMATION  OF  A  MODEL  FOR  PREDICTING  SOIL  MOISTURE 


AND  SOIL  STRENGTH  (SMSP  MODEL) 

PART  It  INTRODUCTION 
Background 

1.  The  study  reported  herein  Is  a  part  of  the  work  to  automate 
performance  prediction  models  of  immediate  military  interest  for  inclu¬ 
sion  in  an  automated  Military  Geographic  Intelligence  data  base  system. 

Completed  performance  prediction  models  include  the  Airfield  Construction 

1  2 
Effort  Model  in  FY  70,  the  Cross-Country  Locomotion  Model  in  FY  71, 

3 

and  the  Helicopter  Landing  Zone  Model  in  FY  72.  This  report  describes 
an  automated  model  for  predicting  soil  moisture  and  soil  strength  (SMSP 
model) . 

2.  Methods  for  predicting  soil  moisture  and  strength  (defined  in 

terms  of  cone  index  (Cl)  and  rating  cone  index  (RCI)^  were  developed 

previously  at  the  U.  S.  Army  Engineer  Waterways  Experiment  Station  (WES) 

for  use  in  predicting  the  trafficability  of  fine-grained  soils  in  multi- 

4  5 

pass  cross-country  vehicular  operations.  ’  For  that  application,  strength 
predictions  were  required  for  the  15-  to  30-cm  layer,  the  critical  layer 
influencing  performance  of  many  military  vehicles.  Moisture  contents  of 
the  0-  to  15-  and  15-  to  30-cm  layers  were  also  provided  for  use  in  the 
prediction  process.  In  a  subsequent  study  of  worldwide  strength  condi¬ 
tions  of  surface  materials,^  soil  strength  relations  were  extrapolated  to 
the  0-  to  15-cm  layer.  These  relations  f’ave  been  incorporated  in  the 
present  model.  An  earlier  version  of  the  model  was  used  to  predict  long¬ 
term  trafficability  conditions  in  the  vicinity  of  Saigon,  South  Vietnam, ^ 
to  provide  data  for  construction  of  soil  moisture  and  strength  records  in 
synthalogous  environments,  to  provide  RCI  data  for  use  in  an  analytical 
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procedure  for  quantitatively  comparing  the  similarity  of  terrain  sites, 

and  to  provide  moisture  content  data  input  to  a  model  that  predicts  induced 

radiation  activity  of  surface  soils  following  nuclear  detonations.***’** 

3.  The  model  presented  herein  has  been  used  to  predict  moisture 

12 

contents  for  the  Seismic  Sensor  Performance  Prediction  Model,  and  strength 

2 

values  for  input  to  the  Cross-Country  Locomotion  Model  and  the  Army  Materiel 

«  e  4  / 

Command  (AMC)  Model  for  Predicting  Cross-Country  Vehicle  Performance.  ’ 

It  can  also  be  used  to  provide  data  that  can  be  converted  to  California 
Bearing  Ratio  (CBR)  and  soil  moisture  condition  for  use  by  the  Airfield 
Construction  Effort  Model*  (see  Appendix  G),  The  model  also  includes 
an  option  to  store  data  for  use  in  another  model  that  predicts  snow  depth, 

snow  density,  and  frost  and  thaw  penetration*^  (herein  referred  to  as 

v  * 

the  Freeze-Thaw  model) . 

Description  of  the  Model 

♦.  The  SMSP  model  is  a  composite  of  soil  moisture  and  soil  strength 
prediction  methods.  The  prediction  method  for  soil  moisture  is  a  daily 
bookkeeping  procedure  wherein  soil  water  is  added  as  a  result  of  precipita¬ 
tion  or  subtracted  as  a  result  of  depletion  (resulting  from  evaporation, 
transpiration ,  and  drainage) .  Predictions  are  influenced  by  terrain  and 
weather  conditions. ^C’  <*’~€’  Soil  strength  values  are  predicted  from 

daily  soil  moisture  predictions  by  using  known  relations  or  relations  that 
are  determined  by  certain  soil  properties.  J’  The  model  combines 

*The  program  for  this  model  is  on  file  at  the  WES. 


2 


SUMMARY 


The  soil  moisture  strength  prediction  (SMSP)  model  is  a  composite 
of  the  methods  developed  at  the  U.  S.  Army  Engineer  Waterways  Experiment 
Station  for  predicting  daily  soil  moisture  contents  and  strengths  (in 
terms  of  cone  index  and  rating  cone  index)  of  soil  layers  at  depths 
of  0-15  and  15-30  cm.  Information  required  by  the  model  includes  soil 
moisture  accretion  and  depletion  relations,  field  maximum  and  minimum 
soil  moisture  contents,  moisture  content  at  start  of  prediction,  soil 
dry  density,  soil  moisture-strength  relation,  daily  rainfall  amounts, 
and  minimum  rainfall  amount  required  f<  \  ,-_-retion.  This  information 
can  be  obtained  from  one  or  more  of  .hi--1  sources:  (a)  directly  from 
measurements  at  a  specific  location;  (b)  indirectly  from  estimated  or 
averaged  data  derived  from  field  measurements,  literature,  or  empirical 
equations  built  into  the  model;  or  (c)  indirectly  from  a  surface  compo¬ 
sition  group  classification  that  closely  follows  the  Unified  Soil 
Classification  System. 

The  computer  program  for  the  model  Is  written  in  Fortran  IV  conver¬ 
sational  mode  for  use  on  a  teletype  connected  to  a  Honeyvell-GE  (General 
Electric)  440  computer.  Output  data  are  stored  in  permanent  files  for 
use  by  other  performance  prediction  models,  for  printing,  or  for  input 
to  plotting  programs. 

The  main  text  of  the  report  includes  a  discussion  of  the  structure, 
operation,  use,  limitations,  and  mathematics  of  the  model. 

Appendixes  A-G  include  detailed  flow  charts  and  listings  of  the  computer 
program;  listings,  organization,  and  format  of  input  data;  examples 
of  prediction  runs  and  graphic  displays  of  results;  and  procedures  for 
converting  output  data  to  terms  required  by  the  airfield  construction 
effort  inode J . 
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both  prediction  methods  and  outputs  daily  soil  moisture  content  and/or 

daily  soil  strength  values  (Cl  and  RCI)  for  soil  layers  at  depths  of  0-15  cm 

* 

(designated  as  layer  1}  and  15-30  cm  (designated  as  layer  2). 

Input  data 

5.  Terrain  and  weather  data  are  required  as  inputs  to  the  model. 

6.  Terrain  data.  The  terrain  data,  which  are  referred  to  as  control 
data  in  the  automated  program,  include  minimum  rainfall  amount  required 
for  accretion  and  the  following  for  each  layer  of  a  given  soil: 

a.  Accretion  relations. 

b.  Depletion  relations. 

jc.  Field  maximum  and  field  minimum  soil  moisture  contents. 

4.  Soil  moisture  content  at  the  start  of  prediction. 

e.  Soil  moisture-strength  relation. 

f.  Dry  density. 

7.  The  terrain  data  may  be  one  or  a  combination  of  three  different 
types  as  follows: 

a.  Specific  data.  These  data  Include  the  above  prediction  rela¬ 
tions  and  soil  property  data  that  are  derived  from  measurements  taken  at 
a  specific  location.  Use  of  these  data  provides  the  most  accurate  predictions 
of  soil  moisture  and  strength  for  a  particular  site  under  a  given  set  of 
terrain  and  weather  conditions.  Procedures  for  derivation  of  specific 


* 


Each  layer  is  actually  15. 24  cm  (6  in.)  thick.  In  converting  data  from 
metric  to  English  units  the  program  uses  the  actual  thickness. 


prediction  relations  are  given  in  reference  Ac;  the  preparation  of  the 
input  data  is  discussed  in  Part  III  of  this  report. 

b.  Estimated  data.  These  data  refer  to  averages  or  estimates 
derived  from  (a)  field  measurements,  (b)  puhlished  literature  containing 
data  on  terrain  analogous  to  the  test  area,  or  (c)  equations  developed 
from  soil  property  data  from  a  large  number  of  sites.  These  equations, 
called  tentative  average  (TA)  relations,  are  incorporated  in  the 

model  and  require  measured  or  estimated  values  of  soil  properties  as  input. 
The  methodology  for  deriving  these  equations  is  discussed  in  references 
&e  and  4j.  The  moisture  equations  and  variables  are  shown  in  table  A  and 
the  strength  equations  and  variables  are  shown  in  tables  6  and  7.  The 
mathematics  of  the  equations  are  discussed  in  Part  IV  of  this  report;  the 
preparation  of  the  input  data  is  discussed  in  Part  III. 

c.  Surface  composition  group  data.  These  data  are  also  averages 
or  estimates  that  have  been  derived  from  sources  similar  to  those  described 
above.  However,  here  a  set  of  control  data  has  been  generated  for  each 

of  a  number  of  surface  composition  groups.  The  groups  correspond  closely 
to  soil  classes  of  the  Unified  Soil  Classification  System  modified  by 
the  addition  of  nonsoil  water,  pavements  and  structures,  and  rock  to  allow 
for  characterization  of  the  entire  surface  area.  Each  soil  group  has  all 
the  Information  needed  for  operation  of  the  model  stored  in  a  file  subject 
to  call  by  a  group  code  number.  Materials  of  the  surface  composition  groups 
and  the  corresponding  group  code  numbers  are  shown  in  table  1.  The  surface 
composition  group  can  be  used  when  specific  data  are  not  available. 

8.  Weather  data.  The  wsather  data  include  daily  precipitation 
amounts  (expressed  as  water  •qulvalent)  and  dates  of  the  beginning  of  seasons 
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when  the  rates  of  depletion  change.  Provision  is  made  to  allow  input  of 
air  temperature  and  snow  depth  for  future  use  with  the  Freeze-Thaw  model.1'* 
Output  data 

9.  Output  data  are  stored  in  permanent  files  for  use  by  other  per¬ 
formance  prediction  models,  for  printout,  or  for  input  to  plotting  programs. 

The  printing  format  is  suitable  for  use  in  reports.  Soil  moisture  content 
can  be  presented  separately  in  units  of  either  centimeters  or  inches  of  water 
per  soil  layer  (15  cm  thick)  or,  if  desired,  by  percentage  of  dry  weight.  Soil 
strength  can  be  written  as  01  and/or  RCI.  The  predictions  can  start  or 
end  at  any  time  during  the  year  as  directed,  and  can  run  a  part  of  a  year 
or  for  consecutive  years.  Predictions  can  be  made  for  one  or  more  selected 
years  of  extreme,  unique,  or  typical  rainfall  distributions.  In  these 
cases  yearly  starts  are  controlled  by  starting  moisture  content  values  Inserted 
at  the  beginning  of  the  year  of  weather  data  input.  When  predictions  for 
another  soil  are  needed,  a  restart  of  the  model  with  new  terrain  input 
information  is  required. 

.  Limitations  of  the  Model 


Layer 

10.  The  model  predicts  the  moisture  contents  and  strengths  of  the 
first  and  second  soil  layers  (0-15  and  15-30  cm  depths).  No  direct  account¬ 
ing  is  made  of  water  from  precipitation  that  is  absorbed  by  vegetation,  runs 
off  the  surface,  or  drains  to  lower  layers.  Depletion  is  considered  a  net 
loss,  and  losses  specifically  from  evaporation,  transpiration,  and  drainage 
are  not  differentiated. 
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Terrain 


—  _  j 

11.  The  model  docs  not  include  all  terrain  attributes  that 

can  affect  soil  moisture  content  and  strength.  It  does  not  consider  the  ! 

direct  influence  of  attributes  such,  as  water  tables,  frost,  snow,  and  soil 
tillage.  The  influence  of  these  attributes  has  been  studied,  but  relations 
have  not  been  developed  sufficiently  for  inclusion  in  the  present  model. 

It  should  be  noted,  however,  that  the  model  was  purposely  designed  in  a 
modular  format  to  allow  for  their  inclusion  at  a  later  time. 

12.  The  model  does  not  provide  daily  moisture  contents  and  strengths  for  ■ 

i 

surface  materials  having  strengths  that  normally  are  not  appreciably  affected 
by  moisture  fluctuations,  such  as  clean  gravels,  cobbles,  boulders,  stones, 
and  rock,  or  materials  in  which  dally  moisture  contents  seldom  vary  signif¬ 
icantly  such  as  very  poorly  drained  (water-logged)  organic  silts  and  clays, 
and  other  organics  (peat  and  muck).  A  general  estimate  of  water  content, 
density,  and  strength  (in  Cl)  for  these  materials  is  shown  in  table  2.*  Also, 
the  model  should  obviously  not  be  used  for  prediction  \n  irrigated  fields,  where 
the  moisture  content  or  strength  of  the  soils  is  not  related  to  precipitation. 

Weather 

13.  The  quality  of  the  weather  data,  like  that  of  the  terrain  data, 
can  have  an  important  bearing  on  the  accuracy  of  soil  moisture  prediction. 

The  prediction  of  moisture  contents  requires  input  of  daily  rainfall  amounts. 

Rainfall,  distributions  may  be  quite  variable  over  an  area,  and  an  intense  1 

rainfall  can  occur  less  than  a  kilometer  3way  from  a  light  shower.  This 


These  data  are  not  included  in  the  automated  program.  They  are  presented  for 
use  by  analysts  desiring  to  characterize  entire  surface  areas. 
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is  especially  true  in  mountainous  terrain,  where  sharp  changes  in  rainfall 
amount  commonly  occur  with  differences  in  prevailing  winds  or  in  elevation. 

In  some  of  the  mountainous  areas  of  Hawaii,  for  example,  the  average  annual 
rainfall  changes  from  50  to  500  cm  over  a  distance  of  10  km.  The  user  of 
the  model  should  recognize  this  limitation  and  try  to  obtain  rainfall  data 
chat  apply  to  the  area  of  prediction. 

14.  The  tentative  average  relations  used  with  estimated  data  (see 
paragraph  7b)  were  developed  from  specific  data  obtained  within  the  continental 

United  States.  Use  of  these  relations,  therefore,  should  be  limited  to 

l 

areas  of  similar  climate. 

Accuracy  of  predictions  for  individual  sites 

15.  Soil  moisture  content.  Prediction  accuracies  were  determined  from 

specific  relations  developed  from  23  test  sites.  The  average  deviation  of 

predicted  values  from  measured  values  for  the  15-  to  30-cm  layer  was 

about  +  1  percent  moisture  content,  dry  weight,  for  the  year  the  prediction 

4d 

relations  were  derived  and  about  +1.5  percent  the  following  year.  The 
natural  variability  of  moisture  content  within  a  site  at  a  given  time  is 
of  the  same  magnitude  (Appendix  D  of  reference  4e) .  From  average  relations 
and  from  rainfall  uv^a  collected  some  distance  from  the  sites,  the  average 
deviation  of  predicted  soil  moisture  contents  from  measured  values  for  the 
15-  to  30-cm  layer  (for  601  test  sites)  was  about  +  4  percent  moisture 
content.  Prediction  accuracy  for  the  0-  to  15-cm  layer  was  somewhat 
poorer  becaure  of  the  greater  variability  in  the  soil  at  the  surface. 

16.  Soil  strength.  Strengths  for  the  15-  to  30-cm  layer  were 
predicted  from  measured  moisture  contents  and  specific  soil  moisture- 
strength  relations.  The  average  deviations  of  these  strengths  from  measured 
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strengths  were  ahout  +  15  Cl  or  RCI  units  under  wet  soli  conditions;  the 

41 

accuracy  decreased  at  lower  moist urt  contents,  4  The  natural  variabilities 

4h 

of  Cl  and  RCI  within  a  site  are  of  about  the  same  magnitude.  Strengths 

for  the  same  layer  were  also  predicted  from  measured  moisture  contents  and 

estimated  soil  moisture-strength  relations.  The  average  deviations  of  these 

strengths  from  measured  strengths  were  +  30  units  for  Cl  and  +  20  units  for 

41 

RCI  under  wet  soil  conditions;  J  again,  the  accuracy  decreased  at  lower 
moisture  contents.  No  tests  were  made  for  prediction  accuracy  of  the  0-  to 
15-cm  layer;  but,  similar  to  moisture  content,  the  prediction  accuracy  is 
assumed  to  be  poorer  because  of  the  greater  soil  variability  at  the  surface. 
Accuracy  of  predictions  for  large  areas 

17.  No  tests  were  made  to  determine  the  accuracy  of  predictions  for 
large  areas  (areas  several  hundred  square  meters  or  more).  The  limitations 
of  time  and  money  required  for  the  development  of  the  computer  program  reported 
herein  precluded  the  design  and  operation  of  a  field  and  desk  study  to 
collect  and  analyze  the  necessary  data.  It  can  be  assumed,  however,  that 
the  range  of  variation  between  predicted  and  measured  values  would  be 
greater  than  for  individual  sites  because  of  the  greater  variability  of  terrain 
and  weather  properties. 
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PART  II;  STRUCTURE  OF  PROGRAM 


Computer  Orientation 

18.  The  program  for  this  model  is  written  in  Fortran  IV  language 
for  the  Honeyvell-GE  (General  Electric)  440*  computer  at  the  WES  and  is 
to  be  used  remotely  by  teletype.  The  system  compiler  Fortran  XFR  is 
used. 

19.  An  average  run  for  two  years  of  daily  precipitation  with  a  complete 
output  (soil  moisture  content  in  terms  of  volume  or  percent.  Cl  and  RCI  for 
the  two  layers)  takes  from  2  to  2-1/2  .min. 

20.  If  graphic  displays  are  required,  the  data  in  output  files  stored 
by  this  model  can  be  plotted  directly  by  using  the  proper  programs.  The 
graphic  display  program  is  discussed  in  paragraph  39. 


Main  Program  and  Overlays 

21 .  The  program  is  in  conversational  mode  and  is  of  modular  con¬ 
struction.  It  consists  of  a  main  program  and  six  overlays, 
and  it  accesses  up  to  four  input  files  and  creates  up  to  nine  output 
files.  If  generated  data  exceed  the  capacity  of  a  file,  the  program 
automatically  opens  another  file  and  names  it  in  sequential  order. 


The  Honeywell  435  was  upgraded  to  the  440  in  early  1972. 
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22.  The  main  program,  overlays,  and  required  input  control 
files  are  on  file  in  the  WES  Automatic  Data  Processing  Center  as  tabu¬ 


lated  below. 


File  or 


Characters 


Name 

Description  or  Function 

Source 

Ob.iect 

Lines 

FTWEMS 

Main  program 

17571 

15333 

472 

COMP MS 

Computation  core 

3234 

2628 

113 

FIL1DA 

Stores'  moisture  content  file 

4128 

3546 

127 

FIL2DA 

Stores  percent  moisture 
content  file 

3222 

2556 

97 

FRODEP 

Stores  output  for  freeze-thaw 
program 

2364 

1626 

62 

STFRSP 

Computes  and  stores  strength 
from  soil  properties 

2142 

2067 

70 

CALSST 

Calculates  soil  strength 

5991 

4875 

173 

DSURGR 

Soil  composition  group 
descriptors  -  data  file 

9906 

122 

DFSTEQ 

Coefficients  of  equations 
to  determine  strength 
relation  constants,  A  and 

B  -  data  file 

1932 

40 

Schematic  Flow  Chart 

23.  A  schematic  flow  chart  shows  the  logic  used  in  assembling  the 
program  (fig.  1).  Detailed  flow  charts  of  the  program  are  given  in 
Appendix  A, 
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OVERLAY  4 
^FRODEP" 

WRITES  DATA 
FILES  OF  DAILY: 
PRECIPITATION 
SNOW  DEPTH 
AIR  TEMPERATURE 
PERCENT  SOIL 
MOISTURE  FOR 
EACH  SOIL 
LAYER 


OVERLAY  5 
"STFRSP" 

CALCULATES 
STRENGTH 
EQUATION 
CONSTANTS 
FROM  SOIL 
PROPERTIES 


OVERLAY  6 
"CALSST" 

CALCULATES 
CONE  INDEX 
AND/OR  RATING 
CONE  INDEX 
FOR  EACH 
SOIL  LAYER 


<D 


Fig.  1.  (Concluded) 


Function  of  Main  Program  and  Overlays 

24.  The  following  outline  lists  the  action  accomplished  in  each  segment 
of  the  program  and  the  order  in  which  each,  if  required,  is  called  into  use: 

a.  Main  program  -  FTWEMS  -  a  program  JF  developed  in  the 

Terrain  Analysis  Branch  of  the  Mobility  and  Environmental 
Systems  Laboratory  of  the  WES  for  predicting  soil  Moisture 
content  and  soil  Strength. 

(1)  Defines  dimensions  of  the  array  of  each  variable  con¬ 
tained  in  "common11  statements. 

(2)  Identifies  initial  data.  - 

(3)  Identifies  formats  of  input  data. 

(4)  Reads  all  input  control  data. 

(5)  Reads  all  weather  data. 

(6)  Converts  input  as  required. 

(7)  Activates  the  proper  overlay  at  the  proper  time. 

(8)  Contains  a  subroutine  called  by  overlays  for 
special  formatting  of  output  when  precipitation 
record  begins  other  than  the  first  day  of  a  week. 

_b.  Overlay  1  -  COMPSM  -  COMP u tea  Soil  Moisture  content. 

(1)  Computes  daily  soil  moisture  contents  for  the  top 
two  soil  layers. 

(2)  Stores  daily  soil  moisture  content  data  in  data 
blocks  of  one  year  (January-December) . 

£.  Overlay  2  -  FIL1DA  -  creates  FILes  for  1st  type  of  output 

DAta.  soil  moisture  contents  by  unit  layer.  Stores 
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soil  moisture  content  data  in  a  separate  file 
for  each  layer.  If  required,  converts  soil 
moisture  contents  to  percent  dry  weight. 


<1.  Overlay  3  -  FIL2DA  -  creates  FILes  for  2nd  type  of  output 
DAta.  soil  moisture  content  by  percent.  Stores  data  in  a 
separate  file  for  each  layer. 

<e.  Overlay  4  -  FRODEP  -  creates  a  file  of  data  for  use  in  models 
that  predict  FROzen  PEPth  or  thawed  depth  of  soil.  Stores 
dally  precipitation,  daily  air  temperature,  measured  daily 
snow  depth  (if  available),  and  dally  percentage  of  soil 
moisture  content  data  by  layers  in  an  output  file.* 
f..  Overlay  5  -  STPRSP  -  calculates  soil  STrength  equation  con¬ 
stants  FRom  soil  moisture  contents  (at  specific  values  of  Cl 
and  RCI)  that  have  been  derived  from  Soil  Property  values. 
These  constants  are  stored  for  use  in  Overlay  6. 

Overlay  6  -  CALSST  -  CALculates  .Soil  STrength.  Calculates 
daily  Cl  and/or  RCI. 

Listings  of  these  program  files  are  in  Appendix  B. 

Input  and  Output  Data 

25.  The  computer  operator  can  enter  the  input  data,  described  in 
paragraphs  5-8»  into  a  designated  input  file  either  by  paper  tape  or 
directly  from  the  teletype  keyboard,  or  he  can  punch  them  on  80-column 


*This  overlay  is  not  required  for  operation  of  the  program.  It  is  included 
for  the  convenience  of  the  user. 
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c<_rds  and  enter  them  in  a  time-sharing  file.  Some  control  data  are 
entered  from  the  teletype  keyboard  only.  This  input  is  concerned  with 
input  and  output  options  in  answer  to  questions  posed  during  the  running 
of  the  program. 

26.  Output  data  generated  by  the  program  are  stored  in  time¬ 
sharing  files,  which  are  formatted  for  printing  by  SYS:MAX  &  WES  media  con¬ 
version  system)  without  line  numbers.  A  built-in  routine  slews  the 
printer  to  the  top  of  the  next  page  when  a  page  of  data  is  complete. 

The  output  files  remain  in  the  time-sharing  storage  until  the  user  unsaves 
them. 

Definitions 

27.  Table  3  contains  a  list  of  device,  variable,  constant,  and 
flag  definitions,  and  units  of  measurements  with  significant  decimal 
places  where  applicable. 
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PART  III:  PROGRAM  OPERATION 


28,  The  successful  operation  of  the  program  depends  on  careful 
preparation  of  the  input  data,  an  understanding  of  the  user  function  in 
operating  the  program,  knowledge  of  the  limitations  and  conventions 
adopted  by  the  program,  and  the  form,  restrictions,  and  location  of  the 
program  output  data. 


Input  Data  Preparation 

29.  Weather  data  preparation  i3  the  same  regardless  of  the  input 
data  source.  All  other  data  are  control  data  (see  paragraphs  6-1).  in 
this  discussion  the  terms  card  and  line  are  used  synonymously  since  data 
from  a  card  fills  a  teletype  data  line. 

Control  data 

30.  Control  data  describe  the  top  two  soil  layers  and  appear  on 
seven  different  card  types  (1-7)  as  shown  in  figs.  C6-C17  of  Appendix  C. 

31.  Specific  data.  If  the  user  wishes  to  run  a  prediction  based 
entirely  on  specific  data  he  must  input  on  the  proper  cards  all  of  the 
control  data  (see  paragraph  6).  The  card  deck  setup  is  shown  in  fig.  C3, 
Appendix  C. 

32.  Surface  composition  group  data.  A  user  selecting  control  data 
input  from  the  surface  composition  groups  needs  only  to  use  the  five 
heading  cards  (type  7)  and  the  cards  with  weather  data  (types  8  and  9)  to 


make  a  prediction  run.  If  he  has  information  concerning  beginning  mois¬ 
ture  content  for  each  layer,  he  should  add  card  type  3.  He  can  also  add 
any  other  card  type  for  which  he  has  specific  data.  The  cards,  including 
card  type  3,  must  precede  cards  of  type  7.  The  card  deck  setup  is  shown 
in  fig.  C4,  Appendix  C.  The  control  data  card  input  for  each  surface  compo¬ 
sition  group  is  stored  in  file  DSURGR.  The  file  is  discussed  in  Appendix  D. 

33.  Estimated  data.  Tentative  average  control  data  stored  in  the 
main  program  are: 

a.  Accretion  relation  constants. 

t>.  Depletion  relation  constants. 

c.  Minimum  precipitation  causing  accretion. 

d_.  Moisture  ranges  of  depletion  by  soil  type,  season,  and 

layer. 

£.  Maximum  and  minimum  soil  moisture  equation  constants. 

Table  4  summarizes  the  relations  from  which  these  constants  are  taken. 

34.  Other  estimated  data  that  are  stored  are  soil  property  equa¬ 
tion  constants  for  calculating  strength  relation  coefficients.  These 
constants  with  codes  that  indicate  the  soil  property  and  its  form  (In, 
lnln,  or  natural)  are  stored  in  file  DFSTEQ  and  are  listed  in  table  D2, 
Appendix  D. 

Weather  data 

35.  Weather  data  appear  on  two  different  card  types  (8  and  9)  as 
illustrated  in  figs.  C18  and  C19  of  Appendix  C. 
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User  Function 


36.  The  user  is  assuaed  to  be  familiar  with  the  remote  operation 
of  the  computer  by  teletype  and  to  have  checked  the  compatibility  of 
the  computer  system  he  is  using  with  the  system  for  which  this  program 
was  written.  He  can  submit  program  card  decks  to  the  computer  operator 
for  loading  the  model  into  the  time-sharing  disc  pack,  or  he  can  enter 
them  by  paper  tape  through  the  teletype.  The  main  program,  all  over¬ 
lays,  and  rll  input  files  must  be  present  to  run  the  program. 

37.  The  user  must  prepare  control  data  and  weather  data  as  described 
in  paragraphs  29-35  and  in  Appendix  C,  and  he  must  be  prepared  to  enter 
file  names  descriptive  of  the  type  of  output  data  desired.  A  file  name 

is  a  six-character  alphanumeric,  the  last  three  characters  of  which  must 
be  001  as  requested  during  the  running  of  the  program.  This  makes  pos¬ 
sible  the  automatic  incrementation  of  output  data  files  in  sequential 
order — XXX001,  XXX002,  etc. — if  necessary. 

38.  Three  examples  of  the  program  operation  using  specific,  surface 
composition  group,  and  estimated  data  are  shown  in  Appendix  E.  Example  1 
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uses  only  specific  data.  Example  2  uses  Che  same  weather  data,  but  a 
surface  composition  group  as  a  control  data  source.  Example  3  uses 
estimated  data  with  tentative  average  relations  for  the  prediction. 

These  examples  also  illustrate  the  user's  response  to  questions  asked 
during  the  running  of  the  program.  A  run  providing  data  for  the  Freeze- 
Thaw  model  is  included  as  Example  4  in  Appendix  E. 

39.  If  plotting  of  the  output  files  is  desired,  a  separate  program, 
FSPMSD,  is  included  for  the  user's  convenience  (Appendix  F) .  If  this 
program  is  used,  a  CALCOMP  drum  (pen)  plotter  and  required  operating 
library  routines  must  be  available.  If  the  required  library  routines 
listed  at  the  beginning  of  the  program  are  not  available  in  his  computer 
system,  the  user  must  write  a  plotting  program  that  will  be  accepted  by 
the  system. 


Program  Conventions 

40.  Program  conventions  that  must  be  understood  by  the  user,  in 
addition  to  those  observed  in  data  preparation,  are  as  follows: 

e.  The  internal  operation  of  the  moisture  prediction  phase  of 
the  program  is  in  inches  of  moisture  per  6-in.  soil  layer. 

b. .  All  data  input  by  the  user,  except  density,  must  be  in  the 

^.ame  system  of  measurement  (either  all  inches  or  all 
centimeters) . 

c. .  Density  is  always  entered  as  g/cc. 


All  72  spaces  on  a  teletype  line  must  be  filled  by  data  or 
accounted  for  by  spaces. 

The  soil  strength  and  freeze-thaw  programs  use  relations 
in  terms  of  percent  soil  moisture  content;  therefore,  the 
percent  option  must  be  selected  for  the  running  of  these 
programs . 

Entry  of  minimum  storm  data  must  be  in  three  decimal  places 
(0.100  in.  or  0.254  cm). 

In  converting  centimeters  to  inches,  the  constant  0.3937 
is  used  as  a  multiplier;  and  in  converting  inches  to  centi¬ 
meters,  the  constant  2.54  is  used. 

In  converting  output  from  centimeters  to  percent,  15.24  cm 
is  used  for  the  thickness  of  the  6-in.  increment. 
Identification  of  the  run  on  the  top  of  each  page  that  is 
printed  from  an  output  file  is  made  from  the  identification 
information  on  card  type  1,  and  not  from  the  card  type  7 
heading  cards  (see  Appendix  C). 

The  last  four  characters  of  the  identification  must  be 
numeric  and  must  agree  with  the  surface  composition  group, 
if  a  group  is  used  as  an  input  data  source. 

If  air  temperature  data  are  included,  both  temperature  and 
snow  data  cards  must  follow  the  last  precipitation  data 
card  for  each  year. 

The  program  requires  a  "RUNBIG"  call  for  its  initiation. 
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m.  This  program  uses  overlays;  thus,  object  files  with  the 
following  names  are  held  under  the  user  work  number  for 
use  in  running  the  program: 


(1) 

FT  WE  MS 

-  0(zero)FTWEM 

(2) 

COMPMS 

-  Q(zero)COMPM 

(3) 

FIL1DA 

0(zero)FILlD 

(4) 

FIL2DA 

-  0(zero)FIL2D 

(5) 

CALSST 

0(zero)CALSS 

(6) 

FRODEP 

0(zero)FR0DE 

(7) 

STFRSF 

-  0(zero)STFR3 

To  save  compilation  time,  the  user  can  call  the  program  "OFTWEM" 
and  initiate  it  with  a  call  to  "RUNBIG. "  This  results  in  a  sav¬ 
ing  of  approximately  30  sec  of  running  time. 

r».  If  the  predicted  moisture  content  on  the  last  day  of  a 
season  is  lower  than  the  minimum  moisture  content  on  the 
first  day  of  the  following  season,  daily  depletion  is  not 
computed  until  soil  moisture  content  has  accreted  above 
the  new  season  minimum. 

£.  If  precipitation  cards  are  out  of  order,  an  error  message 
is  printed  on  the  teletype.  If  a  precipitation  card  is 
omitted,  no  error  is  printed,  but  zeros  are  written  into 
the  output  file  for  the  days  omitted.  The  output  data 
beyond  that  point  will  be  in  error. 

£.  The  program  rounds  values  as  follows: 
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(1)  Soil  moisture  contents  in  inches  or  centimeters  are 
rounded  to  the  nearest  hundredth  of  the  measurement 
unit . 

(2)  Soil  moisture  content  is  rounded  to  the  nearest  hun¬ 
dredth  of  a  percent. 

(3)  Strength  data  are  rounded  to  the  nearest  unit  (Cl  or 
RCI). 

Some  check  data  are  printed  by  the  teletype  during  the 
running  of  the  program.  Should  these  check  data  disagree 
with  those  desired,  the  program  can  be  aborted  by  hitting 
the  stop  button  (S) . 


PART  IV:  MATHEMATICS 


41.  The  mathematics  used  by  the  program  are  simple  straight  line 
and  power  function  relations  derived  empirically  using  regression  tech¬ 
niques.  The  overlays  "COMPSM"  and  "CALSST"  contain  instructions  for  the 
calculation  of  soil  moisture  content  and  soxl  strength  using  the  control 
data  as  defined  in  paragraph  6.  The  program  was  written  primarily  for  use 
with  TA  (estimated)  input  data  with  options  to  substitute  data  supplied 
trom  specific  sources  (paragraph  7a)  or  data  characteristic  of  a  surface 
composition  group  (paragraph  7c) .  This  main  core  of  computations  makes 
possible  the  total  use,  partial  use,  or  the  use  of  a  combination  of  any  of 
the  three  types  of  input  data  (see  Appendix  C).  Provision  is  also  made 
for  the  calculations  of  maximum  and  minimum  soil  moisture  contents  and  of 
the  soil  moisture-strength  relations  if  desired. 

Soil  Moisture  Prediction 


Accretion  relations 

42.  Accretion  (wetting)  of  soil  depends  primarily  upon  the  amount  of 

precipitation  and  amount  of  storage  space  available  in  the  soil  for  absorbing 
Ac  Ad 

water.  *  Precipitation  is  the  critical  variable  when  the  total  precipita¬ 
tion  for  a  storm  is  less  than  the  storage  available  (in  the  top  two  layers 
of  soil).  This  type  of  accretion  is  called  "Class  I."  "Class  II"  accretion 
occurs  when  the  precipitation  is  greater  than  the  available  storage;  then 
accretion  depends  on  the  amount  of  available  storage.  Equations  expressing 
these  relations  are  shown  in  table  4. 
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43.  "RMIN"  (minimum  storm)  defines  the  amount  of  precipitation 
below  which  no  appreciable  wetting  of  the  soil  occurs.  Unless  the  user 
enters  a  value  of  RMIN  on  card  2-1  (see  fig.  C7,  Appendix  C),  the  program 
assumes  an  RMIN  of  0.1  in.  (0.254  cm).  If  precipitation  exactly  equals 
RMIN,  no  accretion  or  depletion  occurs.  For  precipitation  less  than  RMIN, 
depletion  is  determined  in  the  usual  manner. 

44.  Minimum  values  of  accretion  are  set  at  0.03  in.  for  Class  II 
accretion  for  both  layers  and  0.01  in.  for  layer  2-Class  I  accretion. 

Layer  1-Class  I  accretion  for  precipitation  between  RMIN  and  RMIN  +  0.07  in. 
(i.e.  between  0.1  in.  and  0.17  in.  in  fig.  2)  has  a  1:1  relation  to  pre¬ 
cipitation  above  RMIN  when  the  minimum  storm  is  0.1  in.  If  RMIN  is  entered 
as  some  value  other  than  0.10  in.,  the  accretion  value  0.07  in.  must  be 
recomputed.  For  example,  as  shown  in  fig.  2,  a  minimum  storm  (RMIN)  of 
0.06  in.  would  result  in  a  minimum  accretion  of  0.034  in.  The  user  must 
enter  these  minimum  values  on  card  2-1  (see  fig.  C7,  Appendix  C)  if  new 
accretion  equation  constants  are  used.  Accretion  relations  can  be  modi¬ 
fied  by  entering  equation  constants  on  card  2-6  (see  fig.  C12,  Appendix  C). 
The  same  equation  form  must  be  used. 

Depletion  relations 

45.  The  program  accepts  depletion  relation  equations  in  polynomial 
form  up  to  sixth-degree  equations  forced  through  zero.  The  surface  compo¬ 
sition  groups  use  fourth-degree  equations.  The  equation  coefficients  are 

in  file  DSURGR  (see  table  Dl,  Appendix  D).  TA  depletion  equations  (estimated 
data)  are  fourth-degree  polynomials,  except  for  two  that  are  sixth-degree 
polynomials.  The  equations  define  depletion  curves  that  were  determined 
empirically. 
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Fig.  2.  First-layer  Class  I  accretion 


U6.  The  equation  coefficients  were  determined  by  a  multiple  linear 
regression  technique  with  the  intercept  equal  to  zero.  Daily  moisture 
loss  was  considered  dependent  on  soil  moisture  content  above  minimum 
moisture  content.  The  equation  form  Is: 

Y  *  ax  +  bx  +  cx  +  dx  +  ex  +  fx 

47.  For  a  specific  location,  if  the  maximum  and  minimum  soil  moisture 
contents  for  each  layer  can  be  determined,  these  TA  depletion  equations  can 
be  used.  If  this  is  the  case,  the  program  computes  a  factor  by  which  the 

site  is  referenced  to  the  TA  moisture  depletion  daring  summer,  (site  maximum  i 

moisture  -  minimum  moisture)/(TA  maximum  moisture  -  minimum  moisture).  This, 
in  effect,  stretches  or  compresses  the  TA  depletion  curves  to  force  them 
into  agreement  with  the  specific  location  depletion  relation. 

48.  There  are  18  equations  in  the  set,  three  for  each  soil  layer 
for  each  of  three  seasons;  they  are  listed  in  table  4. 

Maximum  and  minimum  soil  moisture  contents 

49.  TA  maximum  and  minimum  soil  moisture  equations  (estimated  data) 
are  empirically  derived  relations.  The  constant  terms  were  derived  by 
regression  analysis.  Maximum  and  minimum  moisture  contents  were  con¬ 
sidered  dependent  on  soil  properties.  The  properties  included  in 

the  program  are  improved  revisions  of  these  equations  and  can  be  used  at 
the  user's  option  to  compute  maximum  and  minimum  moisture  contents  if  the 
soil  properties  for  each  soil  layer  are  available.  The  equations  are 
shown  in  table  4  and  are  of  the  form: 
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Y  =  a  +  bCSPj^)  +  c(SP2)  + . n^sp16> 

where 

Y  ■  maximum  or  minimum  moisture  content 
SP  =  soil  property 

50.  The  user  can  supply  his  own  equations,  if  they  are  of  the  same 
form,  by  entering  soil  properties  and  equation  constants  on  card  type  2-2 
as  shown  in  fig.  C8,  Appendix  C. 

Soil  Strength  Prediction 

51.  Both  the  Cl  and  the  RCI  relations  to  soil  moisture  content  were 

determined  empirically  to  be  straight  line  logarithmic.  Natural  logarithms 

are  used  with  specific,  estimated,  or  surface  composition  group  data  in 

this  program.  Equations  have  the  following  form: 

In  Cl  «  a  +  bln  HC 
In  RCI  *  a  +  bln  MC 

where 

Cl  =  cone  index 
RCI  =  rating  cone  index 
MC  =  soil  moisture  content  in  percent 

Given  relations 

52.  If  surface  composition  groups  are  used  for  prediction,  each  group 
has  an  associated  Cl  and  RCI  equation  for  each  soil  layer.  If  the  user 
does  not  use  a  group,  he  must  supply  the  strength  equation  constants  using 
specific  or  estimated  data.  Information  for  entering  these  constants  can 
be  found  on  card  type  2-4,  as  shown  in  fig.  CIO,  Appendix  C. 
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Relations  computed  from  soil 
properties  (estimated  data) 

53.  In  a  study  analyzing  the  relation  of  soil  strength  to  soil 

4  \ 

properties,  J  six  sets  of  equations  were  presented  to  predict  moisture 
content  at  specific  Cl's  and  six  sets  to  predict  moisture  content  at 
specific  RCl's.  A  multiple  linear  regression  technique  was  used  to 
determine  soil  properties  most  closely  associated  with  moisture  con¬ 
tents  at  Cl  values  of  200  and  300  and  RCI  values  of  100  and  200.  Nine 
soil  properties  were  examined,  and  only  those  that  made  a  significant 
statistical  improvement  in  the  relation  were  actually  used.  Tables  6 
and  7  show  the  various  resulting  combinations  which  were  retained  in 
the  equation  sets.  The  soil  properties  used  are  as  follows: 
a..  Percent  sand 
In  Percent  clay 
£.  Percent  silt 
d_.  Percent  organic  matter 
£ .  Percent  f ines 
f_.  Liquid  limit 
£.  Plastic  limit 
h.  Plasticity  index 
jL  Dry  density 

54.  Once  these  relations  were  determined  for  moisture  contents  at  the 
strengths  stated  in  paragraph  53,  a  straight-line  relation  was  assumed  to 
exist  between  In  Cl  and  In  soil  moisture  content  (see  paragraph  51).  The 
equation  for  this  straight  line  was  then  determined  by  forcing  it  through 
soil  moisture  contents  at  a  Cl  of  300  and  a  Cl  of  200  (fig.  3).  The  same 


28 


Fig.  3.  Cone  index  -  percent  soil  moisture  content  relation 


sets  of  equations  are  considered  applicable  to  both  layers.  Similarly 
straight-line  relations  were  established  for  In  RCI  versus  In  soil 
moisture  content. 

55.  The  partial  regression  coefficients  for  the  six  sets  of  equa- 
tions  used  in  calculating  Cl  coefficients  are  listed  in  table  6.  A  se: 
consists  of  equations  for  calculating  moisture  contents  when  Cl  is  200 
and  when  Cl  is  300.  The  soil  property  coeificients  are  listed  in  the 
order  in  which  they  are  entered  into  the  equation. 

56.  Partial  regression  coefficients  for  equations  used  in  calcu¬ 
lating  moisture  contents  when  RCI  is  100  and  200  are  listed  in  tabic  7 
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Table  1 

Surface  Composition  Groups 


Material 


irganic  Drainage* 
lontent  Potential  G 
Z  Class  C 


Groups  with  Similar  Material  in  0-  to  15-  and  15-  to  30-cm  Layers 


Pavement  and  structures  >25Z  coverage 

Rock,  stones,  boulders,  and  cobbles> 

P.D.  sizes  >6. 07 4am, is  >50Z,  and  >76.2ma, is  >25Z 


Coarse  grained  Gravel,  P.D.  Clean  gravel,  P.D. ' 
P.D.  sizes  sizes  >4.76-  sizes  <0. 074mm. is  <5Z 

>0.074ma, is  76.2mm, is  Gravel  with  fines,  P.D. 

_>50Z  >25Z _ isizes  <0.074mm.is  >5-50Z 

Sand,  P.D.  Clean  sand,  P.D.  sizes 

sizes  >0.074-  <0.074na,is  <5Z _ 

4.76mm, is  Sand  with  fines,  P.D. 

>25Z  sizes  <0. 074mm, is 


Fine  grained,  |Silt,  LL  <35  and  PI  <15 
sizes<0.074n«B, 

is  >50Z  I  Clay,  LL  >35  or  PI  >15 


Organic  silts  and  clays  (plastic)  >7-30  0 _ j  1212 

_ _ _ _ * - i - [1313- 

Peat  (nonplastic) _ ....>30 _ 0 _  [1414  .  . 

Groups  with  Different  Material  in  0-  to  15- and  15-  to  30-cm  Layers  _ 

Sand,  0-15  cm,  over  I  1 _ 0610 

Clay,  15-30  cm _ _ .  ...  |  >Q_7  2  0711  _ 

Silt, "0- 15  cm,  over  !  3  nq,. 

Clay,  15-30  cm  _  _ _ 

^Drainage  potential  classified  by  occurrence  of  water  table  as  follows: 

Class  0  Water  table  occurs  at  surface  90Z  or  more  of  the  time 

Class  1  Water  table  occurs  at  the  surface  less  than  90Z  and  above  120-cm  depth 
10Z  or  more  of  the  time 

Class  2  Water  table  occurs  above  120-cm  depth  less  than  10Z  of  the  time 

**  Gravel  with  sand  matrix  coded  0707;  gravel  with  clay  matrix  coded  1111 
NOTE:  PD  •  Particle  diameter 
LL  *  Liquid  limit 
PI  «  Plasticity  index 
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Table  2 


Soil  Water,  Density,  and  Strength  for 
Surface  Composition  Groups  with 
Constant  Values 


Material 


Group 

Code 


Water 


{  Pavements  &  Structures  0101 

i 

f  Rock,  Stone,  Cobbles,  Boufdeis  0202 


Water *  **  j 

Content  Density  Strength 
.  % . g/cc  pi  ...J 

100  1.00  0  ! 

(liquid)j 
1  2.50  750+  ! 

I 

» 

1  to  5  2.15  750+  ! 


Clean  Gravel  0303  1  2.00  100 


j  Saturated  Organic  Silt-Clay  1212  90f  0.80  25 

f-  ......... 

J  Peat  &  Muck  1414  90t  0.80  25 


*  Percent  on  dry  weight  basis  except  for  water 

**  Dry  density  except  for  water 

+  Represents  an  average  value  estimated  from  a  small  number  of  samples. 
Water  contents  are  highly  variable  and  increase  with  an  increase  in 
the  percent  organic  matter  of  the  material. 
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Table  3 


BD(2)  Bulk  density  (dry  weight)  g/cc 

Cl (l6,2 )  Individual  soil  properties  variable 

C2(4)  Maximum  and  minimum  in./6  in.  or 

moisture  content  cm/15  cm 

IDC1(35)  Storage  location  for  rereading 

input  data 

I HD (20)  Storage  for  one  line  of  heading 

IMCA(4,12)  Names  of  months 

ISINUM(4)  Four  three-character  words  making 

up  a  site  identification  number 
ISL£W(3)  Code  for  slewing  to  top  of  next  page  - 

when  printing  output  by  the  WES  SYS:MAX 


media  conversion  system 


2 

3 

2 
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Table  3  (continued) 


Name 

KIND(3) 


MAXD(12) 

NMAXD(2) 


Mi£CIP(12,35) 


ADF(2) 

C3(2,2) 

C4(2,2) 

C5(l6,4,2) 

CMAX(2) 

CMIN(2) 

C0EF(108) 


Decimal 

Description _ Units  Places 

Code  for  tv  e  of  surface  material  and 
layers  one  and  two  as  needed  by  Freeze-Thaw 
model 

1  -  coarse 

2  -  fine 

3  -  other 

Number  of  days  in  each  month 

(1)  Number  of  last  month  having  data  -  •  - 

for  the  year 

(2)  Day  of  month  for  last  data  for 
the  year 

Storage  of  daily  precipitation  for  in.  2 

one  year  by  month  ami  day  or  cm 


Moisture  Content  Prediction 
Inverse  depletion  factor 
Accretion  equation  coefficients 
Accretion  equation  constants 
Maximum  and  minimum  equation  constants 
Maximum  soil  moisture  content 

Minimum  soil  moisture  content  modified 
by  the  depletion  factor 
Tentative  average  depletion  equation 
coefficients 


4 

4 


in./6  in.  2 
or  cm/15  cm 


» 
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Table  3  (continued) 


Name _  Description _ 

CONEW(3,2,6)  Depletion  equation  coefficient:!  for 

specific  relations  (polynomial  -  up 
to  sixth  degree) 

CXMIN(2)  Minimum  soil  moisture  content 


DWMCt'2,7) 

PMC(2) 

RANG(2 ,3,3) 

SDF(2) 

STGRMC(12,31) 

STORMD(12,31) 

X(2) 

XRANG(2,3) 


Soil  moisture  content  in  percent 
dry  weight 

Present  (prior  to  prediction)  soil 
moisture  content 

Maximum  and  minimum  soil  moisture 
content  limits  of  the  tentative  average 
depletion  curves 
Site  depletion  factor 

First  layer  soil  moisture  content  storage 
for  one  year 

Second  layer  soil  moisture  content 
storage  for  one  year 
Available  storage  for  soil  moisture 
content  (maximum  moisture  content  less 
present  moisture  content) 

Soil  moisture  content  limits  (maximum 
and  minimum)  for  three  climatic  seasons 
for  specific  depletion  equations 


Decimal 
Units  Places 

7 

Same  as  2 

CMAX 

Z  2 

Same  as  2 

CMAX 

Same  as  2 

CMAX 

Same  as  2 

CMAX 

Same  as  2 

CMAX 

Same  as  2 

CMAX 

Same  as  2 

CMAX 
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'able  3  (continued) 


Name 


ABAR(2) 
B(13,2 ,2) 

BBAR(2,2) 

C0N1(2,2) 

D(2),E(2) 

IMOD(12,3) 

STC0N(2,4) 


STC(12,31) 

STRLI(2,4) 


SNC(12,35) 
TEMPER (12, 35) 


Decimal 

Description _ Units  Places 

Strength  Prediction 

Calculated  equation  constants  for  cone  index 
Equation  coefficients  for  moisture-strength 
from  soil  properties 

Calculated  equation  constants  for  rating 
cone  index 

Natural  logarithms  of  200,  100,  300,  200 
Intermediate  summations  in  equations 
determining  strength  equation  constants 
Form  of  soil  property  as  it  appears  in  - 

the  equation 

Soil  strength  equation  constants  -  3 

(1.4)  cone  index 

(2.4)  rating  cone  index 

Daily  soil  strengths  for  one  year  by  months  -  0 

Maximum  and  minimum  strength  limits  -  0 

(1.4)  cone  index 

(2.4)  rating  cone  index 

Freeze-Thaw  Model  Input 

Daily  snow  depths  for  one  year  by  months  m  2 

Daily  temperatures  for  one  year  by  months  degC  2 
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Table  3  (continued) 


Name 

Other  Variables,  Constants  and  Flags 

Description 

Units 

Decimal 

Places 

AJBN,JBN,JN 

Input-Output 

Beginning  day  number  divided  by  neven 

. 

DMCLAl , DMCLA2 

Variables  holding  unit  volume  moisture 

- 

- 

DPMLA1.DPMLA2 

content  output  file  names 

Variables  holding  percent  moisture  content 

_ 

— 

DCILA1,DCILA2 

output  file  names 

Variables  holding  cone  index  output  file 

— 

— 

DRCLA1 , DRCLA2 

names 

Variables  holding  rating  cone  index 

_ 

— 

DFRDEP 

output  file  names 

Variables  holding  Freeze-Thaw  model 

_ 

_ 

IBD.JB 

input  name 

Beginning  day  of  data 

_ 

— 

IBM,K 

Beginning  month  of  data 

- 

- 

IBY 

Beginning  year  of  data 

- 

- 

IDAB 

Beginning  day  of  week 

- 

- 

I  DAE 

Ending  day  of  week 

- 

- 

IND 

Computed  format  no. 

- 

- 

ITEM1 , ITEM2 , IDEV1 

Additional  reread  storage 

- 

- 

IIT 

Nisaber  of  soil  properties  to  be  read 

- 

- 

ISD 

Day  of  new  season 

- 

- 
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Table  3  (continued) 


ISM 

Month  of  new  season 

1S1 

Code  for  climatic  season 

1  Summer 

2  Transition  (fall  or  spring) 

3  Winter 

ISV 

Year  of  new  season 

ITYPE 

Code  indicating  data  type 

IUNMEA- 

Code  for  input  unit  of  measurement 

JEND 

Last  day  of  data  for  year 

JNK 

Index  for  weeks,  1-5 

LK 

Last  month  of  data  for  year 

MEAOU1 

Code  for  output  unit  of  measurement 

MEAOU2 

Code  for  output  in  percent  moisture  content 

newyr.nyr 

New  year  of  data 

NSOUR 

Code  for  input  source 

NTYPE 

Code  for  subtype  data  identification 

SECT1. . . • 

SECT10  Sector  locators 

Moisture  Content  Prediction 

ACCR 

Amount  of  day’8  accretion 

in. 

C1CGF 

0.07  (point  at  which  accretion  equation 

is  modified) 

in. 

C1ROF 

Amount  precipitation  exceeds  minimum  rain 
when  accretion  is  less  than  0.07  in. 

in. 
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Table  3  (concluded) 


Naive 


ACCR2 

DEP 

IDA,MC,1YR 
I COUNT 

1TEXA.ITEXB 


RMIN 


NCI 

NEWLIM 

NRCI 

NSOST 


NCODE 

NSKOW 

NTEMP 


Decimal 

Descript  ion _  Units  Places 

0.03  (point  at  which  accretion  equation  is  in.  2 

modified) 

Amount  of  day's  depletion  in.  2 

Present  day,  month,  year 

Number  of  days  in  current  manth  for 

moisture  calculation 

Soil  texture  for  layer  1  and  layer  2  - 

C  coarse  (sands) 

M  medium  (silts) 

F  fine  (clays) 

Minimum  precipitation  causing  accretion 
Strength  Prediction 

Number  identifying  cone  index  equation 
desired 

Code  for  new  soil  strength  limits 

Number  identifying  rating  cone  index  equation  - 

Code  for  strength  equation  constants 


Freeze-Thuw  Model  Input 
Code  for  air  temperature  and/or  snow 
depth  data 

Code  for  sncw  depth  data 
Code  for  air  temperature  data 
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Table  4 


Built-In  Program  Control  0| 
ACCRETION  RKl.ATTnvo  ISStatlve  Average  RelaMnn., 


Class  I 
Class  II 

where 


0~  to  15-ca  Layer 

Y  *  0.47X  -  0.01 

Y  -  0.752  -  0.05 


to  30—  ca  Laver 

Y  -  0.22X  -  0.01 

Y  -  0.60Z  -  0.02 


A 

Y  -  predicted  accretion,  m.  per  6-in.  layer 
X  -  rainfall,  in. 


arauBDie  storage 


~ -  wmc  OLdl L 


depletion  relations 

Sumner: 

Sand  - 

First  layer  moisture  range,  1.26 

Hx)  «  0. 15970948x  -  0.28827749x2  +  0.37098284*3 

Second  layer  moisture  range,  1.00 

f(x)  -  0. 20865303*  -  0.64984646x2  +  0.87560732*3 

Silt  - 


0.13067845X* 

0.30605206x4 


First  layer  moisture  range,  1.87 
f(x)  -  0.15190803x  -  0.077 20229x2  - 
Second  layer  moisture  range,  1.63 
f(x)  -  0.09612290*  +  0.01073670x2  - 
Clay  - 


0.01280603x3  +  0.01851488x4 
0.10720756*3  +  0.05548077x* 


First  layer  moisture  range,  I.47 
f(x)  *  0.11641490*  +  0.06704479*2  _ 
Second  layer  moisture  range,  1.08 
f(x)  »  0.14304453x  -  0.10102288x2  - 


0.09112914x3  +  0.03120642x4 
0.02103815x3  +  0.10071059x4 
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Table  4  (continued) 


DEPLETION  RELATIONS  (con.) 

Transition: 

Sand  - 

First  layer  moisture  range,  1.05 

f(x)  -  0.08401220x  -  0.01562037x2  -  0.10949078x3  +  0.12228349x4 
Second  layer  moisture  range,  0.72 

f(x)  =  0.12713529X  -  0.53982006x2  +  1.12524850x3  -  0.58561447x4 
Silt  - 

First  layer  moisture  range,  1.55 

f (x)  -  0.1138717x  -  0.6558373X2  +  1.8700644x3  -  2.3253802x4  + 
1.2952910x5  -  0.26246686x6 
Second  layer  moisture  range,  1.25 

£(x)  -  0.02789788x  +  0.24453131x2  -  0.86060142x3  +  1.2998033x4  - 
0.097359290X5  +  0.29132112x® 


Clay  - 

First  layer  moisture  range,  1.43 

f  (x)  °  0.083218ilx  •<*  0.11716883X2  -  0.14710818x3  +  0.04910032X4 
Second  layer  moisture  range,  0.93 

f(x)  =  0.03492386x  +  0.36829518X2  -  0.90742210x3  +  0.61582018x4 

Winter: 

Sand  - 

First  layer  moisture  range,  0.54 

f (x)  -  0.12882001x  -  0.32751050x2  +  1.0509886x3  -  1.0462256x4 
Second  layer  moisture  range,  0.35 

f(x)  -  0.09068948x  +  0.37645999x2  -  0.79837781x3  +  1.1306980x4 
Silt  - 

First  layer  moisture  range,  0.42 

f(x)  -  0.16554710X  +  0.32448206x2  -  3.5931164x3  +  6.0923182x4 
Second  layer  moisture  range,  0.30 

f (x)  -  0.16854730x  -  1.2164306x2  +  5.0226953x3  -  5.9828479x4 
Clay  - 

First  layer  moisture  range, 0.92 

f (x)  -  0.13851029x  -  0.32294633x2  +  0.51095728x3  -  0.23644523x4 
Second  layer  moisture  range,  0.46 

f(x)  -  0.14694158x  -  0.50856800x2  +  1.2942848x3  -  0.74715639x4 
where 

x  ■  soil  moisture  content  (in.  per  6-in.  layer)  above  minimum 

moisture  content  2  of  3  sheets 
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Table  4  (concluded) 


MAXIMUM  SOIL  MOISTURE  COMTEK!  AND  MINIMUM  SOI?,  KJiSTURE  CONTENT 
USDA  Soil  Textural  Classification  Tonus 
Field  maximum 
First  layer 

F.  max  -  2.01  -  0.013  S  +  0.132  OM  +  0.189  WI 
Second  layer 

F.  max  *  2.01  -  0.15  S  +  0.008  C  +  0.215  WI 
Field  minimum 
First  layer 

F.  min  -  -0.121  +  0.018  C  +  0.101  OM  +  0.105  WI 
Second  layer 

F.  min  -  0.170  +  0.025  C  +  0.013  0M  +  0.061  WI 
Unified  Soil  Classification  System  Terms 
Field  maximua 
First  layer 

F.  max  -  0.602  +  0.011  F  +  0.014  LL  +  0.181  WI 
Second  layer 

F.  max  -  0.221  +  0.014  F  +  0.011  LL  +  0.224  WI 
Field  minimum 
First  layer 

F.  min  -  -0.354  ••  0.001  F  +  0.024  LL  +  0.092  WI 
Second  layer 

F.  min  «  -0.094  -  0.001  F  +  0.025  LL  +  0.044  WI 


where 

F  -  fines  in  percent 
LL  «  liquid  liid.t 
WI  -  wetness  index  (see  table  5) 
S  *  sand  in  percent 
OM  ■  organic  matter  in  percent 
C  "  clay  in  percent 
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Table  5 

Classification  of  Sites  by  Wetness  Index 


Wetness 

Index 

Potential 

Wetness 

Depth  to 
Water  Table 

Depth  of 
Wetting 

General  Characteristics 
of  Sites* 

0 

Arid 

Indeterminable 

Less 

than 

1  ft 
(0.3  m) 

Located  in  desert  regions 

1 

Dry 

Indeterminable 

1-4  ft 
(0.3  to 
1.2  m) 

Steeply  sloping,  denuded,  or 
severely  eroded  and  gul¬ 
lied.  Mostly  semiarid  to 
arid  regions 

2 

Average 

More  than  4  ft 
(1.2  m) 

More 

than 

4  ft 

(1.2  m) 

Well-drained  soil  with  no 
restricted  layers  or  pans; 
fair  to  good  internal  and 
external  drainage .  Slope 
may  be  flat  to  steep 

3 

Wet 

1-4  ft 
(0.3-1. 2m) 

To  water 
table 

Soil  not  well  drained.  Re¬ 
stricted  layers  or  deep 
pans  may  be  present.  May 
occur  at  base  of  slopes, 
on  terraces,  upland  flats, 
or  bottomlands 

4 

Saturated 

Less  than  1  ft 
(0.3  o) 

To  water 
table 

Sites'  waterlogged  or  flooded 
at  least  part  of  year. 
Bottomlands  subject  to 
frequent  overflow.  Upland 
flats  with  poor  internal 
drainage  or  shallow  pans. 
Slopes  with  very  poor  in¬ 
ternal  drainage 

*  For  use  in  classification  when  water  table  and  wetting  depths  are  not 
measured. 
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APPENDIX  A:  DETAILED  FLOW  CHARTS  OF  COMPUTER  PROGRAM 


Appendix  A  contains  detailed  flow  charts  of  the  main  program  and 
subroutines  produced  by  the  computer  program,  FLOWGEN/F-II ,  which  is 
on  file  in  the  WES  Automatic  Data  Processing  Center.  Each  segment 
(subroutine)  is  detailed  separately.  The  flow  charts  of  the  main  program 
and  subroutines  are  as  follows:  main  program,  fig.  Al;  subroutine 
WRIBWK  (part  of  main  program),  fig.  A2;  subroutine  COMPMS,  fig.  A3; 
subroutine  FIL1DA,  fig.  A4;  subroutine  FIL2DA,  fig.  AS;  subroutine 
FRODEP,  fig.  A6;  subroutine  STFRSP,  fig.  A7;  and  subroutine  CALSST, 
fig.  A8. 


Fig.  Al.  Detailed  Flow  Chart  of  Main  Program 
(31*  pages) 
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QUIM2 


_ _ 35 _ 5 

r 

_ _ i 

7 

] FORMAT) 20X.22HS0IL  MOISTURE  CM/1S  CM)  | 

_ 37  _ i 

I  FORMAT I24X.14H0  •  6  IN  LAYER./ I  I 

- 1 

Ifl _  j 

I - 

£. _ 

I  FORMAT  I 24X  .  ISHS  -  12  IN  LAYER./ >  1 

1 

3.  j 

i 

1  F0RMATI2A3.IX.2A3)  1 
- ^ - 

. .  -i.  -  _  - - 

PRINT.  FOR  INFORMATION  CN  R'JNMINO  THIS  MODEL  SEE  USER  S  MANUAL  • ' 
PRINT.  . 'WHEN  ANSWERING  TES  OR  NO  QUESTIONS.' 

PRINT.1  USE  I  FOR  TES  ANO  2  FOR  NO.' 

PRINT.1 IS  INPUT  RAINFALL  IN  CM. 

READ .  tUNMcA 

PRINT.1 IF  YOU  NANT  OUTPUT  MOISTURE  CONTENT  IN  ' 

PRINT.1  INCHES  ENTER  A  1 <  CM  ENTER  A  2 

REAO.MEAOU: 

~~.=  ~  ~  s  _ r~ 

PRINT, -NAME  LAYERS  1  AND  2  SOIL  MOISTURE  CONTENT  PILES. 
PRINT.1  EXAMPLE  XXI001.XX200I 

REAO  39.  NMCLAI  .NMCLA2 

PRINT. -00  YOU  NEEO  OUTPUT  HOISTURC  CONTENT  IN  PERCENT. 
REA0.MEA0U2  _ 


CONT.  ON  PO 


56 


fll 


I  GO  TO  SO  > 


PRtNT. 'ENTER  THE  NflHE  OF  TOUR  INPUT  FILE'  | 

_  -  ' - 

I  READ  32.SBH  | 

CONT.  ON  PO  8 


'  q  7  n 


5 


RE BQ  STRENOTH  EOUBTIOK  COEFFICIENTS' 


CONT.  ON  PO  16 

PP  IS  OF  34 


65 


soft  1 )=RANO( 1 .11T.1  l/( CHflXI I  l-CXHINt  I  )  I 


I  AOFt 1 >= 1 /SOF (  I  )  | 


22  - 


215 _ 4 - , 

I  CONT  INUE  I 

31 


i=i 
jUG=  1 
JUB=  I 

'“I" 


C0T0160 


•READ  PRECIPITATION  OR  AIR  TEMPERATURE  OR  SNOH  CARD* 


STB 

|  REAOdOl  .7.£N0  =  305)  1QA.HO.NCOOE  1 
1  NQA= IQA«6  I 


lFtNTEMP.EO.NTOS.ANO.NSNOH.EQ.NTOS) 


V 

CONI .  OH  PG  27 


76 


78 


tO-lt.-V 


1F1HEA0U2-EQ.2) 


F 


COMPUTE  STRENGTH  equation  constants  frot 


81 


PC  3:  OF  34 


Fig.  A2.  Detailed  Flow  Chart  of  Subroutine  V7RIEWK 

(6  pages) 


r 
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j  FORMAT!  1X.12.1H-.12.4X.SI7H 

i.2F7.on 

_ { 

2 _ _ _ 

| FORMAT! IX.12.1H-.12.4X.6I7H 

I.F7.2I  1 

_ _ { 

_ _ _ _ _ 

|  FORMAT!  1X.12.1H-  .12.4X.CI7H 

1.F7.01  I 

I  =  JN I 


WRITE! 3.61 1IDHB . 1  DOE . ( STQRM01 I . JJ I . JJ; JNK . I  ORE I 


WR1 m  g  ,S3  » IQftB . IDRE .1 STORMCt I . JJ  )  .JJ-JNK .IDBE  ) 


8 


? 


t 

Fig.  13.  Detailed  Flow  Chart  of  Subroutine  COMM 

(11  pages) 
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SUBROUTINE  CQKPflS 


612 


61(j. 


JEN  =  2 
0QT0613 


SU — 2 - - 

|flCCR=0.  1 

HO  TQ  619  >c{jg 


? 


oifl _ & 


CONTINUE 

~1HZ 


ACCRETION  CALCULATION* 


X(  1  UCHAXI  1  1-PflCI  I  I 
X ( 2  );CHAX(2 1 - PMC  I  2 1 


CONI.  ON  PG 


I  SQTQ6iT~>- 

RaT1 


621, _ & 


111=1  I 

lJEH=4  | 


AS 


m — . - , 

XI  I  liRFI 


tfiCCRETIQNICLflSS  1  AND  CLASS  2l« 


A6 


JUJ _ 2. 


[  accwscac  u .  n»xm«c4i  n .  i » 


CONI.  ON  PC  5 


*JLJL-Jts 


9 


STORNDIHO.ITEHP 


CONTINUE 


Fig.  A4.  Detailed  Flow  Chart  of  Subroutine  FIL1DA 

(15  pages) 
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1  HBmfe.iin 


CON'l .  ON  PO  4 


107 


L.  M 


IF  i  o  i  ORMC  (  I.KKt.EO-O.Ol 


!  114 

t 


1 


CONT.  ON  PO  13 


pQ._n.ai... 


116 


lFinODlNYtSEO.18) -EQ.Ol 


X 

1 

f 

1  ISEO=I 

SEO*  t  1 

t 

rt  -- 

m 

I  NHCLfll 1 2  ):NUH8« 1SEQ/10«2S6«MQD1 iSEQ.ID) 


I  N(1CL«2(2)=NHCLflU21 
InEHFII  rO  I 


[ return  1 


I  END 
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Fig.  A5-  Detailed  Flow  Chart  of  Subroutine  FIL2DA 

(10  pages) 
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F ORHflT  < 14 > . 1 4hQ  ■  6  IS  1 


_ 38  6 _ 

|FOR!1flT I 24X , J5h6  -  T2  IN  LATFR./l  i 
- x - -1 


•test  output  checks  and  open  files« 


QmtP.EO.  IBr  .Oft.NEMF  II  .eojj_ 
T 


1 00  TQ  700  >- 

a - 


fCBU  OPENFl2.OPm.fllT] 
*- - 

£ 

fCfllL  SETSC_T1?<SECT3)  I 


■WRITE  HEflOlNOSi 


-ZCfl.-g 


MR] TE 1 2 .27 ) . ISINUH  | 
- — 

IwRiTEi^TzaTl 


HEflOU!  .EQ.»J^>-t>{|^j 


I  Hfl ! TE < 2 .37  l  I 


CONT.  ON  PC  3 
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123 


1088  = 


I0«E=7 

H*IT£(2.3I)  linoflllI.il.  11  =  1. 41.  ITS 
WRITE! 3.31  I  I  IH0«< 1 1  .1  1.11  =  1  .«  I.JYR 


Bm: 


HRPg  1  2  .27 » .  IS  IHUW 


IFlHEflOUl  .£0*2 ) 


rwR~n£l3.38) 


L  A 


Fig.  A 6.  Detailed  Flow  Chart  of  Subroutine  FRODEP 

(5  pages) 


131 


SUBBOUT  1NE  FRQQEP 


.... .WRITE  RIB  TEMPERATURE  ANO  SNOW  FILES* 

.....FORMAT  STATEMENTS*................... 


ju  ft 


[fORM  ATI  18X .24H1 DENT 1FICAT ION  NUMBER'  .40371 


22 _ ft 


FORMAT  I 3A3) 


_ 31  _ 

FORMAT  1 28 X . 3B3  .RI .X .2H19 . 1? 1  I 


32 _ ft - 

I  FORMAT! 2 A3  1 


32 _ 4 


FORMAT I6X.25HDEN5ITY'  0-15  CM  LAYER.  ,F5.2.23H0/CNi  15-30 

CM  LAYER.  .F5.2.4H0/CHI 


3A _ ft 


FORMAT! ; IX. 49HFR0ZEN  DEPTH  OF  SOIL ( CRREL - 1960 1  DATA  INPUT 
FILE./)  _ 


35  ft _ 

| FORMAT! 1H  ./I  \ 


32 _ 4 


j FORMAT! 26X.19HMEASURE0  SNOW  DEPTH.//)  | 


CONT.  ON  PO  2 


PC  1  ..OF _ i. 


132 


FQRHRTl  11X.4HSN0M.3X  .THPRECIP-  ,3X  . 3Hfl  1R  ,6X  .  1 3HSCI L  MOIS^RE  .6* 
13MS0II  MOISTURE) 


IF  I  |  .EQ.NMRXUl 1 t  ! 

l 


JENu=NHRXO I  2  I 


wRUE  ORTR  ON  FILES 


CONT.  ON  PO  5 
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Fig.  A7.  Detailed  Flow  Chart  of  Subroutine  STFRSP 

(5  pages) 
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CAUL  SE T SC T ( 1 .SECT  10  1 


CONT.  ON  PC  3 


Fig.  A8-  Detailed  Flow  Chart  of  Subroutine  CALSST 

(21  pages) 
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SUBROUTINE  CfllSST 

. '4 


....■COMPUTE  CONE  INDEX  AND  RATING  CONE  INDEX* 
.... .FORMAT  STATEMENTS*. 


Ji _ £ 


!  rQRHRTl?0X.?3HI0ENTiFICAT10N  NUMBER1  .4A3) 


29  _1 

7 

| FORMAT 

1 3 A3  1  | 
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r 
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:.3A3.A1  .X.2M19.I21  | 

_ _ _ 32 _ 3 

l 
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f 
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' 

| FORMAT ( 23X • 16H IS 
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31  5 

>  . 

1  FORMAT  1 24X • 14H0 

-  6  IN  LAYER./ 1  | 

CONI.  ON  PC  2 


PG  1  OF  21 
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FORMAT  1 203 . 1 X .203 . 1 X .203 . 1 X .203  I 


if 1 1 yr .gt - lay  i 


GO  *0  7S0 


TO  DETERMINE  SOIL  STRENGTH  REQUEST.. 

1 

L " '  

PRINT.  00  TOU  MONT  CONE  INDEX  FOR  FIRST  SOIL  LAYER.  SECOND  lOyER. 


PRINT.  RCI  FOR  FIRST  COYER.  SECOND  COYER.  SEPARATE  ANSWERS 
PRINT.  UJTH  COMMAS.  EXAMPlEFOR  RCI  ONLY  FOR  BOTH  LAYERS  :2  .2 . 1  . 1 
REAO.INSTEQi  I  > .1:1 .41 

PRINT.  ENTER  NAMES  OF  SOIL  STRENGTH  FILES  TO  BE  USED. ENTER 
PRINT.  WJOATO  IF  0  FILE  IS  NOT  TO  BE  USEO.- 

PRINT.  EXAMPLE  FOR  RCI  FOR  0GTH  LAYERS  sNOQATO .NQORTR ,xx 1  GO  I . XK2CQ I 
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continue 


.■•••test  OUTPUT  CHECKS  AND  OPEN  FUES« 
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CONT .  ON  PG 
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,g 


CONI.  OH  *0  S 
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1  ! 
i  i 

i  * 


CALL  SETSCTO.SECTei 
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F)4 


-e*- 


JZS 2 _ ft 


I  XLlti]=STKI.It  1  .1  1  I 

. . £ . 

I  XllH2  =  STRU<2.H  | 

<  00  753  J=K.LK  > 
- 

UBrl  | 


IF (  J.EQ.lBM.flNO.lYR.EQ.IBY  I 

T 


i 


CONT.  ON  PO  6 
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CO  TO  7S4 


ET 


L 


f * 


|  SlQiJ.l  UEKPiSTCON)  1 . 1  l«5TC0NI  2 , 1  UflLuCl  STQRWC  l'j  .l  .  i  i 


CONT.  ON  PC  7 
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ST  01  J  .1  )  :  XL  1 H 1 


I 


IFIHOOi  I.2I.EQ.0I 
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CONT  .  OK  PG  9 


151 


WRITE (NO. 27  l  ISINUH 


( ogjgr.TB?  „>c{jjj 


CONT.  ON  PO  10 


152 


-ffl 


CONT.  ON  PG  J3 

B..5  11.2  *.. 


1SS 


CBIL  SWVSCH2.SECT7I 


CONI-  ON  PC  19 


161 


CALL  SAVSC'I 2 .SECT7  ) 


APPENDIX  B:  COMPUTER  LISTINGS 

I.  This  appendix  contains  the  computer  listings  with  line  numbers  of 
the  main  program  and  subroutines  as  follows:  main  program,  table  Bl;  subroutine 
COMPMS ,  table  B2;  subroutine  FIL1DA,  table  S3 ;  subroutine  FIL2DA,  table 
B4;  subroutine  FRUDEP,  table  B5;  subroutine  STFRSP,  table  B6;  and  sub¬ 
routine  CALSST,  table  B7. 
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Table  B1 

Lletlm  of  Hein  Protrea.  FTWEHS 


10O*0VR.C0MPhS 

llOSOVR.FlLlOA 

l2u*0VR.FlL2DA 

130$0VR.r»0DEP 

1AOSo»/R»CAlSST 

150*OVR»STFRSP 

160*RPc 

i65C.».«. common  variables********************************************* 

1?‘j  COMMON  IDCi(35>»PMC<2),ST0(12«31>.nSTE0{4),RF{7), IHD«20>«CXH1NC2>» 
18,i4T£mPER(12, 35) , SN0( 12# 35) < C0NEW<3<2<  6) ,XR ANSI 2,3 ) »STCON (2.4 ) , 
l9.;4sTRLM2,4),NMAXO<2>.lSlNUM<4),C3<2.2>.C4<2,2>,DWMC<2.7>,COEF<i08>, 
20o4haxD(12)«  ImO* (4,12)»BD(2)»RANG(2»3.3).X(2)»C1(16«2),C5<18«4,2). 
5lP4CM*X(2),CHIN(2),C2(4),STORMC(l2.3l).STORMO(12,3l),AOF(2)»SDr(2), 
c2si4P  <2),KJnD<3),  1SLEW<3).DUM<25).PRECIP(12.35).B<13.2*?>»B1U2), 
23o4COn1(2»2)»D(2)»E(2>. ImOD ( 12 , 2 ) , ABAR ( 2,2 ) <BBAR(2 , 2) < NmCLAI ( 2 ) , 
2Ar)«NMCLA2(2)»NPMLAl(2),NPMLA2(2).NrRDEPt2).-NClLAl<2).NCiLA2(2), 
2!>u«NRCLAl(2)*NRCLA2(2).NS0UR. IUN.mEA . NTE«P< NSNOW.NYR.  I TEXA . 

26  c  4 1 'EXB,RMIN,1TyPE,ITEM1,ITFM2, II.IDEVl.ISl.ISD.ISK.lSY.llT.lDA, 

2 7 04m0 , IYR. 1TEhP,ICOUNT.ACCR»OEP,C1COF,C1ROF,NCODE,ACCR1.ACCP2»HOK. 
2Bo4jUl» JiL. JUG«LUB.J0K.M08# JEM, JAZ, JEL. J£B» 
36G4nTyPE,NEhLJM.NEWC0F,NS0ST.NR,LL» JUB. IBD, I0H. IBY.NEWYR, 
3/04MfcAOUl.MEAOU2,K,LK, JB, JEND.NTOS, A JBN. JBN. JN.KPMC.SECT1 « SEcTZ, 
38o4StCT3.SECT4. SECTS. SECT6.SECT7, SECTS. SECT9.SECT10. NCI. NRCl. 

3904SAM. 1ND. IDaB, IDAE, JNK. JNI .NO.STDN.NUMB. ISE0.NY1SE0.NY2SE0 

. . .  VARIABLES . •••••••*• . . . . 

45u  EOuIVALENCE  { NMCLAl . DMCLA1 ) . ( NMCLA2. OHCl A2 > . (NPMLA1 . nPML*i > , 

46,)  4  < NPMlA2» DPMLA2) . (NFROEP. DPRDEP ).(NCILAl.DClL*l).<NCIl  A2.DCILA2). 
47o4(NrClA1 .ORCLAl) . (NRCLA2 . DRCL A2  > 

475C*«.*.Data  INITIALIZATION*** . *..*•••••••• . ••••••• 

48q  DATA  NUMB. I SEO.NY1SEO.NY2SE0.NEMFIL.ACCR1.ACCR2/-001-. 0.0,0, 

49040. O.dl. 0.03/ 

500  DATAIN0A/12HJANUARY  . 12HFE8RUARY  .12HMARCH  . 

510412HAPRIL  .12HMAY  « 12HJUNE 

52c4i2hJuLY  , 12HAUGUST  , 12hSEPT£MBER 

53o4i?hOCToBER  , 12HN0VEMBER  ,12hDECEm9ER  / 

54,i  D*TAMAxD/3l,2fl.31.30.31.30.31.31.30.3l.30.3l/ 

550  DATA  COEF/Cr. 15970948.-0. 28827749.0.37098284,-. 13067845. 0..0.. 

5 7o4u- 20865303. -0.64984646. 0-67560732.-. 30*05206. 0..0.. 

5904. 15190803. -.07720229.-. 01280603. .P1851468.0. ,0.. 

6104. 0961229c. .01073670, -.10720756. .05548077,0. .0.. 

63040. 1164190. 0.06704479. -0.09112914,0. 03120642, 0..0., 

65o4. 143u4453.-. 10102288. -. 02103815, .10071059, 0..0.. 

6 7040.08401220, -0.01562037, -0. 10949073. 0.12228349,0. .0. . 

69o4o.  12713529, -0.53982006, 1.1252485,-0. 58561447. 0..0., 

7104. 11387171, -.65583730. 1.8700644. -2. 3253802. 

72041- 2952910. -.26246686. 

73040. 02789788.1. 24453131. -0.8606ni42, 1,2998033. 

7404-. 97359291'.  .29132112. 

7504.08321811. .11716883. -.14 710816.. 04910032. 0..0., 

7706.03492386. . 36829518. -.90 742210 >.61582018. 0.,0.. 

79oij. 12882001. -0* 3275105,1. 0509886, -j . 0462256 . 0 . • 0. . 

81040  *  O9 068948. 0.37645999,-0. 79837781 . 1.13C 696*0. <0.. 

03o4.165547iy. .32448206, -3. 593116*. 6. 0923182. 0..0.. 

85o4. 1685473. -1.2164306, 5. 0226953, -5. 9826479. 0..0., 
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Table  B1  (continued) 
Uetlns  of  Main  Program.  FTWEMS 


87u4. 13851029, -.32294633. .51095728.-. 23644523, 0 .  .0  . . 

89(j«.  14694:58. -.508568. 1.2942848. -.74715635.0.  .0./ 

91  i:  D*T*  BhIn. SfRL I ,  ClCOF. NEwCOF /0 .1. 750. #C . • 750 . .0. ,300. .0 . <300 . . 0 . . 
9158u , 0 / , 0/ , 

V204S ToRmC, STORED, TEMPER, SN0.ST0/372*0.. 372*0., 420a0,,42oe0..372ao./. 
93u4NrR.NMAX0/1.0.0/.C3/,47, .75, .22. .60/.C4/-.01,-.05,--,nl.-.02/ 

97n  0*T*  RANG/i. 26. 1.00,1. 67. 1.63, 1.47. 1.08. 1.05. 0.72, 1,55. 

98041.25.1  43,0.93.0.54, j. 35. 0.42, 0.30. 0.92. 0.46/ 

990  DATA  RF.PRECIP/.-'O  . , 420«C . / . MOK, JUL,  JIL.  JUG, LOB.  JOK.MOB,  JEH,  JAZ. 
10C04JPL. JEB/O.O.0.0  •  O.C.O.O.O/, JU8.KlND/0.3*0/, ISLEw/6193023, 

101038355711. 8355711.  ,  )C0F, NEWL I H.NSOST/O , 0 . 0/.NTOS/1H9/ 

1050  UaTA  C5 /o . 602. 0 . , u. • &. .0.011,0. 034. 0.»C. ,0.,0.,0.l3i ,0..0.,0..0.,0.. 
106080,221,0. ,0. .0. ,0.014,0.011,0. .0. , 0 . . 0 . . 0 . 224, 0 . , 0 . ,0 . . 0 . , 0 , , 

10704“0*  354, o.,  i*., o.,— 0*001. 0.024, 0«.r.*0.,0.,0, 0  92 ,0.,0.,0.*0*,0»* 
lOfl08'Q. u94, 3, , u. ,0. >-0.001,0. 025, 0.,t>.<0.,0., 0, 044, C,,0.«0«.0.<0.» 
109062. b I, *0.  (■  l3. 0. >0. >0..G*>0..0.«0.1 32, 0>,0. 189, 0.,0.,0>. 0*. 0,, 
lit.  08/. li, -0.0:5, 0., 0.008, 0..0..0..0. ,0.,0*. 0.215, 0.f0.,0.,0.,0., 

Ill U« - 0. l2l. L. . C. , u. 018. 0. , 0. . 0. . 0.. 0.101,0. >0.1 05. 0..0..0..0..0., 
112040.1 70.0. ,0 ..0. 025.0.. 0. , 0 . . 0 . . 0 . 0x3, 0 . . 0 . 061, 0. . 0. , 0 . . 0 . . 0,/ 
112bC“»**«FORMAT  STATEMENTS* ••**••*•••*••*•*••*••••••*•••••••• •****•••••* 

1130  1  FORmAT(2X,4a3.2X.2(A1»2X)»2F10.0) 

1140  2  FoRMAT(1X.F9, <  7F1 0 .  ) 

1150  3  FoRnATdX.  I2.2U.4X.6F7.3) 

1160  4  F o««AT(8f 7.3) 

11)0  *5  F CRHA  T  t  x  *  2C *3  ) 

11S0  6  F OPhAT ( x* 1 1 • 3 ( X. 12) ) 

1191  /  FoRMAT«iX.212.2X.Ai) 

i2. u  o  format ( a i > 

l2l 0  9  FoRMAT(8X.7F6.2) 

1220  10  FORMAT(ll) 

123')  ll  FORMA)  («X. 7F6. 2) 
l24»  1?  rOBMAT(lH  ) 

1250  ) 3  FORMAT(7x.l2.11.6Fl0.fl) 

1260  14  F0RmaT(4x,2(F5.3.x).4(2F6.0J) 

1270  15  f ORMAT (10x«3(2F5.2) J 
I2b0  i 6  f0RHAT(5x.4F7.3.5x.4F7,3) 

1290  17  FORMAT (3x,4(2F8.3) ) 

13'’”  18  FORMAT  (3x,12.|1.6F10.3) 

1310  19  FORMAT ( 10x. 14) 

1320  20  FORMAT (14» 

1330  /I  FORMAT(X. 12) 

134.)  2?  FORMAT  (5h1SITE.2X,4A3) 

135)  23  FORMAT!  6H  TEXTURE ,8X. Al, 4X.A1 ) 

1360  x4  F0RmAT(8h  DENSiTY.7x.F6. 2. 4X.F6. 2) 

137‘-  25  FGRMaT(7H  BEG  MC.4X.2F10. 2) 

:36(-  26  F0RmaT(16h  SOIL  PROPERT I ES » 5X, 5F10 . 3/2X , 7F10 . 3/2X.6F10 . 3 ) 

1390  .'7  FORMAT  ( llH  H*X-MlN  NC .  6x « 4F10 . 3  ) 

1400  28  F  ORHaT ( 3x . 2  5  2 ) 

1 4 : 29  F0RmaT(13H  COEFFICIENTS, 8x,4Fiu.3/2X.6Fl0.3/2X,6F10.3i 
14?<J  3p  FORMAK1X.F9.0.F1C.  ) 

14)9  3:  F0RMAT(1X.3I2»X.7F6.2) 

1440  )?  FORMA) ( A6 ) 

14^0  33  FORMA) (24X.15HC  -  l5  CM  LAYER,/) 

24- u  34  FORMAT (23x«16Hl5  -  30  Cm  LAYER./) 
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Table  B1  (continued) 

Listing  of  Main  Program.  FTWBMS 

1470  35  FORMAT (20X • 21HS0IL  MOISTURE  lN/6  IN) 

1480  36  FORMAT ( 20X « 22W3Q I L  MOISTURE  CH/15  CH> 

1450  3 7  F ORMAT <  24x » 14H0  -  6  IN  LAYER,/) 
lboO  3 8  FORMAT (24X,15H6  -  12  IN  LAYER./) 

1 t 1 0  39  FORMAT (2A3.1X.2A3) 

lbi5C‘-«***ADDlTlONAL  CONTROL  INPUT*********************************** 

1520  PRINT, "FOR  INFORMATION  ON  RUNNING  THIS  MODEL  SEE  USFR'S  MANUAL." 

1539  print. ♦."when  answering  yes  or  no  questions." 

1&40  print,"  use  1  for  yes  and  2  for  no," 

1550  Print, "IS  INPUT  RAINFALL  IN  CM?" 

1560  READ.IUNMEA 

1570  PRINT, "jF  YOU  WANT  OUTPUT  hGjSTURE  CONTENT  INI" 

i5po  print,"  inches  enter  a  1;  cm  enter  a  2" 

1590  READ.MEA0U1 

16P0  print, "NAME  layers  1  AND  2  SOIL  MOISTURE  CONTENT  FIlES." 

16)0  print,"  example;xxiooi.xx20ci" 

1620  read  39,  NMCLA1.NMCLA2 

1630  PRINT, "DO  YOU  NEED  OUTPUT  MOISTURE  CONTENT  IN  PERCENT?" 

1640  READ.MEA0U2 

1650  I f ( MfcAOUt . EQ . 2  )  00  TO  50 

1660  Print, "NAME  LAYERS  1  AND  2  PERCENT  MOISTURE  FILES." 

1670  PRINT,"  EXAmPLE:XX1001#XX2001" 

168"  READ  39,  NPMLA1 » NPMLA2 
1650  ?o  I OC1 ( 1 ) = 1 01 

l/oo  PRINT, "IF  YOU  HAVE  AIR  TEMPERATURE,  ENTER  A  T;  IF  Nt)T  A  9" 

1'/10  READ  8.NTEMP 
1/20  JOB-1 

1/30  PrInT."iF  you  HAVE  SNOW  DEPTh(M), ENTER  an  Sn;  IF  NOT  A  9" 

1740  READ  8 , nSNOW 

i 75  0  I f ( nSNOW . eQ. 3hS  )K I ND( 1 ) *10 

1760  print. "WHAT  IS  YOUR  INPUT  SOURCE?" 

1770  PRInT. "ENTER!  1  FoR  SURFACE  COMPOSITION  GROUP.  OTHERWISE  A  2" 

1780  HEAD.NSOUR 

1790  I F  f  NSOUR • EO . 1 )  GO  TO  60 

J.aoo  PRINT, "ENTER  THE  name  OF  YOUR  INPUT  FILE'- 
1810  «e*d  32. Sam 

1820  print, »aRe  max  and  min  soil  moisture  equations  included  in  the" 

1830  Print, "INPUT  DATA  F1lE7» 

1840  READ.MMDD 

let-o  call  openfu.sak) 
i860  go  TO  160 

1870  60  call  OPENF ( 1 , "DSURGR" ) 

i860  Print, "ENTER  THE  NUMBER  OF  YOUR  SURFACE  GROUP" 

1890  HEAD  20.NSCGNU 
iV.’U  Do  7C  1=1. 50C 
1910  HEAD ( 1 , 1 0 ) 1  TYPE 

1920  iFdTYFfc.EQ.J  )RE AD ( 101 » 19 )NCKSGN 
1930  I F ( ( NSCGNU-NCKSGN) > EQ>  0 )  «0  TO  115 
1940  /C  CONTINUE 

. . READ  SITE  AND  SOIL  IDENTIFICATION**"**********#"* . 

1950  U5  MOKsl 
i960  ITEmIsIDCKI) 

1970  READ < I TEM1, 1 )  ISINUH, ITEXA, ITEXB,BD<1>,BD<2> 
i960  PRINT22, ISINUM 
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199U  PRINT  23. ITEXA. ITfiXB 
2000  lrdTEXA.E0.3HN  >K1ND(2>»10 
20i0  lrdTEX8.E0.3HN  >KIND<3>«10 
20?0  PrInT24.8D(1).8D<2) 

2030  !rdTEXA-3HC  >122.120.122 

2040  120  lT(l>*KIND(2)sl 
2050  GO  TO  128 

2060  1 22  lr ( I TEXA-3HH  >125.124,125 
2070  J24  1 T{ 1  )sKIND(2 ) z2 
2080  GO  TO  128 

2090  12b  IP;  ITEXA-3HP  >128,127,128 
2100  127  MND<2)i2 
2U0  I  Td  >  =3 

2170  128  IP ( I  TEXB-3HC  >132.130,132 
2130  130  1  T  ( 2 ) =K !ND<  3 ) zl 
2140  GO  TO  150 

2150  132  IP ( I TEX8-3HH  >135,134.135 
2160  134  KIND<3)=2 
2170  I T (2  )  =  2 
2180  GO  TO  150 

2190  135  lr ( ITEXB-3HP  >150,137,150 

22y0  137  KIND<3>*2 

22l 0  d(2>=3 

22 20  150  JUL -1 

2230  JoKsl 

2235C*«»»«REaD  DATA  FROM  INPUT  FlLE***«*«***»****e*»*«****»««***»*** 

2240  160  REAOd.lOHTYPE 

2250  GOTO (115, 170, 180, 190, 190, 200, 210, 230, 240). I  TYPE 

2260  170  READ(l0l.21>NTYPE 

2270  Go  T0(172,174.176.176,179,175)NTYPE 

22 75C 3 ••••REaD  LIMITS»*»**»»**«»*»»*». •*••«••••*•••***••••*•••••«••• 

220'J  172  READ(101,14)  RHlN.ClCOF.  ( (STRLI  d  .J) ,  I»1.2>,  J*l,'4> 

2290  NEHLIM=NTYPE 

23u0  IP ( IUNMEA.EQ.l . AND.NS0UR.EQ.2)  RMIN«RHIN*0.3937 
23i0  GO  TO  160 

23i5c^.«*«READ  DEPLETION  EQUATION  . . ••••••• 

2320  174  lr (NSoUR.EQ.2)  GO  TO  1175 

2330  READ(10l,l3)  NS.NL, < CONEH ( NS . NL > JJ > » J J- 1 , 6 ) 

2340  NEHCOrsNTYPE 
2350  Go  TO  160 

2360  1175  REaD(101.16)  NS.NL>  <CONEH( NS.NL* 3.'?, JJ»1.6) 

2370  NE«COr=NTYPE 
2380  GO  TO  160 
2390  175  CONTINUE 

2400  READ(101,16X(C3(N»! >.N*1.2>.l«l,2>.((C4(N,l>,N*l»2>»Isl»2) 

2<20  Go  TO  160 

2425C««»**REaD  MOISTURE  RANGE  BY  LAYER  AND  SEASON*»»»«»*«»»»»«»*»**» 
2430  176  READ< 101,15) ( (XRANQC 1,J),1>1.2>,J*1,3) 

2440  Go  TO  160 

2445C* ••••RE AD  STRENGTH  EQUATION  COEFr ICJENTS*»»»«*** . ••*• 

2450  178  READdOl.l?)  <  (STCONt  I  *  JJ .  I*-i,2> ,  J=i» 4 > 

2460  NSOSTsNTYPE 

. . READ  IDENtIPICAtION  CODE  FOr  EQUATION  TO  COMPUTE  STRENGTH* 

2466C*»»»*  EQUATION  constant$»**»*»»»**»*»«**»****«******»*«***«»*« 
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2470  i7«  REA0(Ull»2fl)KCI.NRCl 
2480  NSOST=NTYpE 
2«90  r,o  TO  160 

2495C*****READ  BEGINNING  SOU.  MOISTURE •••*•••••••* 

2^00  180  REAO(101,30)  PMCU)  ,PMC<2> 

2»10  Print  25»PHC(l),PHC(2) 

2520  IF( IUNHEA.E0. 2. OR. NSOUR.EO.l)  GO  TO  184 

25. TO  Do  162  I«l*2 

2540  182  PMC( 1 )*PHC( I )*0 .3937 

2559  184  JUt*2 

2560  GOT0160 

. . READ  AND  PRINT  max  AND  mIN  SOU.  moisture 

2570  190  READ(lH»2)(C2(I).  1=1.4) 

2580  IF  < I  TYPE ”5) 192. 191.192 
2590  191CONT INUE 
26'.)0  U0  195  1=1,2 

26.0  Ji I *2 

2620  CMAx(l)=6.»BD(l)*C2(I)/100. 

2610  CmIN< I )=6.«BD( I )*C2( JJ/100 . 

2640  195C0NTINUE 

2650  192  lFUUNHEA.EO. 2. OR. NSOUR.EO.l)  GO  TO  3  94 

2660  Do  193  1*1.4 

2670  19J  C2( I )=C2( I )*G .3937 

2680  194  PrinT  27. (C2( I >. lsl.4) 

2690  Uo  196  1=1.2 

2700  Js 1*2 

2710  CmAX<1)=C2(I) 

2720  CxHINU  >*C2(J> 

2730  196  CONTINUE 

2740  If<NS0UR.EQ.2>  GO  TO  197 

2750  CALL  CLOSEf  <1> 

2760  Print, "ENTER  NAME  OF  RAINFALL  FILE" 

2770  READ  32,SaM 

27eo  call  openfu.sami 
2790  197  Go  TO  < 198 , 160 ) , JUL 
2800  198C0NT1NUE 
28)0  PmC(1)*CMAX(1) 

2820  PmC (2 ) sCHAX (2 ) 

2830  JuL=2 
2840  G0TOI6C 

2845c1 ••••read  soil  properties********************, 

2850  20OCONTINUE 

2860  I 0Eyl=10l 

2870  GO  TO  (201. H50), MOM 

2880  201C0NTINUE 

2890  MoK=2 

29‘.0  Oo  206  h*l»2 

29 1 0  ItEh2=IdEv1 

2920  ReAD<ITEM2.21)  II 

2930  lFUl.GE,8)G0  TO  202 

2940  READU01.3)  1 1  .NSCS.  J.  (Cl( l .  J) ,  I  =1,  II ) 

2950  GO  TO  2o3 

2960  202  HE AD ( 1^1 . 3 )  1 1 . NSCS. J, (Cl ( 1 . j ) . 1  =  1 , ft ) 

2970  2o3  iDEVlsl 
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“““"g  or  Program,  FTWEMS 

2980  lr( J 1-8)2 05# 205 #204 

2990  c04  READ<1.4>(C1(I,J),I=9»II) 

30.  0  >05  pRlNT  26.jCKI.mj.U1. II, 

30 1 0  206C0NT1NUE 

3020  Ir<MMDD.E0.2>  GO  TO  1150 

30.50  G0T0I6O 

30-iU  1 150CONT  lNUE 

3050  uo  1190  M  =  l,4 

3060  I IcM3= IDEV1 

3070  Ir<H.E0.1.0R.M,EQ,2 )N=M 

3080  ! f ( M . EQ . 3 )N=1 

3090  ir<H.EQ.4)Ns2 

31"0  If {hHDD.Eq,2)G0  TO  1180 

ll«  "MS™BE  "*x  *" 

3130  Ir(  I  I-8)H7o, 1170.1160 

^5°  REaD(1#4  )  (C5  < !  #  M#  NSCS )  #  I  =9#  1 1 ) 

3150  U7u  PRINT  29. <C5< I.M.NSCS). 1*1. 1 1 ) 

3160  I181)  CONTINUE 
3170  C2(H)=0. 

3180  DO  1190  Isl.JI 

3p’°  pi9!lnC2‘S,i:C2<M,*C5<I,M*NSCS)»Cl(I,N) 

3200  G°  TO  ( 194 , 160 ) JUL 

32o5C®.*..REaD  AND  WRITE  HEADINGS... ••••*..... 

3220  i-1  READ<l01»5,(lHD<J).J=l,20)  *  *  ******* 

3230  Go  TO  217 
3240  Do  220  1=1.3 

RE*D(1.5)(IHD(J).Js1,20, 

326?  21/CONTINUE 

3<«0  wRlTfc(2,5)  ( IHD( J) , J*l,20 ) 

33-0  <IHD(J).J=i,20) 

33..,)  22j  continue 

33 1  ■’  If  (MhA0Ul.EQ.2j  00  TO  231 
3320  WRITE(2.35) 

3330  WRITE(3.35) 

3340  wp l TE ( 2 » 37 ) 

3350  WRITE (3. 38) 

33o0  wr!TE(2»12) 

3370  wrItE(3.12) 

33? 0  Go  TO  232 
3390  231  WRI TE ( 2, 36  ) 

3400  wrITE<3.36) 

34i0  WRITE (2, 33) 

34?0  WRITE(3,34  ) 

34.30  wrJtE (2. 12)  , 

3<40  WrITE(3.12) 

3450  232  GO  TO16O 

3455C®. *..rEaO  SEASON  change  CARO.. * 

?On  '01'6>  ISl.ISD.ISM.JSY  . 

3470  IF  CJUB.EQ.DGO  TO  235 
3480  IBDsJDACKs JSD 
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Table  81  (continued) 
Listing  of  Main  Program,  FTWSMS 


34gi.  A  I BD*  !  8D 

35.  ;u  A  j8n*A18D/ '  . 

3b-u  irh*mck=Ism 
3520  NEWTS*  JfaYs  I SV 
3b^o  235  Continue 
3b -10  Do  215  1*1.2 

3b^i-  1 1  T s I T  ( 1  ) 

3560  IrtN6KC0F.LT. 2- OR. NEWCOF .GT.2)  GO  TO  225 

35  7')  SqF  =1 . 

35?'!  adFsI. 

359  :  Go  TO  2l5 

36. J  225  SDF ( I ) =R ANG t I » J IT  » 1 )/(CMAX ( I ) -CXH IN( I ) ) 

3<> :  u  *0F  <  I ) s l / sCF  (  I  ) 

36^-j  2i5  Continue 

3630  Irl 
3640  JuGsl 
3651)  Ju8  =  l 
36o0  G0TOI6O 

36c5C***«*REaD  precipitation  or  AIR  temperature  or  snow  carom 
3670  240  RE AD ( 11 1 »  7 . EnDs3  5)  IDA.MO.NCODE 
36flO  NdAs1da*6 

36*.n  If(nTemP.EO.NTOS.ANO.NSNOH.EO.NTOS)  go  to  2240 

3 7 .. . 0  lr!NTEMP.£Q.NTOS>GO  TO  2249 
3?; U  iFtNCoDE-NTEMP)  2249.2245.2249 
3/;-.  2249  I F  ( NCODE  -NSNOH  J  2240.2244,2240 
3730  2244  READU01.9).  (SNOtMO.  JJ).  JJ=!DA'NDA> 

374  Go  TO  3  ,0 

3 750  2245  READtl01.il). t TEMPER t MO . J J > . JJS I  DA* NDA > 

37.. C  Go  TO  3  0 

3//'  224,.  READtlOl.31)  IDA.MO. IYR. tRFt JJ>.JJ=1«7) 

3/60  1  r I  IBM ,E0 .MO . AND , !8Y .EQ . I YR. AND . 1 SD.GT. IDA ) IDACK= IDa 

3/9"  I  F  t  I  YR , l T . NEKYR )  GO  TO  498 

38j0  IF  l  IYR. CO. NEmYR.AnD.MO.LT. NCKJQO  TO  498 

38-  1  lFtIYR.E0.NEHYR.AND.H0.EQ.MCK.AND.5DA.LT.IDACK)G0  Tn  496 

3B;-o  kcXsMO 

383)  IdacK=1DA 

384';  call  eof Tst t i » mm > 

38S0  i F ( MM . EO . 2  >  GO  TO  24l 

38to  I F  t I YR , gT . NEHYR )  NMAXD 1 1 ) =NMAXD t 2 ) *0 
38-*0  I F  f  I  YR ,  GT  .  NEwYR  )  l8DrlBM  =  l 
3080  «U1  iFt 1UNMEA.E0.2)  GO  TO  247 
38=n  Do  246  Jjrl.7 
39n'3  RF(JJ)*rF(2J) *0.3937 

3V  5  246  lF(RFtJJ).GT.30.AND.RF<JJ>.LT.35.)RFfJJ)*77.77 
39  •  <)  247  IF  (MO.EO.2)  GO  TO  248 

39 v g  lFtMO.EQ.4.OR.MO.EQ.6.0R.M0.EO.9.0R.MO.E0.11)  GO  TO  249 

3930  I  COUNT r  31 

3940  Go  TO  250 

3950  243  I  COUNT  =  28 

3960  GO  TO  250 

39 7U  249  I  COUNT  =30 

39{- U  250  CONTINUE 

3990  N  TDA  =  I  0/ 

40n0  Do  25l  NJ=1.7 
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Table  B1  (continued) 
Listing  of  Main  Program,  FTWEMS 


40 i 0  PRtCiP<H0,NTPA)=RF(Nj)»2.5< 

40't'O  461  MDAsnT0A*1 

4025C  ■••••COHPoTfc  SOIL  MO  I STURE* •••••• »•••••••*•••*••*••*•••«•*•••< 

40.70  t*LL  L  I  NK  <  5 ,  "OCOMPM"  J 

4040  call  comPms 

40o0c». 

4060  3.J0  IF(JA7.E0.2>  GO  TO  301 
40  ;0  JgBsl 

«0fl()  301  GO  TO  (420.302). JA2 
4090  3(;2C0NTINuE 
41.  0  j:',  =  1 

4110  IF(MO.EO.NMAxDd)  .0R.M0.EQ.12)  GO  TO  303 
41,0  GO  TO  410 

41;-u  oo3  Go  TO  (304.304.3:-5)  JOB 
4 1 4 u  304  IF(nTeMP.EO.NTOS)  GO  TO  305 
4150  Jo8=JOB*1 
4lt0  Go  TO  166 

41a5C  3  ****WR I  TE  SOIL  moisture  FILES** *•••**»**••*••«•••••••••••••• 

4170  iOi>  CALL  LlNK<5.«0FIL10M> 

4ion  call  pilida 

419"  lF<MfcA0U2.E0.1)CALL  LlNK(5,"0FlI.2D"> 

4200  iF(MfcA0u2.E0.1)CALL  FIL2DA 

42iUC***** WRITE  FILES  FOR  INPUT  To  FREEZE-THAw  MODEL************** 
422U  IfInTEMP.GT.nTOS.OR.NTEMP.LT.NTOSJCALL  LlNK(5."0FR0nE") 

4230  If<NTEMP.GT.NTOS.OR<NTEmP.LT.NTOSICALL  frodep 
424"  JoBsl 

4246  I F (MEA0U2 . ED . 2 )  GO  TO  400 

424/C«****C0hPiITE  STRENGTH  EQUATION  CONSTANTS  FROM  SOIL  PROPERTIES 
4250  1f(NS0ST,eQ.5.AND.IYR.EQ.I8Y)  CAll  LINK<5."0STfRS"> 

4260  1F(NS0ST.EQ.5.AND.IYR.E0.!BY)  CALL  STFRSP 

. . COMPUTE  SOIL  STRENGTH*** . . . ••••••••••• 

4270  IF<NS0ST.E0.4.0R.NS0ST.EQ.5)CALL  lINK«5."0CALSS") 

4280  IF<NS0ST,E0.4.0R.NS0ST.EQ.5)CALL  CALSST 

4290  4 g o  JEL  =  1 

4300  JEB=2 

43)0  GOT0450 

43?U  *10  JEB=3 

4330  4?o  READ ( 1 » 10  I  ITYPE 

434"  Mo8=l 

4350  call  eofTsTu, iend) 

4360  IF< IEND-2>460.430.46C 
4370  4j0  M0B=2 
4300  GoTq470 

439"  460  IF(ITyPE-1)40O, 470.480 
44oO  470  GO  TO  ( 302* 495. 490 ) » JEB 
4*10  4e0  1 F <  ITYPE. EO. 3)  Go  TO  180 
4420  IF< ITYPE~8J240.230.240 
4430  490  JEL»2 
4440  450  GO  TO ( 420 » 495  ) . MOB 
4450  495  QO  TO  (115.5C0J.M0B 

4460  498  PRINT, "PRECIPITATION  FILE  HAS  DATA  OUT  OF  ORDER." 

4470  Print. "CHECK  D*TA IDA ■ MO. I YR 
448"  600  STOP 
4490  fcNo 
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Table  B1  (concluded) 
Listing  of  Main  Program.  FTUEMS 


4492  Subroutine  hribwk 

449 JC» •••••*••••••••*• 

449t>c». •••WRITE  FIRST  MEEK  OF  OUTPUT  DAT*  FILES  WITH  BEGINNING  DAY 
4496C“»«»*  OTHER  THAN  THE  FIRST  DAY  OF  WEEK****************** 

4550  61  FORMAT ( lx> I2> 1H-* 12, 4X«  7H 
4560  HI  FORMATllX. 12, 1H-, 12. 4X.7H 
4570  62  FORHATdX, I2dH-.I2»4X,2<7K 
4580  162  f 0RMAt<1x*I2.1H-,I2.4x.2(7H 
4590  63  FORMAT (IX . 12. 1H-, 12* 4X.3 ( 7H 
4600  163  FORHATdX.  12.1H-,  I2.4X.3<7H 
46i 0  64  FORHATdX.  12. 1N-.  12. 4X.4(7H 
4620  164  FORHATdX, 12, 1H-, 12. 4X,A(7H 
46.30  65  FORHATdX,  12. 1H-.  12, 4X.5I7H 
4640  165  FORMATdX»I2dH-,!2.4x,5(7H 
4650  66  FORHATdX,  I2.1H-.  I2.4X.6(7h 
4660  166  FORHATdX. 12. 1H-.  12. 4X.6(7H 
4670  UJNI 

4600  iFtSTDN.EO.l. )  GO  TO  8010 
4690  lr(!ND.EQ.6l)  hR1TE(2.61)IDAB.IDAE. <ST0RMC<! . JJ> » JJaJNK. IDaE) 
4/00  IF ( IND.E0.6l)  WRITE <3.61)10*6. IDAE. <STORHD( I. JJ), JJsJNK* iOAE) 
47io  lFdND.EQ.62)  WRITE <2» 62)10*8. IDAE. (STORMCI I , JJ> , JJaJNK. IOAE) 
4  720  lFdND.E0.62)  WR 1 TE < 3.62 )  IDAB.  I  DAE.  < STORMDd  .  JJ ),  JJsJNK .  IDAE) 
4  730  !FdND.EQ.63)  WRI  TE  <  2 . 63 )  I  DAB,  I  DAE .  (STORMCd  ,  JJ) ,  JJsJNK.  IDAE) 
4  74  0  IrdN0.EQ.63)  WR  I  TE  <  3 . 63 )  1  DAB » I  DAE ,  <  STORHDd  ,  JJ )  .  JJ«  JNK .  I  DAE) 
4 /SO  lFdND.EQ.64)  WR  I  TE  <  2. 64 )  I  DAB .  IDAE .  I  STORMCd  ,  JJ) ,  „Js  JNK.  IDAE ) 
4/6  0  lFdND.EQ.64)  WR I TE < 3.64 ) I D*B. ! DAE , < STORHD d » JJ ) , JJs JNK. IDAE ) 
4/7  0  lFdND.E0.65)  WRITE <2. 65)  IDAS.  IDAE.  (STORMCd  ,  JJ) ,  JJsJNK.  IDaE) 
4/80  lFdN0.EQ.65)  WR  I  TE  (  3 .65 )  1D*B .  I  DAE ,  (STQRMD  ( I  .  JJ )  .  JJsJNK.  I OAE ) 
4790  lFdND.EQ.66)  WR  I  TE  ( 2. 66 )  1 0*6.  IOAE ,  ( STORhC  ( 2 ,  JJ ) ,  JJa  JNK.  IDAE ) 
4800  JFdND.E0.66)  WRI TE (3.66) IDAB, I OAE • (STORHD( I . JJ) .JJsJNK. I CAE) 
48iU  Go  TO  8011 
4820  OolO  CONTINUE 

4830  lFdND.E0.6l)  WR I TE < NO. 161 ) ID*B. I D*E. <STO ( I . JJ) . JJ* JNK. ID*E> 
4840  lFdND.EQ.62)  WR  I  TE  (NO,  162 ) .  I  DAB ,  !D*E»  (STOd  .  JJ )  .JJsJNK.  IDAE ) 
4850  lFdND.EQ.63)  WR  I  TE  ( NC.  163)  IDAB ,  IDAE ,  ($Tod  ,JJ) ,  JJ*jNK ,  ID*E ) 
48ft  0  IFdND.EQ.64)  WR  I  TE  <  NO,  164 )  I  DAB.  1  DAE ,  (STO  (I .  JJ  > .  JJ»jNK,  ID*E ) 
4870  lFdND.EQ.65)  WR I TE (NO, 165 ) 1Ca8. IDAE, (STO ( I . JJ > . JJ«JNK, ID*E ) 
4880  1F(1ND.EQ.66)  WRI TE (NO . 166 ) IDAB, I  DAE. <STO< I .JJ) ,JJ*jNK. ID*E) 
4890  Boll  RETURN 
49i)0  EnD 


» 6F7, 2 ) 

, 6F7 , 0 ) 
J.5F7.2) 

) »5F7 . 0 ) 
),4F7,2) 

) .4F7.0) 
>,3F7.2l 
1.3F7.0) 
) « 2F7 ,2 ) 
)*2F7.0) 
>,F7.2) 

) »F7 , 0  ) 
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Table  B2 

Listing  of  Subroutine  CO HP MS 


60ou  subroutine  cohphs 

60 1  DC* •••*•• *••••••*«• 

6020cj«»»#CAlCuLAT10N  C0RE*******»«*» ••••••••••••••••«••••• 

6030  Do  6f>0  J*1.7 
60b0  I F  (RF  (J)«77.77)601,6'12,603 
6070  6ol  GO  TO  (604,605)  JOK 
608U  lF< ISk-mO>606,6U6,60? 

6090  fcu6  I F ( I S0“ 1  DA- J*1 )  608. 608.607 
6lu0  608  GO  TO  (640.605)  JOK 

61:0  °C4  JOK =2 

612'J  005  1  S3  ’Si 
6130  00  609  1*1,2 
614H  l  l  TsJ  T( I  ) 

6150  IF(NEwCOF.LT.2.Or.NEwC0F.GT,2)  GO  TO  610 
6160  CMlNt 1 )=CHAX( 1 )-XRANG( 1 , IS) 

6170  Go  TO  6  .*9 

618.'  6 1 0  CHIM  I  )=CHAX  <  1  )-RANG<  I » 1  IT,  !S>*ADF ( i  ) 

6190  609  CONTINUE 
620”  JuG*l 
62-0  o07  JCH*1 
6220  00  611  1*1.2 
6230  DgHC ( I , j ) =0 , 

6240  Go  TO  (612.613, 614, 615J.JEH 
6250  6i2  1F(RF( J)-RMIN)616. 617,618 
62oU  616  JEN=2 
6270  G0T0613 
6280  6 i 7  ACCR«0. 

6290  Go  TO  619 
63u0  6i8C0NT1NuE 

63: 0C=. •••ACCRETION  calCul at t on**»*< ••••••••••••••••• •••••• 

6300  X(l)*CNAX(l)-PMC(l) 

6350  X(2)*CMAX(2)-PMC(2) 

636  0  Ic(RF  (J)-X(D-X (2)  >621,620.620 

6370  620  H,2 

6380  JEM=3 

6390  G0T0614 

6400  621  11*1 

64 i 0  JEH=4 

6420  6l5  X {  l ) *RF ( j ) 

64 3 0c«.«»» ACCRETION! CLASS  1  AND  CLASS  2>***»* 

645”  6X4  ACCR*C3( 11. 1 >»XCI  )*C4{ 1 ! . I ) 

6460  IF(ACCR>622.622.623 

6470  622  lF(PMc(I)-CHAX(l))624,625,625 

6480  624  1F( 1 1 ,E0.l)G0  TO  626 

6490  PHC(l)*pHC(l)*ACCR2 

6500  GO  TO  627 

6510  626  PhC< I )=PhC( 1 )*ACCR1 
6520  G0T0627 

6530  623  fF(l!.EO.l.ANO.ACCR.lT,CiCOF.AND.I,EQ.l>  00  TO  628 

65«0  Go  TO  629 

6550  628  ClROF *RF ( J )-RHlN 

6560  IF (ClROF .LT.CICOF. AND. I .EO.l)  00  TO  630 
6570  Go  TO  629 

65aOC»»»«*ACCRE  riGN'fONE  TO  ONE  RELATION)***** 
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Table  B2  (continued) 
Listing  of  Subroutine  CO HP US 


660 630  ACCR  =  RF  (j)-RHIN 

6610  629  iF(ll.EO.l)  GO  TO  619 

662CC***** ACCRETION (CONST ANT )•••*• 

6640  IF( ACCR.GE. ACCR2)GO  TO  619 

6650  PhC( I >spMC< I )*ACCR2 

6660  Go  TO  627 

6670c  Ngfc  MOISTURE  CONTENT 

6690  61?  PHC( I )sPMC( I )*ACCR 

6 7()0  627  1F(PHCU)-CMAX(I))  631.631.625 

6/1O  62b  PMC( !  )=CMAX( I ) 

6720  631  IF < hEaDU2 . EC . 2 )  GO  TO  611 

6730C*****C0NVERT  MOISTURE  CONTENT  TO  X  DRY  HE  I GHT***e . . . 

6750  UhMC( I . J)=(PMC< I  >/{6.*8D< I >>  )*100. 

6^60  Go  TO  611 

6/80C«»»«»DEpLETION  CAlCuL AT ION*** •*••••••*••••••••••••**••• •••••••• 

680O  6i3  lF(PHC(l)-CM!N<ln  631.632.632 
6«10  632  IF (nEhCOF .NE.2)  GO  TO  633 
6820  DEPsO.O 
68 3  )  Uo  634  ls1 . 6 

6840  634  DEPrDEP*CONEH(!S.l.L)*(PMC(l)-CMlN(l))**L 

6850  IF(DEP)  631,631.635 

6860c 

6870C  I JKL  s  I  ♦MAxI*(j-1)*MAxI*MAXJ*(K-1>*MAXI*MAXJ*SAXK*(L-1) 
6880C 

6890  633  1 JKL*1*6* ( l”l)*12*{iTf| )  —  1 ) *36* < I S— 1 > 

6V0O  DEP=0.0 
69:0  Do  636  L*1 . 6 

. . TENTATIVE  AVERAGE  DEPLETION . 

6940  OE?=DEP«COEF( IJKL)*(SDF(I )*iPMC<I J-CMINU) ) >«*L 
6950  IjKL*IJKL*1 
6960  636CONT INuE 

6970C*****TRANSF0RME0  DEPLET ION***** 

6990  UEP*OEP/SDF(I) 

70oO  I F ( DEP )  631.631.635 

70i0C*****NEh  MOISTURE  CONTENT************************************** 

7030  635  PHC< I )=PMC( I l-DEP 

7040  IF(P«C< I)-CHIN(I ))637, 631.631 

7050  637PMC< I )=CHIN( I ) 

7060  G0T0631 
7070  6jl  CONTINUE 

7oao  6o2  ItehPiIda*j-i 

7090  IF(RF ( J) .LE.RHIN JGO  TO  638 
7100  STORMC(MO. ITEMPJiO.O 
7110  I F ( RF ( j ) , lT . 25 . 0 )  GO  TO  638 
7120  STORMO(MO, ITEMP)=0.0 
7130  Go  TO  639 

7140  63B  ST0RHC(M0.ITEKP)=PHC(11 
7150  St0RMD(H0. ITEMP)=PMC(2) 

7160  639  IF ( 1  TEMP- 1  COUNT )  600.640.640 

7170  640  JA2=2 

7180  Go  TO  645 

7190  60O  CONTINUE 

7200  JAZxl 

72:0  Go  TO  645 
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tfcble  B2  (concluded) 

Lis  cine  of  Subroutine  COMP US 


7210  OoJ  iTCMPstD*'* j-1 
72 ii)  Nh*XD!x  :*M‘' 

7240  Nh*XD(2)4j’E«P-l 
7250  JAZ»2 
72til  645  Rr.fURN 

7270  end 


Table  B3 

Hating  of  Subroutine  FIL1DA 


8o<-o  Subroutine  fjuda 

80 1  0Ca***** . •  •••• 

SOlbC-.^e-WRtTE  SOIL  MOISTURE  PILES . •••••••*< 

80} 6c°. •••FORMAT  STATEMENTS* •*••••: •••••*••••••••#• 4 

3020  i 7  FORMAT (20X.23H IDENTIFICATION  NUMBER:  ,4a3) 
8030  29  FORMAT ( 3 A3 ) 

8040  30  FORMAT  !  4  ( J.H  ./)> 

8050  3i  FORmaT(2x,3hDAY.24X,3A3»Ai,X.2h19, J2> 

BOfcO  32  F0RMATllX.I2,lH-.i2,4X,7F7.2> 

8070  33  FORMAT (24x»15hc  -  15  CM  LAYER./) 

8080  34  FORMAT (23x.16h15  -30  CM  LAYER./) 

8090  35  FORMAT (20X.22HS01L  MOISTURE  IN/  6  IN> 

8100  36  FORMAT (20X.22HS01L  MOISTURE  CM/15  CM) 

SUO  37  FORMAT  (24X.14H0  -  6  IN  LAyER.Z) 

8120  38  FORMAT (24x>15M6  -  12  IN  LAYER./) 
812t>C»***»TEsT  OUTPUT  CHECKS  AND  OPEN  FILES********' 
8130  IF(NEHFIL.EO.I)  GO  TO  650 
8140  IF(  IYR.gT.NEHYR)  GO  TO  651 
8150  Go  TO  652 

8160  65i  Call  openf<2,dmclad 
8170  call  setscT(2.sectd 

81 75C****» WRITE  headings***************************' 

8180  650  MRITE(2.27) ISINUM 
8190  t F (MEAOUl . EQ . 1  )WR I TE(2.35 ) 

8200  lFtMEAOUl.EQ.il  WRITE(2.37) 

82i 0  1 f ( MEAOUl . EQ , 2  )  WR!TE(2.36) 

8220  IF  (MEAOUl. EQ. 2)  WR1TE<2.33> 

8230  mr1TE(2.3'j) 

8240  I F l NEwF I L . EQ . 1 )  GO  TO  653 
8250  Call  0PENF(3,DMCla2i 
8260  Call  SETScT(3.SECT2) 

8270  653  HRITE(3.27)ISINUM 
8260  iF(MEAOUl.EQ.l)  wR!TE<3.35) 

8290  IF(mEaouI.EQ.I)  WRITE (3.38) 

83oO  IF(MEaou1.E0.2>  WRITE (3.36) 

83i 0  If(MEaOU1.EQ.2)  WRITE(3.34) 

83?0  wr I TE  c  3 »  3 o  > 

8330  NEWYRslYR 
8340  652  K=lBM 
8350  LK-MO 
8360  DO  654  I'K.LK 
8370  JB®1 

8380  IF(IBD.GT.1.anD.I.EQ.I8M.AND.IYR.EQ.IBY)JSsIBD 
8390  I0AB*1 
8400  IDAE«7 

8410  I F ( I «EQ.NmAXD< 1 > )  JEND*NMAXD( 2 ) 

8420  IF ( 1 1 EQ. NM*XD( 1 ) )  GO  TO  655 
8430  JENO«MAXD(l) 

8440  655  IF (MEAOUI.EQ.1)  GO  TO  656 

8450  DO  657  JsjB.jEND 

8460  STORMCt l » J)*ST0RmC( I . J)*2.54 

8470  65?  ST0RMD ( I . J )*STORMD( I . J)*2 <54 

8460  656  WRITE(2.31)  (IMOA(H.I),  11*1.4),  I YR 

8490  WRlTE(3.3l)  ( I  MO A { 1 1 . I ) >  11*1. 4), IYR 
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85qU  JRN»A JBN 

85i 0  Uo  698  NH*1,5 

8520  lrUTR.GT.I8Y)  GO  TO  659 

8630  InKjO 

8540  IFUJBn.GE.NW.AnO.I.EO.IBHJGo  TO  660 
8550  IrCJB.EO •1iOR.JBN.LT. (NW-1) )  GO  TO  659 

8555C •••NR 1 Tg  OUTPUT  HHEN  BEGINNING  WEEK  IS  ONLY  PARTIALLY  TILLED 
8560  Uq  661  nK=1DA8.1DAE 
8570  iNKaINK*). 

8580  If<STORhC{ l#NK).EO«0.0)G0  TO  661 
8590  1nDc6C*|NK*1 
86o 0  JnKsNK 

86it>  JMl>l 

8620  Call  hriBnk 

8630  Go  TO  660 
8640  661  CONTINUE 
8650  Go  TO  660 

8660  659  lr(I>EQ.2.AND.NW.E0.4)NWsNW*l 
8670  IT(NN.EO«5)10AE=HAXD(I) 

8660  JU=jEND/7 

8690  lTU.EO.NHAXD(l>.ANO.JN.EQ.Nu)IDAE>UENO 

8696C»»*«# WRITE  DATA  ON  T I LES**»*«*'« **•••••••••••••••••*••••••••••• 

t7o«  HRlTE(2.32)  IDAB.IDAE.(STCRKC(I.JJ)»  JJ-IDAB.  1DAE) 

8 7 i 0  wrItE(3.32)  IDAS. IDAE. (STORNDt i. JJ).  JJ*IDaB. IDAE) 

8720  1 r c I . EQ.NhAxD( 1 ) . AND . JN.EQ.NW)GC  TO  662 
87j0  660  IDAB«IDAB*7 
8/40  658  I0AE»IDAE»7 

8750  662  IFU.E0.6.0R.1.EQ.12)  NR!TE<2.29)  JSlEH 

876o  IF(I<E0.6.CR.I.EQ.12)  MRITE(3.29)  ISlEm 

8770  lrU.LT. 6. OR.  J.GT.6)  GO  TO  663 

8780  hRlTE<2,27).ISINUM 

8790  1 F  <  Ht  AOul .E0.1)HR1TE(2<35) 

681,0  1f(HEA0U1.E0.1)UR1TE(2>37) 

8B1U  I F  < Mt AOUl . EO . 2 )  HR  I TE ( 2 » 36 ) 

8B2O  ;fimEA0u1-E0.2)  uRlTE(2>33) 

8630  wrITE(2.3u> 

8840  wrITE(3.27)  ISINuh 

8850  lF«HEA0Ul.E0.1)NRITE(3.35) 

6860  I ~ ( ME AOUl • EO . 1 )  HR  I TE ( 3 . 36 ) 

8670  1F(HEA0U1.E0.2)  HR  I TE ( 3, 36 ) 

8880  IF (Ht AOul . E0.2 )  HRITE<3.34) 

88911  wRlTfc(3,3C) 

8900  663  IF (MEA0U2.E0.2)  GO  TO  654 
89;0  1F<MEaOU1.EQ.2>  GO  TO  664 
8920  UO  665  J«JB.JEND 

8930  STO««C< ! . JJ=STORHCU • J)»100./(BD(1 )*6.  > 

8940  665  ST0RMDU.J)sST0R*D<I.J>»inC./(BD(Z)»6. > 

8?5u  Go  TO  654 

8960  664  continue 

8970  DO  666  j«jB. JEND 

8980  STORHCU  .  J)=STORMCU.  J)«100./(B0(1)«15.24) 

8’90  t&6  STORHUU. J)=STORHD(I»J)»100./(BO(2)»15.24) 

90oO  654  CONTINUE 

90:0  call  savSctiz.secu) 
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9020  1F(IYR,E0.1BY.0P.NEwFIL.EQ.1)GG  TO  667 
9030  CALL  CLOSEft?) 

9040  GO  TO  668 

9050  667  CALL  CLOSED (2.DHCLA1) 

9060  668  Call  S*VSCT(3.SECT2) 

9070  IF(  ITR.E0.  I8y.0R.NEwFIL.E0.DG0  TO  469 
9080  Call  closef  <  3  > 

9090  Go  TO  670 

91  O'*  069  CALL  CLOSEF  <  3 ,  DHCLA2 ) 

91 1 1*  670  NY1SEOSNY1SEO*1 

9ll5c»**»*NEED  FOR  NEW  CONTINUATION  OUTPUT  F 1LE ••••••• ••••••••••••• 

9120  1f<HOD<NY«SEQ.18> .EQ.O)  1SEQ=ISEQ*1 

9130  1F(N0D(NY1SE0.1S) .E0.0)NHCLAi(2>=NUHB*!SEQ/10*256*MoD(ISE0.10) 
9140  NmClA2(2)=NHCLA1(2) 

9160  NEWFlL  =0 

9160  IF t MOOt NYlSEQ. ie) . EQ.O )NEHF1l=1 
9170  RETURN 
9i«o  end 
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lOnOj  SUBROUTINE  FJL2i)a 
IOoIOC. «•••••••••••..•• 

19015C. ••••WRITE  PERCENT  MOISTURE  FILES*********..*..*.*...**.**.... 
.0U6C.. ...FORMAT  STATEMENTS**********. •••<•••••••••». •••••••. •••••• 

10020  27  FORM At ( 20 X.23H IDENTIFICATION  number:  .4A3) 

1C0A0  26  F ORmAt C  19x » 24HPERCENT  MOISTURE  CONTENT; 

10f4i3  29  F  ORMAT  t  3  A  3  > 
lo-jbt  30  format ( 4 c ih  ,/n 

1016c  31  F0RM*Tt2X.3MDAY.24X.3A3.Al.X.2H19.I2) 

10>j7i  32  F0RmAt(1x.I?.1H-.12,4X.7F7 .2) 

10n»u  33  f ORmAt ( 2*x. 15h0  -  15  CM  LAYER./) 

ICoVu  3*  F0RHAT(23X.16H15  -30  CM  LAYER./) 

1010.  3/  F0RMAT(24x. 14HO  -  6  IN  LAYER./) 
lOlle  38  F0RMATS24X.15H6  -  12  IN  LAYER,/) 

10U5C...*.TEST  OUTPUT  CHECKS  AND  OPEN  F IlES***.********.**. •*••«•«• 
10i2o  IFcIyr.Eq.IBv-OR.NEwFIl.EO.DGO  TO  700 
10l3o  CALL  OPEnF (2.DPMLA1) 

10-4u  CALL  SEtSCT ( 2 » SECT3 ) 

101 45C ..•••WRITE  HEADINGS* •*...•••...*.*.•••*..••. ••*••••••••••••••• 

lOibu  70u  WRl Tg 5  2,27 ) , 1SINUM 
10i6u  WRITEC2.28) 

101/0  lFcMEAOUl.EO.il  wRl T£  c2 , 37 ) 

IO180  IF  CMEAOUl .EQ . 2 )  UR1TEC2.33) 

IO1V0  WRITEC2.3C) 

10/Ot-  iFclYR. EO.IBY-OR.NEwFlL.EO.il  GO  TO  701 

10211  CALL  OPEnF C3.DPHLA21 

10>2‘j  CALL  StTsCT  c3«SECT4i 

10230  70i  WR I Tg C 3, 27 ) , 1 S INUM 

1024k  mRI TE C 3 . 28  1 

10,50  IF CMEAoUl.Eo.ll  wR!T£c3.3B) 

10260  If CMEAOUl. EQ. 21  WR1TEC3.34) 

102/0  WRHEC3.301 

1028.1  DO  702  1=K.LK 
10290  j8=l 

10300  iFf IbD.Gt.i.AND.I .EO.1BM.ANO.IYR.EO. 1BY)JB*IBD 
10310  IDaB=1 
10320  IDAE:7 

19330  WRITEC2.31)  CIMOACII.D,  II  =  1.41,IYR 
1034(1  wR  I  T  E  C  3*  31 1  ClH0ACl!.l>.n*1.4),IYR 
lOjby  DO  7 03  mvjsi.5 
1036c,  IFclYR. GT.IBY1G0  TO  704 
103/O  !NK*0 

10380  JFCAJBN.GE.NW. AND. I. EO. IBM)  GO  TO  705 
1039--  IF ( JB. EO. 1 .OR.JBN.LT . CNW-1) )G0  TO  704 

10395C... ..WRITE  OUTPUT  WHEN  BEGINNING  WEEK  IS  ONLY  PARTIALLY  FILLED 
1 0  -s  0  »j  no  706  Nk=  I  DAB.  I  DAE 

1041.1  INK»INK*1 

1042o  JF  c S 1 ORMC < I , NK 1 . EO . C . 0 )G0  To  706 

10<3(!  IND«60*lNK-l 

1 0  4  4 .  jNK  =  NK 

lC45o  jN I s I 

10460  CALL  WRiBWK 

104  Z i'  GO  To  7C.5 

10481,  706  CONTINUE 
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104V.'  GO  TO  70b 

lOjO,-  704  IF ( 1 .E0.2.AND.NW.EQ.4)  NH«NW*1 
1 0 1  <  [F  (Nw.E0.5)  IOAEsMAxD(I) 

10’>2ii  jNrJENO/7 

10541  IF  (I  .EQ.NHaxOU)  .  and . jn .eq . nw j idaes jend 

106450»4**«WR1t£  DATA  ON  FILES**** . . 

1054o  wRIT£(2,32)  IDAB,  IDAE,  (STORMCd.JJ).  JJs l DAB, I  DAE ) 
l05b>.  uR I TE (3,32 )  I  DAB ,  1  DAE ,  ( STORnDd  ,  JJ  > ,  JJdDAB,  IDAE) 
1 0 5 6 ij  IF(1.E0,NMAXD(1).AND.JN.E0.NH)G0  TO  707 

105  7n  705  IDAB=1DAB*7 
1058m  703  1DaE:IDAE*7 

lO'-’Vo  707  1F(  I  ,  EQ.6.0R.  I  .  EQ.12)  WRITE(2,29)  ISLEw 
lOftOi.  |F  ( I . E3.6.0R.  I .E0.1?)MRITE(3,29)  1SLEW 

1061. :  IF(I.LT.6.0R.I.GT.6)  GO  TO  702 
1062u  hRITE(2«27), IslNUM 

lOt4o  wRIT£<2»26) 

1 0 6 4 ij  jF  (hEAoUl  >Eq.  1  )UR!TE(2, 37) 

1065.. .  JF  (MEA0Ul.E0.2)HRlTE(2*33> 

10660  wR|TE(2,3i‘) 

1 0 6  / (j  wR  I  T£  (  3 1  2  7  )  I S  I  fJUM 
10680  wRITE(4,28) 

10690  IF (HEAOU1 . EQ . 1 > WR l TE ( 3 • 38 > 

10700  iF(HEA0U1.E0.2>HRITE(3,34) 

1071:  wRITE(3,3U 

1072>  702  CONTINUE 

1073, >  c*LL  SAVSCT ( 2 , SECT3 ) 

10  74..  iFdYR.EO.  JBY.OR.NEMFIL.EO.1)  GO  TO  708 
10/bi.  CALL  CLOSEF  (  2  ) 

1076..  G0  To  709 

107/,,  708  CALL  CL0SEFC2.DPMLA1) 

10/8.J  709  CALL  SAVSCT  ( 3 .  SECT4  > 

10’9.,  IFdYR.EO.  Tfl'Y. OR. NEHFIL. EG. l)GO  TO  710 
10/ 'O'.  C*LL  ClOsEF  ( 3  ) 

10810  GO  To  711 

1062.  71(1  CALL  CLOSEF  (3.DPHLA2) 

. . TEST  NEED  FOR  NEW  CONTINUATION  FILE . . 

I064u  711  NPMLAl(2)=NMCLAl<2) 

1 0 P 4 1,  nPMLA2(2)=NMCLA1(2) 

10c5l.  nEhFil  =  o 

1 0 B 6 iF<HOD<NY1SEO«18).EO.O)  NEWFIL=1 
108/0  RETURN 
lope,-  end 
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120U0  SUBROUTINE  FRODEP 
120U1C. •••••• •••*•«•••• 

12002C. ••••WRITE  AIR  TEMPERATURE  AND  SNOW  F ILES*»««*.***. ••**«*••*.*•••* 
12.030 . format  statements******************************************** 

12()1m  2>  FORMAT(10X.24H1DENTIF!CAT1ON  NUMBER:  ,4A3) 

1202  i  29  F  ORMAT  <  3a3  ) 

l2l)3j  31  F0RmAM28X.3A3,A1,X.2H19.  12) 

1204.,  32  F  ORMAT ( 2A3  ) 

l2obt  33  F0RmAT(6x.?5H0ENSITy:  0-l5  CM  LAYER.  ,  F5 . 2 . 23HG/C>'. )  15-30 

1 2 0 6 •'! 4  CM  LAYER.  .F5.2.4HG/CM) 

120/..  34  F0RMAT(iiX.49HFR0ZEN  DEPTH  OF  S01L<CHREL-1968>  DATA  INPUT 
I208u4  FJLE./) 

120V0  35  FORMATIIH  ./) 

1 2 1 0 1.  36  F0RMAT<26X.19HMEASURED  SNOW  DEPTH.//) 

1 2 1  lo  3/  rORMAT(llX.4HSNOW.3X,7HPREClP-.3X.3HAIR,6X,l3HSf)IL  M0ISTURE.6X, 
l2l2o*l3HSoiL  MOISTURE) 

I2l3u  38  F0RMAT(6x.  1 2 , 3 < 2x » F6 . 2 ) , 8X , F7 . 2. llX , F 7 . 2 > 

1 2 1 4 y  39  FORMAT(6x.3HDAY.2X.5HDEPTH.2x,7HITATION.3X.4HTEMP.4X,i5HO-l5 
I2l5n&  Cm  LAYER. X.3X.17H  15-30  CM  LAYER.*) 

I2l6y  iF(IyR.Ot.IBY)  GO  TO  850 

12165C.****A0D)TJ0NAL  INPUT  CONTROL************************************* 

l2l >\j  PRINT, "ENTER  names  of  snow, frozen  soil  depth  model  input  file" 

I2l8j  PRINT,"  EXAMPLE:  FRD001" 

1 2 1 9 u  READ  32 , NFRBEP 

12iV5C.****TEST  OUTPUT  CHECK  AND  OPEN  F I LES*** ••»*••*•*•**•**•*••**••*•• 

12280  850  1F( IyR.EQ. IBy.OR.NEWFIL.EQ.1 )  GO  TO  851 
l2c?l  )  CALL  0PEnF(4.DFRDEP) 

1222j  CALL  SETSCT ( 4 , SECT9 ) 

1223  >  851  CONTINUE 

12235C. ••*. WRITE  HEADINGS*********************************************** 

12240  DO  852  I=K,LK 
1225U  WRITE(4,27). ISINUM 
1226,1  wR!TE<4.34) 

122/ij  lF(KINDtl).EO.O)  GO  TO  853 

1228.,  uR  1 1  E ( 4 ,36 ) 

1229.,  853  WR 1  TE  ( 4 , 31 )  ( IMOA(  11,1), 11=1. 4), lYR 
1230-'  wRITE(4,33).BD<1),BD{2) 

1 2 *  1 1*  WRITE(4,35) 

1232m  wR)TE(4,37) 

1233o  wfi I TE ( 4 , 39 ) 

12340  jBi IBD 
l2’5n  jEND=MAXD(I ) 

12360  IF  ( l  .EQ.NHAXDU)  )  JENO=NMAXD( 2 ) 

12365C.***«WR1 Tfc  DAT*  ON  :'ILES*********** *•*•*••••••**•*•••**••*••••••*• 

123/fJ  DO  854  JsJB.JENO 

1 2 3 8 (.  854  HR ! T£ (4, 38) J.SNOC I , J) ,PRECIP( 1 , J) , TEMPER 1 1 . J) »STORMC< I . J) . 

l239r.4sToHMD(  1,  j) 

124 0 J  852HR I TE ( 4, 31 )  ISLEW 
12410  C*LL  SAVSCT ( * . SECT / ) 

I242y  IF(IyR.EQ.IBY.OR.NEwFIL.EQ.I)  GO  TO  855 
12430  CALL  ClOSEFM) 

1244U  GO  TO  856 

1245C  855  CALL  CLOSEF t 4 , OFRDEP ) 

12460  856  NY2SE0sNY2SEQ*l 
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124/u  If (M0D(NY2SE0»3) .  EO.O) 1 SEQe 1 SEQ+l 

I24d0  |f (MOD(NY2SEO#3).EO.O)HFROep(2)=NUM8*ISEQ/10»256*MOD(ISEQ»10> 

124D5C . TEST  NEED  FOR  NEW  CONTINUATION  F1LE««»»« . •  . 

12490  NEwflLsO 

12b0u  If (MOD(Ny2SEO.3).EQ.0)  NEwftl«l 
12510  RETURN 
1252ii  END 
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l5oOo  SUBROUTINE  STPRSP 
l5ulOC. •*••»••••••••*.*. 

. . COMPUTE  STERNGTH  EQUATION  CONSTANTS  PROM  SOIL  PROPERTIES 

■  l5o2'J  34  FORMAT (5x»21H INCORRECT  Cl  MOD  CODE) 
l5«3o  CONl  ( 1 » 1 )  3  A(.OG  ( 200  . ) 

1&0<0  CONlU.2)=A|_OG(100. ) 
l5i;5j  CONI  ( 2 . 1 ;  iAlOG  ( 300  .  ) 
l5o6u  CONI.  ( 2  >  2 )  sAtOG  ( 200  . ) 
l5n/o  CALL  SAVSCT ( 1  ,  SECT10 ) 
l5od(,  CALL  CLOSEP(l) 

15u85C.**»«sEArCH  POR  STRENGTH  EQUATION  CONSTANTS****************** 
l5o9o  CALL  0PENF(1,"DFST£Q") 

13i  on  Do  800  1 1  si  >  NC! 

lSll.t  READ  (1.)  (  I  MO  D  (  I « 1)  < 1*1, 12 ) 

15i2(i  Do  800  J*1.2 

li>lde  READ  (1.)  (B<I,J.l),I=l,l2> 

15140  800  CONTINUE 

I5l50  I2*7-II*NRCI 

l5l6i)  Do  801  11=1.12 

lSl'O  READ  (1.)  <1M0D<I,2>,!=1,12> 

15i8o  DO  801  Jsl > 2 

15190  READ  (1.)  (B( I > J,2) . 1*1(12 ) 

15200  801  CONTINUE 
15210  CALL  CLOSEP(l) 

1522(1  CALL  OPENF (1,  SAM) 

15230  CALL  SETSCT(l.SECTlO) 

lS:J2C.#*N  INDICATES  LAYER. LI  IS  CONE  INDEX, L2  IS  RATING  CONE******. 

15240  DO  802  N=1.2 

l5?5(i  dO803L  =  1.2 

1526(i  D<L)=B(1.1.L)-B(1,2.L) 

152/C  E<L)=B(1.1.L) 

1528u  008031=2.12 
15290  jTEMPslMOD(I-l.L) 

15300  TEMPsCKI.N) 

15310  1P(JTEMP)604.80o.805 

15320  805CONTINUE 

15330  IP (J TEMP-1 >804,806,807 

15340  807CONTINUE 

1535u  IF (JTEMP-2 >804,808, 809 

1536u  808C0NTINUE 

15J65C.»»*. COMPUTE  LN  OP  VARI ABLE********************************** 

153/0  TEMP=AL0G<TEMP) 

15380  G0T0806 

15390  809  IF (JTEMP-3 >804,810, 804 

15395C. ••••COMPUTE  LNLN  OP  VARI A6lE*»*****«*< *•*••••••*••••«••••*•• 

15400  810  TEMP=ALOG(ALOG(TEMP>> 

1541,i  GOT0806 
1542,,  804  PRINT  34 

1543j  806  D(L)=0{L)*(B<I.l.L)-B<l,2»L)>»TEMP 
1544o  E<L)=E(L)*8(I.1,L)*TEMP 
1545i.<  803C0NT I NUE 
15460  NN»N 

15465C.*»». COMPUTE  A  AND  B  OP  STRENGTH  E0UATtON********************< 
154  /ij  d08J2L=1,2 
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l54  8o  aBaR ( L  <  N )  "CONI  •  l.L)*(C0Nl{2,L)-C0Nl(l»U  )«E(D/0(L) 
154 Vu  bBAR (L  »N)  =  (CONI ( 1 . L )-C0Nl(2, L ) >/D(L ) 
l5bt)o  PRINT, "LAYER  r  "  ,  N 
1 55 '1 1  1 P  ( L  .  EO 1 2  )G0  TO  820 

155U2  PRINT,"  CONE  INDEX  EO.  A«  ",AB/R(L,N> 

155J4  PRINT,"  B*  ",BBAR(L»N> 

15505  GO  To  630 

15506  82u  PRINT,"  RATING  CONE  EO.  A«  ",ABAR<L,N> 

15500  PRINT,"  Ba  ",BBAR<L, N) 

15510  83d  STCON  (1*NN)=ABAR(L»N> 

l552u  STC0n(2*NN)=BBAR(L,N) 

1553u  N'NrNN*2 
15540  802  CONTINUE 
1555.,  RETURN 
15560  END 
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1-JuOO  SUBROUTINE  calcst 

l3ol‘>C  ••••••••  •*••••  •• 

13.U5C.***. COMPUTE  CONE  INDEX  AND  RATING  CONE  INDEX**********************. 
1 J..16C** •••FORMAT  STATEMENTS********* ••«•*• ••••••**•«•••*»•*•••••••••*•••< 

lJj/.i  2 7  F ORmAT(20X.23HIDENTIFICATION  NUMBER:  ,4A3> 
l3i;3ij  29  F  ORMAT  <  3  A 3  ) 

I3'i4ij  30  F0RmAT(4(1H  ./)) 

13.. 5,  31  F0RmAT<2x.3HDAY.?4x  3A3. Ai » X, 2H19, 1 2 > 

13o6.j  32  FORHATdX.  I2.1H-,  12.4X.7F7.0) 

1 3  o  /  o  33  FORMAK24X.15H0  -  15  CH  LAYER,/) 

1 3 1 •  0 u  34  F0Rh*T(23x«16H15  -  30  CH  LAYER./) 

13:>9.;  37  FQR.lATI24X.14H0  -  6  IN  LAYER,/) 

13 1 0.i  38  F0RMAT(24X.15H6  -  12  IN  LAYER,/) 

1 3 1 1 v  39  FoRmAT(26X.10HC0N£  INDEX) 

13)  40  F0RmAT(23X.17HRATING  CONE  INDEX) 

1 3 )  3  •  j  41  F0RmAT(2A3,1X»2A3,1X.2A3,1X,2A3> 

13  i  4 11  1 F  (  IYR.GT.  IBY)GO  TO  75G 

13i45C#*«*«adDiTionAL  CONTROL  INPUT  To  DETERMINE  SOIL  STRENGTH  REQUEST** 
13;5.>  PRINT, "Do  YOU  want  CONE  INDEX  FOP  FIRST  SOIL  LAYER?  SECOND  LAYER?" 
l3l 60  PRINT, "Rcl  FOP  FIRST  LAYER?  SECOND  LAYER?  SEPARATE  ANSHEPS" 

1 3 1 7 0  PRINT. "WITH  COMMAS.  EXAHPLEiFOR  RCI  ONLY  FOR  BOTH  i.AYERS)2,2,l,l« 

13  |  READ,  <nSTE0(!).1*1.4> 

13 1 9i.  print, "Enter  names  of  soil  strength  files  to  be  usfd. enter" 

1320.,  PRINT,"  NOOATA  IF  A  FILE  IS  NOT  TO  BE  USED." 

132 lo  PRINT, "ExamplE:F0R  RCI  FOR  BOTH  lAYERS»NODATA, NODATA, XXi3oi,XX2001" 
1372.  READ  4i,  NCILA1.NCILA2.NRCLA1.NRCLA2 
13,  3-j  750  CONTINUE 

13235C.».*.TEST  OUTPUT  CHECKS  AND  OPEN  F I LES*. **••**•***. ••****«*••••••*• 

1324.  DO  /73  Id. 4 

13/5, ■  I F ( NSTEQ ( 1  ) . EO • 2 )  GO  TO  773 

13260  IF ( lYR.EO. IBY.OR.NEwFIL.EQ.1)  GO  TO  752 

1 3 2  7 v  I F  < 1  . EO . 1  ) CALL  OPENF (2.0CILA1) 

1328'  I F ( I , EO < 1 )  CALL  SETSCT ( 2 . SECTS ) 

1379c  I F ( I , EQ , 2 ) CALL  OPENF ( 3 . DC  I LA2 ) 

13»0..  I F  <  I  .  EQ ,  2  )  CALL  SETSCT  ( 3 ,  SECT6  ) 

13.11  •  I F  ( I  ,  EQ •  3  )  CALL  0PENFI2.DRCLA1) 

1332,  I F ( I  . EQ , 3 )  CALL  SETSCT ( 2. SECT7 ) 
l3?3o  1 F ( I , fcO , 4 )  CALL  OPEnF (3,DRCl*2) 

1 3 3 4 u  I F  1 1  . EO . 4  »  CALL  SETSCT (3»SECT8) 

13,00  752  XLIMi=STRLI(1.I) 

1336o  XLIM2SSTRLI (2, 1 ) 

133/'.  DO  753  J  =  K,LK 
1338o  jB  =  l 

1339o  I T  c j.fQ, IBM. AND. IYR.EQ. IBY)JB»IBD 
134  00  iF(J.EQ.NMAXD(1))  JEND  =  NMAXD(2) 

1341.)  I F  { J .  EQ ,  NMAxD  <  1 ) )  GO  TO  754 

1342.,  jEnD  =  MAXD<J) 

1343(1  754  CONTINUE 

1314.,  oO  753  LsJB.JEND 

1345.,  1 F {  I  .EQ.2.0R.  I  ,E0.4)G0  TO  755 

1346.,  iF (STORHC(J.L) .EO.O. )STO< J»L)»9. 

134/i,  |  F  (STORMCC  J,L  )  ,E0, 0. )  GO  TO  753 

134BJ  sTO ( Ji L ) =EXP <STCON< 1 , 1 ) ♦STC0N(2. t )*ALOG(STORMC<J»L> > ) 

1349.,  iF(ST0(J,L).GT.XLIM1)ST0(J»L)=XLIM1 
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Table  B7  (continued) 

Listing  of  Subroutine  CALS3T 
1-5*1 0»  iF(ST0( J.L)  .LT.XLlH2)ST0( J»L)*xUM2 

1-551*’  GO  To  753 

1352^  755  IF(STORMD(J.L> .E0.0. )ST0<J.L)=0. 
l5S.Su  iF<STORMD(J.L>  .EO.O.  >GO  TO  753 

1354  •  sTOl J.l)=EXP(STCON(l. I )*STC0N<2, 1 ) • ALOG ( STORHDl J. L ) ) ) 

13550  I F ( S  TO ( J.L) .  GT  .  XL  1  HD  STO  <  J»  L  >  3XL  1  Ml 

1356.J  jF(ST0(J,L)  .LT.XLtMZ)ST0<J.L)*XUM2 

155/u  753  CONTINUE 

135U0  I F ( HOO ( I  <  2 )  . EO , 0 )  N0  =  3 

135V.J  lf(HOD(I.2).EO.O)GO  TO  756 

136**0  nO  =  2 

136  05C  <>••••  WRITE  HE AD  I NGS*»« ••••«••*•• •••*••* •*•••»••••••••«•*••••**« 

13610  756  HRITE1N0.27)  1SIHUM 
13620  lF(l.E0.3.0R.!.E0.4)  GO  TO  757 
I3n3ij  yR  1  T£  (NO .39 ) 

1364ft  iF(I.E0.1)GO  TO  758 

13650  lF{hfcA0Ul.EQ.l)WR!T6(N0.38) 

1366u  JF (MEAOUi.E0.2)KRlTE(N0.34  ) 

136 ! U  GO  To  759 
13680  757  WP.lTE<N0.40> 

156Vu  IF { l  .E0.3JG0  TO  758 
1 3 7 0 u  iF (HEAoUl.EQ.l )HRITE(N0*38> 

1371  .  IF(HEA0U1.EQ.2)WRITE(N0.34> 

1 3 7 2 0  r,0  To  759 

1 3 7 3 (!  758  IF(HEAOUl.EQ.l)  HRITE(N0.37) 

1 3  7  4  i,  jF  (HEAOUltEO-2)wRITE(NO»33) 
l3/5(,  759  KR  I  TE  ( NO.  30  ) 

1 3 7 6 1,  DO  76?  LOUTsK.LK 
1 3  7  / 1*  jBsl 

1 3 ; 8 .j  IF(  Ipd.Gt.1. AND. LOUT. EO. IBM. AND. IYR.EQ. I8y>JB=IBD 
1 3  7 V i)  i  DAb  =  l 
1  3 8 *) .J  J  DAE:  7 

1381(  WRITE (NO .31) ( 1M0A< II .LOUT  >»11— 1.4J.JYR 

1302, i  oO  7ol  Nw  =  l»5 

13H3'  iFdYR.GT.  IBYJGO  TO  762 

1564c  iNK*0 

1385.,  ]F < Ajbn.GE.NW. AND. LOUT. EO. 1BM)G0  TO  763 
1 3 8 6 (i  jF  ( JB.  EO .  1  <0R  .  JBN.LT  . (NW-1 ) )G0  TO  762 

13865C.»»». WRITE  OUTPUT  WhEN  BEGINNING  WEEK  IS  ONLY  PARTIALLY  PILLED 

1387..  [)0  764  NKDDAB.IDAE 
1588i.  JNK*INK«1 

l3pVu  lF(STO(LOUT,NK).EO.C.O)GO  TO  764 
I390i  J ND*60» INK-1 
I391i.  jNK*NK 
13929  JNJSUOUT 
1 3 9 d sTDN:l. 

1394  ,  c*LL  WRIBWK 
13<iS  .  sTDN  =  0  . 

13960  GO  To  763 
139/.,  76-.  CONTINUE 
1398?  G0  To  763 

139Vu  762  IF (LOUT . EQ . 2 , AND « NW , EQ, 4 )NM*NW41 
1 4 r, 0 . j  iF,Nw.E0.5)  IDAEaMAxD(LOUT) 

1 4 U 1 u  jN=jend/7 
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Table  B7  (continued) 

Listing  of  Subroutine  CALSST 

l4  02u  if  (LOUT.  EO. NHAXD (1 ) . AND . JN. £0. NW )  IDAEsJEND 

l4n25C.»«*«<»WRlTfc  DATA  ON  F ILES««««»»*«»»5«*«*«»»»**»»*«**»**#»#**»***» 

I4y3,j  hR|TE(NG.32)  IOAB . I  DAE. < STOILOUT, JJ> . JJ= I  DAB. 1DAE) 

1404U  1F(L0UT.E0.NMAXD(J) .AND.JN.EO.NWJGO  TO  765 
l4i!bu  763  JDa8zI0AB*7 
14»6’i  761  10AE^IDAE*7 

14o/j  765  IF (LOUT .EQ.6.0R .LOUT.EQ. 12)  WR!TE<N0»29) ISLEW 
1 4 rj 0 1)  !F(LOUT.lT.6.0R.L0UT.GT.6)  GO  TO  760 
l4i)Vu  lFd.EO.3-OR. I. EO. 4)  GO  TO  766 
14uV5C»*««»HR1tE  hid-year  headings 
1 4 1 0 ii  hRiTE(NO,27)ISINUH 
14110  wR I T£ ( NO . 39 ) 

1 4 1 2 o  IF (I .EQ.1.AND.HEA0U1.E0.1)WRITE(N0»37> 

14i3Ci  | F  <  I  . EO. 1 . ANO.MEAOUl .EQ.2)HRITE(N0»33  ) 
l4l4u  | F < I  .E0.2.AND.MEA0U1 .E0.1)WRlTE(N0»38l 
1 4 1 1. <j  IF  (I  ,E0.2.AND.HEA0U1.E0,2)HRITE(N0.34  i 

1 4 1 6 u  GO  TO  767 
1 4 1  / 1*  766  WRITE(N0.27)!SINUM 
I4l8n  hRITE(N0,40) 

14i9u  IF ( 1 .E0.3.AND.KEA0Ul.Etl.l)WR!TE<N0»37> 

1<200  I F  ( I  .E0.3.AND.HEA0U1 . EO .2 ) WR I TE < NO « 33 ) 

14210  I F ( I .E0.4.AND.MEA0U1 . EO . 1 > HR  I TE ( NO. 38 > 

1^220  lF{l,EQ.4.AND.HEA0Ul.EO.2)WRITE(N0»34) 

1423u  767  HR  1  TE  <  NO  >  3(J ) 

14240  76u  CONTINUE 
l42btj  I F  ( NEWr  1 L  .EO .  1 )  GO  TO  768 
14260  IF  ( I  YR.GT  .  IBY  IliO  TO  769 
142/0  768  IF ( I -2 )  770.7'1.772 
1428u  7/o  CALL  SAVSCT  < 2 . SECT5 ) 

1429.)  CALL  Cl0sEF(2.DCILA1) 

14300  GO  TO  773 

14310  7  7l  CALL  SAVSCT ( 3 . SECT6  J 
1 4 3 2 o  CALL  CL0SEFC3.DCILA2) 

1433.)  GO  TO  773 

1434-  7/2  1 F ( I . GT . 3 )  GO  TO  774 
14390  c*LL  SAVSCT ( 2 . SECT7 ) 

14360  C*LL  CL0SEF(2.DRCLA1) 

1 4 3 / •»  GO  TO  773 

1438.1  774  CALL  SAVSCT  (  3  .SECT8  ) 

1439.  CALL  CL0SEF(3.DRCLA2) 

14400  G0  To  773 

1441,  769  IF(I-2)775.776.777 

14421,  775  CALL  SAVSCT ( 2 , SECTS > 

14430  CALL  CLOSEF ( 2 ) 

1 4 4 4 o  GO  TO  773 

14450  776  CALL  SAVSCT ( 3 . SECT6 ) 

l44t'i  CALL  CLOSEF  {  3  ) 

1447.  gO  To  773 

1448,,  777  I F  ( I  .GT.3)G0  TO  778 

14490  C*LL  SAVSCT ( 2 . SECT7 ) 

1«50(.  CALL  CLOSEF ( 2 ) 

14510  GO  TO  773 

l452o  778  CALL  SAvSCT (3.SECT8 ) 

1453(1  CALL  CLOSEF  ( 3  ) 
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Table  R7  (concluded) 
Listing  of  Subroutine  CALSST 

I4i>4.,  Vi  C0Nt!NUE 
l4bb.-  nC1LA1(2)*NMCLA1(2) 

1456o  nCILA2(2)sNHCLAk2) 

145/0  NRCLAl<2)=NMCLAl(2) 

14500  nRCLA2(2)=NMCLA1(2) 

1 4 5 V u  .MEwriL  =  0 

l45V5C.o*».TEST  NEED  FOR  NEW  CONTINUATION  F|LE«**»« 
14600  IF(HOD(NY1SEO»18).EO.O)NEHFIL  =1 
14610  RETURN 
1462,,  gND 


APPENDIX  C:  INPUT  DATA  ORGANIZATION  AND  FORMAT 


1.  The  format  for  preparing  control  and  weather  data  on  cards  is 
described  in  this  appendix  (see  Part  III  of  the  main  report).  Figure  Cl 
shows  the  texture  code  for  soils  of  the  U.  S.  Department  of  Agriculture 
textural  classification  system  that  is  entered  on  type  1  card  (see  fig.  C6). 
Figure  C2  shows  a  complete  card  deck  setup;  and^figs.  C3-C5  show  the  card 
deck  setup  for  each  of  three  examples  illustrating  the  use  of  specific,  sur¬ 
face  composition  group  and  estimated  (tentative  average)  data,  respectively. 

2.  Figures  C6-C19  show  the  positions  of  the  variables  on  each  of  the 
cards  with  a  brief  description  of  the  variables.  The  decimal  point  as  shown 
is  preferred;  however,  it  may  appear  anywhere  in  the  area  reserved  for  the 
variable.  It  should  be  remembered  that  the  decimal  point  as  entered  takes 
precedence  over  the  format  call  by  the  program.  Double-headed  arrows  show 
the  column  limits  of  the  area  reserved  for  the  variable.  If  there  are  no 
arrows,  the  variable  must  appear  in  that  column  only.  A  "A"  indicates  a 
space  that  must  be  left  in  the  column. 
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r.  C2.  Card  deck  setup  for  a  time-sharing  file 


NOTE:  These  data  can  be  entered  by  computer  operator  or  by  user  through 
the  teletype  (keyboard  or  punched  paper  tape^.  One  card 
corresponds  to  one  lino.  An  example  using  this  setup  is 
presented  in  Appendix  E,  sample  1, 

Fig.  C3.  Card  deck  setup  using  specific  data 
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NOTE: 

R  -  REQUIRED  IN  PROGRAM 
0  -  OPTIONAL  IN  PROGRAM 


WEATHER 

DATA 


SUPPLEMENTARY 
TO  SURFACE 
COMPOSITION 
CROUP  CONTROL 
DATA  SELECTED 
FROM  FILE 
"DSURGR" 


An  example  using  this  setup  is  presented  in  Appendix  E,  Example  2. 


Fig.  C4 .  Card  deck  setup  using  surface 
composition  group  data 
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THESE  IS  — L 
REQUIRED 


NOTE: 

R-REQUIRED  IN  PROGRAM 
O-OPTIONAL  IN  PROGRAM 

An  example  u3ing  this  setup 


is  presented  in  Appendix  E 


CONTROL  DATA  AND 
AREA  IDENTIFICATION 


example  3. 


Fig.  C5.  Card  deck  setup  using  estimated 
(tentative  average)  data 
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Fig.  C13.  Type  3  Card  -  Beginning  Soli  Moisture  Content 
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Uumber  of  the  last  coefficient  entered  as  prescribed  by  the 
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Coefficients  of  other  appropriate  soil  property  terms  in 
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APPENDIX  D:  STORED  INPUT  DATA  FILE  LISTINGS 


1.  This  appendix  includes  listings  of  two  files  containing  control 
information  for  the  surface  composition  groups  and  tentative  average  soil 
strength  relations,  respectively.  Table  Dl,  file  "DSURGR",  contains  the 
required  control  data  input  to  the  program  for  each  of  the  surface 
composition  groups.  The  first  set  of  figures  on  each  line  is  the  line 
numher.  The  remainder  of  the  line  follows  the  organization  of  input  data 
as  described  in  Appendix  C.  For  example,  on  the  first  line,  501,  the 
first  number,  1,  identifies  card  type  1  with  the  following  information: 
the  surface  composition  group  is  0505;  the  two  "C's"  indicate  a  storage 
location  within  the  program  for  holding  the  depletion  equation  coef¬ 
ficients  for  the  two  soil  layers.  The  remaining  two  values,  1.50  and  1.50, 
are  the  soil  densities  for  the  two  layers. 

2.  Table  D2,  file  DFSTEQ,  contains  data  for  the  variables  in  the 
tentative  average  soil  strength  relations.  The  information  to  the  left 
of  the  line  number  is  not  in  the  file.  It  is  included  for  identification 
of  the  information  only  and  matches  the  number  of  the  equation  set  as 
defined  in  tables  6  and  7  of  the  main  text.  The  numbers  on  the  first  line 
of  each  set  to  the  right  of  the  line  number  identify  the  form  for  each 

of  the  variables  required  in  the  two  equations  of  the  set.  For  example, 
equation  set  4  for  RCI  (table  7,  main  text)  provides  data  on  variables 
for  liquid  limit  and  plastic  limit,  respectively.  (The  data,  soil 
property  values,  are  input  on  card  6,  see  fig*  Cl6a,  Appendix  C.)  On 
line  370  of  the  file  a  3  in  the  fifth  slot  indicates  the  In  In  form  for 
liquid  limit  and  a  1  in  the  sixth  slot  indicates  the  natural  form  for 
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plastic  limit.  Zeros  for  the  other  variables  indicate  that  they  are  not 
used  in  the  equation.  The  appropriate  sets  of  equations  are  identified 
by  numbers  1-6  entered  on  card  2-5  illustrated  in  fig.  Cll,  Appendix  C. 

Line  380  of  file  DFSTEQ  contains  the  constants  of  the  equation  for  calculat¬ 
ing  the  moisture  content  of  a  soil  with  an  RCI  of  100  when  the  liquid  and 
plastic  limits  are  known.  In  this  case,  the  equation  is: 

A 

MC  (of  a  soil  with  an  RCI  of  100)  =  0.783718  +  1.498687  (lnlnLL)  +  0.026353(PL) 


ii 

5 
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501  1  0505  C  C  1.50  1.50 

502  2  1  .1  0.0/  200.  75.  0.  0.  0. 

50/  ?  1  11  0.110707-0-0231589  -0.105756  0. 141049 
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APPENDIX  E:  EXAMPLES  OF  PREDICTION  RUNS 


1.  The  purpose  of  the  examples  in  this  appendix  is  to  illustrate 
the  input  data  setup  required  for  prediction  of  soil  moisture  content 
and  soil  strength  and  for  input  to  the  Free2e-Thaw  model.  Example  1 
illustrates  a  prediction  run  using  specific  (measured)  data  and  relations 
derived  from  measured  data  and  example  2  shows  the  setup  for  a  run  using 
the  surface  composition  group  data.  In  example  3  the  only  soil  and 
terrain  data  available  are  those  that  have  been  derived  from  estimated 

or  averaged  values  and  from  weather  records  from  a  nearby  weather  station. 
Ir.  this  case  tentative  average  relations,  which  are  built  into  the  program, 
are  used.  Example  4  shows  input  and  output  data  required  for  use  in  the 
Freeze-Thaw  model. 

2.  Appendix  F  contains  graphic  displays  of  the  results  of  runs  for 
the  first  three  examples. 


Example  1,  Run  Using  Specific  Data 


3.  Assume  that  daily  prediction  of  soil  moisture  content  in  inches  and 
percent  and  Cl  and  RCI  for  the  two  layers  is  needed  at  a  site  for  the 
period  13  April  1951  -  31  December  1952.  Records  are  available  for  data 
as  follows: 


a.  Physical  soil  properties 


X  Organic 

X  Sand  X  Silt  Z  Clay  Matter 

Layer  15  78  17  1.71 

Layer  24  78  18  1.24 


Dry 

Density,  g/cc 
1.35 
1.53 


El 
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b.  The  site  is  located  on  the  levee  of  a  bottomland,  which  is 
classified  as  wetness  index  4  (table  5,  main  text). 

£.  The  soil  is  a  Collins  silt  loam. 

d.  The  following  relations  have  been  determined  from  measured 
data: 


l 

5 

i 


(1)  Accretion  relations 

Layer  1 
Class  I  0.7X-0 .04 
Class  II  0.96Z-0.06 


Layer  2 
0.44X-0.10 
0.76Z-0 .02 


Minimum  storm  -  0.10  in.  Layer  1-Class  I  relation 


changes  to  a  1:1  relation  when  the  soil  moisture 


content  prior  to  accretion  is  0.10  in. 


(2) 

X 

Depletion 

x2 

relations 

X3 

X4 

X5 

X6 

Summer 

Layer  1 

0.112577 

0.892538 

-2.41564 

2.34487 

-1.01126 

0.163761 

Layer  2 

0.0875938 

0.148577 

-0.303785 

0.127866 

0 

0 

Transition 

Layer  1 

0.233207 

-0.454052 

0.359594 

-0.0938728 

0 

0 

Layer  2 

0.404195 

-2.2826 

5.68009 

-6.27261 

2.52196 

0 

Winter 

Layer  1 

0.198281 

-1.61755 

7.46858 

-10.673 

0 

0 

Layer  2 

1.01464 

-24.9913 

255.935 

-1102.69 

1693.1 

0 
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(3)  Strength  relations 


_ Cl _  _ RCI _ 

Layer  2  InCI  -  10.72  lnMC-1.64  InRCI  »  23.864  lnMC-5.815 

e.  Limits  of  soli  moisture  content 

Maximum  _ Minimum _ 

_ By  Seasons _ 

Summer  Transition  Winter 

Layer  1  2.52  in.  0,66  in.  1.02  in.  2.13  in. 

Layer  2  2.43  in.  0.90  in.  1.53  ic,  2.19  in. 

f.  Measured  soil  moisture  contents  on  13  April  1951  were 

2.12  in.  for  the  first  layer  and  2.26  for  the  second  layer. 

Beginning  of  season  dates 

Spring  Summer  Fall  Winter 

10  Mar  10  May  21  Sept  24  Oct 

h.  Precipitation  as  listed  in  data  file  (table  El)  beginning  with 
line  20. 

4.  From  this  information  the  required  control  data  (lines  1-13  of 
table  El)  and  weather  records  (beginning  line  19  of  table  El)  are  set  up 
for  a  daily  soil  moisture-scil  strength  prediction.  A  card  deck  would  be  set 
up  in  accordance  with  fig.  C3,  Appendix  C.  The  teletype  irput  during  the 
running  of  the  program  is  shown  in  table  E2.  Partial  lists  of  predicted 
soil-moisture  contents  in  inches  and  percent  and  Cl  for  the  6-  to  12-in. 

(15-  to  30-cm)  layer  are  shown  in  tables  E3-E5. 


Example  2.  Run  Using  Surface  Composition  Group  Data 


5.  Suppose  all  that  is  known  about  the  area  in  example  1  is  that: 
(a)  it  is  on  the  levee  of  a  bottomland;  (b)  the  water  table  occurs  at 


E3 
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the  surface  less  than  90  percent  of  the  time  and  above  the  120-cm 
depth  more  than  10  percent  of  the  time;  (c)  the  soil  is  silt  with  an 
organic  content  less  than  7  percent;  and  (d)  the  starting  moisture  con¬ 
tents  for  the  two  layers  and  the  weather  data  (daily  precipitation  and 
beginning  of  season  dates)  are  known.  Daily  predictions  of  soil  mois¬ 
ture  content  and  soil  strength  can  be  made  by  using  the  program  controls  from 
file  DSURGR  (table  Dl,  Appendix  D) .  To  determine  the  surface  composition  group 
for  this  area,  a  search  is  made  of  table  1  of  the  main  text,  and  Group 
0808  is  selected  for  obtaining  control  data.  The  card  deck  setup  is  shown  in 
fig.  C4,  Appendix  C.  Table  E6  lists  the  input  data  file  printout  for 
Group  0808,  and  table  E7  shows  the  teletype  input  during  the  run.  Partial 
outputs  of  soil  moisture  content  in  inches  and  percent  and  Cl  are  shown 
in  tables  E8-E10. 


Run  Using  Estimated  Data 


6.  Assume  that  all  that  is  known  about  the  area  in  example  1  is 
that  it  is  a  bottomland  with  an  estimates  wetness  index  of  4.  Percent¬ 
ages  of  sand,  silt,  clay,  and  organic  matter,  and  values  of  dry  density 
are  obtained  from  soil  samples.  U.'.th  this  information  and  weather  data, 
an  input  deck  is  set  up  and  a  data  file  is  output  as  shown  in  table  Ell. 

The  program  then  computes  maximum  and  minimum  soil  moisture  contents. 

The  maximum  moisture  content  is  used  as  beginning  moisture  cr— .tent.  The 
program  proceeds  with  daily  predictions  by  using  tentative  average  relations 
for  accretion  and  depletion.  (See  fig.  C5,  Appendix  C,  for  the  deck  agtpn.) 
The  teletype  input  and  response  during  run  time  is  shown  in  table  E12. 
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Tables  E13-E15  show  partial  outputs  of  soil  moisture  content  in  inches 
and  percent  and  Cl  predictions. 

Example  4,  Run  Providing  Data  for  Freeze-Thaw  Model 

7.  An  example  of  a  data  file  with  air  temperature  and  snow  depth 
input  data  is  shown  in  table  E16.  The  weather  data  setup  is  shown  in 
fins.  C2  and  C3,  Appendix  C.  The  teletype  input  and  response  for  the 
Freeze-Thaw  model  are  presented  in  table  E17 .  Table  E18  is  a  listing 
of  the  required  soil  moisture  content  and  weather  data  output  by  the 
program. 
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Table  El 


Example  1.  Data  File  Printout 


00001 

1 

RIFLE 

TS0002  MM  1 

.35  1.53 

0  0  U  0  2 

2 

1 

O.lOj  0.100  750. 

0.  750.  0.  300. 

0 

.  300. 

0. 

00003 

2 

2 

11 

.112577  .892538 

-2.41564  2.34487 

-1- 

01126 

163761 

00004 

2 

2 

12 

.1875938  .148577 

-.303785  .127866 

o. 

0. 

O0UO5 

2 

2 

21 

.233207  -.454052 

.359594  -.0938728 

o. 

0. 

00006 

2 

2 

22 

.404195  -2*2626 

5.68009  -6.27261 

2- 

52196 

0. 

00007 

2 

2 

31 

.198281  -1.61755 

7.46858  -10.673 

0. 

0. 

000n8 

2 

2 

32 

1.01464  -24.9913 

255.935  -1102,69 

1693.1 

0. 

00009 

2 

3 

1.86  1.53  1.50 

0.90  0.39  0,24 

0  0  0 1 0 

2 

4 

■,.0  0.0  10.72 

-1.64  0.0 

0.0 

23.864 

-5.015 

0  0 1 1 

2 

6 

0.700  0.960  0.440 

0.760  -0.04 

-0.06 

-0.10 

-0.02 

00012 

3 

2 . 

12  2.26 

00013 

4 

2. 

52  2.43  0 

.66  0.90 

00014 

7 

SPECIFIC  RELATIONS 

00015 

7 

LOCATION  -  RIFLE 

RANGE-  VICKSBURG.  MISS. 

0  0  y  1 6 

7 

SOIL  -  COLLINS  SILT  LOAM 

00017 

7 

VEGETATION  -  HERBACEOUS 

00016 

7 

RECORD  -  13 

ApR  1951-31  DEC  1952 

0  0  0 1 9 

62 

13 

4  51 

00020 

9 

8 

451 

77.77  77-77  77-77 

77-77  77-77 

.02 

00021 

915 

451 

.63 

2.69 

00^22 
00023 
00024 
U0|!25 
00026 
00»1<>  7 
00028 
o  o  n  i?  9 
0  0  (l  3  0 
00i)3l 
00(132 
00(133 
00034 
00035 
00Q36 
0  0  ii  J  7 
00.138 
001)39 
0  0  (1 4  ij 
0  0  (s  4  l 
00042 
00043 

00i)44 


922  45i 
9 29  45j 
9  1  55  l 
81  i.  5 
9  8  551 
915  55.1 
522  55: 
929  55: 

9  1  65i 
9  8  651 
915  65i 
922  65: 
929  651 


.0? 


5l 


.03 


9  1 
9  8 
915 


751 

751 

751 


922  751 
929  751 


9 
9 
915 


1  851 
8  851 
851 


.40 

.81 


.28 


,  4 ; 


.43 
.4  / 

.24 

.72 


.05 


.22 


1.79 


.40 


.06 


.01 


,02 

.13 


.01 


,06 


.64 

.88 


2.66 


1 .61 


922  851 
929  851 


2Z6 


Table  E2 


Exaaple  1.  Teletype  Input 


FOR  INFORMATION  ON  RUNNING  THIS  MODEL  SEE  USER'S  MANUAL. 

WHEN  ANSWERING  YES  OR  NO  QUESTIONS. 

USE  1  FOR  YES  AND  2  FOR  NO. 

IS  INPUT  RAINFALL  IN  CM’ 

INPUT*  01 560 
’  2 

IF  YOU  WANT  OUTPUT  MOISTURE  CONTENT  IN* 

INCHES  ENTER  A  H  CM  ENTER  A  2 
INPUT*0 1 590 
’  1 

NAME  LAYERS  1  AND  2  SOIL  MOISTURE  CONTENT  FILES. 

EXAMPLE  *  XX 1 0  0 1 . XX200 1 
INPUT*0 1620 
’  AMI 001 » AM2001 

DO  YOU  NEED  OUTPUT  MOISTURE  CONTENT  IN  PERCENT’ 
INPUT*01640 
’  1 

NAME  LAYERS  1  AND  2  PERCENT  MOISTURE  FILES. 

EXAMPLE  *XX1 00 1 .XX200 1 
INPUT*0 1680 
’  AP1001 . AP2001 

IF  YOU  HAVE  AIR  TEMPERATURE.  ENTER  A  TJ  IF  NOT  A  9 
INPUT*  0 1710 
’  9 

IF  YOU  HAUE  SNOW  DEPTH CM). ENTER  AN  SN)  IF  NOT  A  9 
INPUT*  0 1 7A0 
’  9 

WHAT  IS  YOUR  INPUT  SOURCE’ 

ENTER*  1  FOR  SURFACE  COMPLEX  GROUPS.  OTHERWISE  A  2 
INPUT*01780 
’  2 
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Table  E2  (concluded) 


ENTER  THE  NAME  OF  YOUR  INPUT  FILE 
INFUTSO 1810 
’  UARIFH 


ARE  MAX  AND  MIN  SOIL  MOISTURE  EQUATIONS  INCLUDED  IN  THE 
INPUT  DATA  FILE’ 

INPUT*  0 1 840 
’  2 


SITE  RIFLE  TS0002 


TEXTURE 
DENSITY 
BEG  MC 


H  M 

1.35 

2.12 


1.53 

2.26 


MAX-MIN  MC  2.520  2.430  0.660  0.900 

SPECIFIC  RELATIONS 

LOCATION  -  RIFLE  RANGE.  VICKSBURG.  MISS. 

SOIL  -  COLLINS  SILT  LOAM 
VEGETATION  -  HERBACEOUS 
RECORD  -  13  APR  1951-31  DEC  1952 
DO  YOU  WANT  CONE  INDEX  FOR  FIRST  SOIL  LAYER’  SECOND  LAYER’ 
RCi  FOR  FIRST  LAYER’  SECOND  LAYER’  SEPARATE  ANSWERS 
WITH  COMMAS.  EXAMPLESFOR  PCI  ONLY  FOR  BOTH  LAYERS? 2.2. 1 .  1 
INPUTS  13180 


’  1.1. 1.1 


ENTER  NAMES  OF  SOIL  STRENGTH  FILES  TO  BE  USED. ENTER 
NODATA  IF  A  FILE  IS  NOT  TO  BE  USED. 

EXAMPLESFOR  RCI  FOR  BOTH  LAYERS? NODATA. NO DATA. XX 1 00 1 .XX200 1 
INPUTS  13220 

’  AC1001.AC2001.AR1001.AR2001 


STOP 


RUNNING  TIME? 


78*1  SECS  I/O  TIME  s  42. 6  SECS 
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Example  1.  Predicted  Soil  Moisture  Content  In  Inches 


SPECIFIC  RELATIONS 

LOCATION  -  RIFLE  ,.ANGE.  VICKSBURG,  MISS. 
SOIL  -  COLLINS  SILT  LOAM 
VEGETATION  -  HERBACEOUS 
RECORD  -  13  APR  1951-31  DEC  1952 


C*Y 

8-14 

15-21 

2.17 

2?-2b 

2.33 

2.12 

cay 

1-  / 

2. 16 

8-14 

1.88 

15-21 

1.69 

22-28 

1.30 

29-31 

l.f  8 

CAY 

1-  / 

1.-  3 

8-14 

1.96 

15-21 

1.67 

22-28 

1.62 

29-30 

1.  70 

SOIL  MOISTURE  IN/6 
6  -  12  IN  LAYER 


APRIL 


2.14 

2.11 

2.08 

2.31 

2.28 

2.25 

2.09 

MAY 

2.03 

2.01 

1.98 

1.96 

1.99 

1 . 93 

1.63 

1.57 

1.53 

1.26 

1.22 

1.16 

1.06 

j  .04 

JUNE 

1.02 

1.00 

0.99 

1.05 

1.74 

1.68 

1.78 

1.72 

1.79 

1.56 

1.50 

1.44 

1.71 

N 

1951 


2.23 

2.20 

2.25 

2.22 

2.36 

2.22 

2.19 

2.16 

1951 

1.95 

1.93 

1.90 

1.87 

1.81 

1-75 

1.45 

1.40 

1.34 

1.15 

1.13 

1>10 

1951 

0.98 

0.98 

0-97 

3  .62 

1.55 

1.74 

1.80 

1.74 

1-68 

1.39 

1.34 

1.63 
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Predicted  Soil  Moisture  Content  In  Percent 


IDENTIFICATION  NUMBER:  RIFLE  TS0002 
PERCENT  moisture  content 

6  -  12  IN  LAYER 


D*Y 

JANUARY 

1952 

1-  / 

24.18 

24.03 

23.91 

25.64 

25.35 

25.12 

24.94 

0-14 

24.78 

24.80 

24.66 

24 . 52 

24.30 

24.2? 

24.07 

15-21 

23.94 

23.87 

23.86 

23.66 

23.86 

24  .59 

24.45 

22'20 

24.30 

24.14 

24.00 

23.89 

23.86 

25.63 

25*34 

29-31 

25.11 

24.93 

24.78 

day 

FEBRUARY 

1952 

1-  / 

25.  es 

25.54 

25.27 

25.38 

25.14 

24.95 

24.80 

0-14 

24.66 

24.52 

24.38 

24.2  2 

25.71 

25.42 

25.17 

33-21 

25.94 

26.27 

25.96 

25.65 

25.36 

25.99 

25-68 

22-25 

26.  t'6 

25.76 

25.46 

26. Cl 

26.34 

26.02 

25.71 

PAY 

MARCH 

1952 

1-  / 

25.42 

25.17 

25.94 

25.63 

25.35 

25,12 

24-93 

0-14 

24.70 

24.00 

25.87 

25.50 

25.28 

24 , 97 

24.66 

15-21 

24 ,  *3 

24.00 

23.66 

25.58 

25.26 

24.97 

24  66 

22-20 

25.  bo 

25.52 

25.22 

24.91 

24.60 

24.27 

23-94 

29-31 

23.60 

23.26 

22.92 

r*Y 

APRIL 

1952 

1-  / 

22.59 

22.27 

21.96 

23.26 

22.93 

22.60 

22*27 

0-14 

21.96 

21.66 

21.77. 

21-47 

24.22 

23.88 

23-54 

15-21 

23.20 

22.87 

22.54 

22.2? 

21.91 

22.02 

21-71 

22-20 

21.42 

25.04 

24.73 

24.41 

24.08 

23.74 

23*40 

29-31' 

23.  '16 

22.73 

DAY 

may 

1952 

1-  / 

22.40 

24.86 

24.54 

24.21 

23.86 

23.54 

23*20 

$-14 

22.87 

22.54 

23.27 

22.61 

21.97 

21.33 

20.68 

15-?1 

20.03 

19.36 

18.69 

23.78 

25.61 

24.75 

24-01 

22*20 

24.11 

25.69 

26.06 

25.13 

24.34 

23.63 

25-57 

29-31 

26.o4 

25.10 

24.3? 

C*Y 

JUNE 

1952 

1-  • 

23.61 

22.95 

22.30 

21.66 

21.01 

20.36 

19.70 

ft-14 

l9.:-3 

18.36 

17.68 

17.01 

16.37 

15.74 

15-16 

15-21 

14.61 

14.10 

13.64 

13.23 

12.86 

12.53 

12*24 

22-20 

11.97 

11.74 

11.54 

11.36 

11.20 

11.06 

10*93 

29-30 

10.82 

10.72 

230 


Table  E5 


I  DENT  1 F 1  CAT  I  ON  NUMBER  •  R1TLE  TS0Q02 
CONE  INDEX 
6  -  12  IN  LAYER 


cay 

JANUARY 

1952 

3  '  ' 

244. 

247. 

249. 

222. 

226. 

229. 

232 

H-14 

234. 

234. 

236. 

239. 

241. 

243. 

246 

15-21 

248. 

249. 

250. 

250. 

250. 

238. 

24Q 

22- 2» 

242. 

245. 

247. 

249. 

250. 

222. 

226 

29-31 

2  5;,. 

2  52. 

235. 

cay 

FEBRUARY 

1952 

l-  7 

2l9. 

223. 

227. 

226. 

229. 

232. 

234 

0-1 A 

236. 

239. 

241. 

243. 

221. 

225. 

229 

15-21 

2lU. 

213. 

217. 

22  2. 

?26. 

217. 

22l 

22-28 

2l6. 

220. 

224. 

217. 

212. 

217. 

221 

i:  ay 

MARCH 

1952 

232 

•  -  / 

225. 

229. 

218. 

222. 

226. 

229. 

6-14 

235. 

233. 

219. 

223. 

227. 

232. 

236 

15-21 

24  2. 

247. 

253. 

223. 

227. 

232. 

236 

22*28 

219. 

224. 

228. 

233. 

237. 

243. 

248 

29-31 

254. 

260. 

266. 

C  *  Y 

APRIL 

1952 

v-  / 

273. 

279. 

286. 

260. 

266. 

273. 

279 

8-14 

286. 

292. 

290. 

297. 

244. 

249. 

255 

15-21 

261. 

267. 

274. 

280. 

287. 

285. 

29i 

2?-2« 

296. 

231. 

235. 

240. 

246. 

252. 

258 

29-30 

264. 

270. 

D*Y 

MAY 

1952 

1-  / 

277. 

233. 

238. 

244. 

249. 

255. 

26i 

8-14 

268. 

274. 

260. 

272. 

286. 

300. 

315 

15-21 

332. 

351. 

372. 

251. 

222. 

235. 

247 

22-28 

245. 

221. 

216. 

229. 

242. 

254. 

223 

29-31 

2l6. 

230. 

242. 

cay 

JUNE 

1952 

1-  ’ 

254. 

266. 

279. 

292. 

307. 

323. 

34l 

8-14 

361. 

383. 

408. 

434. 

463. 

493. 

525 

15-21 

557. 

590. 

623. 

656. 

687. 

717. 

745 

2?-28 

750. 

750. 

750. 

750. 

750. 

750. 

75o 

29-30 

75c- . 

750. 

£31 


Table  E6 


Example  2.  Data  File  Printout 


00012 

3  2.12 

2.26 

00014 

7 

SURFACE  COMPOSITION  GROUP  0808 

00015 

7 

LOCATION  -  RIFLE  RANGE.  VICKSBURG.  MISS. 

00016 

7  SILT.  FINE  GRAINED.  ORGANIC  MATTER  0-7  PERcSN 

00017 

7 

INTERNAL  DRAlNiGE-CLASS  1 

00016 

7 

RECORD  -  13  APR  1951-31  DEC  1952 

00019 

62  13  4 

51 

77.77  77.77  77.77  77.77  .02 

00020 

9  8  451 

77.77 

00021 

915  451 

.63  2.69 

00022 

922  451 

00023 

929  45l 

.02 

.03 

00024 

9  1  551 

00025 

81  10  5 

51 

00026 

9  6  551 

.47  .22 

00027 

915  551 

.01 

00026 

922  551 

00029 

929  551 

OOO-SO 

9  1  651 

.43  1.79  .02  .64 

000-51 

9  8  651 

000-32 

915  651 

.47  .*0  .13 

00033 

922  651 

.  88 

000-54 

929  65l 

.40 

.24 

00035 

9  1  751 

.81 

.72 

00036 

9  8  751 

00037 

915  751 

.05  .06  .01  2.68 

00036 

922  751 

.28 

00039 

929  751 

000*0 

9  1  851 

000*1 

9  8  851 

.06  1.61 

000*2 

915  851 

000*3 

922  851 

000*4 

929  851 

.58  .03 

000*5 

9  1  951 

000*6 

9  6  951 

.71  .07  2,14 

000*7 

82  21  9 

51 

.03 

000*8 

915  951 

000*9 

922  951 

.  89 

3.82 

00050 

929  951 

.04 

00051 

9  11051 

00052 

9  81051 

00053 

9151051 

00054 

83  24  10 

51 

00055 

9221051 

.84 

00056 

9291051 

.01  .33 
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FOR  INFORMATION  ON  RUNNING  THIS  MODEL  SEE  USER’S  MANUAL. 

WHEN  ANSWERING  YES  OR  NO  QUESTIONS/ 

USE  1  FOR  YES  AND  2  FOR  NO. 

IS  INPUT  RAINFALL  IN  CM? 

INPUT: 01 562 
?  2 

IF  YOU  WANT  OUTPUT  MOISTURE  CONTENT  IN: 

INCHES  ENTER  A  1J  CM  ENTER  A  2 

INPUT:01 590 
?  1 

NAME  LAYERS  1  AND  2  SOIL  MOISTURE  CONTENT  FILES. 

EXAMPLE : XXI 00 1 / XX200 1 

INPUT :0 1620 
?  MC1001.MC2001 

DO  YOU  NEED  OUTPUT  MOISTURE  CONTENT  IN  PERCENT? 

INPUT J01 640 
?  1 

NAME  LAYERS  1  AND  2  PERCENT  MOISTURE  FILES. 

EXAMPLE: XXI 00 1/XX2001 

INPUT:0I680 
?  PM 1001/ PM200 1 

IF  YOU  HA,,C:  AIR  TEMPERATURE/  ENTER  A  T>  IF  NOT  A  9 

INPUT : 0 1 7 1 0 
?  9 

IF  YOU  HAVE  SNOW  DEPTH <M )/ ENTER  AN  SNI  IF  NOT  A  9 

INPUT: 01740 
?  9 

WHAT  IS  YOUR  INPUT  SOURCE? 

ENTER:  1  FOR  SURFACE  COMPOSITION  GROUP/  OTHERWISE  A  2 

INPUT: 01 780 
?  1 
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Table  E7  (concluded) 


ENTER  THE  NUMBER  OF  YOUR  SURFACE  GROUP 
INPUTJ0 1 890 
?  0808 


SITE 
TEXTURE 
DENSITY 
MAX -MIN  MC 
ENTER  NAME  OF 
INPUT *02770 


0808 
C  C 

1.34  1.41 

2.900 
RAINFALL  FILE 


2.800 


7  DA080K 


0.650  0.700 


BEG  MC  2.12  2.26 

SURFACE  COMPOSITION  GROUP  0806 
LOCATION  -  RIFLE  RANGE.  VICKSBURG.  MISS. 

SILT.  FiNE  GRAINED.  ORGANIC  MATTER  0-7  PERCENT 
INTERNAL  DRAINAGE-CLASS  I 
RECORD  -  1J  APR  19S1-31  DEC  1952 
DO  YOU  WANT  CONE  INDEX  FOR  FIRST  SOIL  LAYER’  SECOND  LAYER’ 
RC1  FOR  FIRST  LAYER’  SECOND  LAYER’  SEPARATE  ANSWERS 
WITH  COMMAS.  EXAMPLE* FOR  RCI  ONLY  FOR  BOTH  LAYERS; 2.2. 1 , 1 
INPUT* 1 31 RO 
’  1.1.  1.  1 


ENTER  NAMES  OF  SOIL  STRENGTH  FILES  TO  BE  USED. ENTER 
NODATA  IF  A  FILE  IS  NOT  TO  BE  USED. 

EXAMPLE* FOR  RCI  FOR  BOTH  LAYERS I NO DATA. NODATA. XXI 001 .XX20G1 
INPUT* 13220 

’  AC 1 00 1 . AC200 1 . AR1 CO  1 . AR200 1 


STOP 


RUNNING  TIME*  81.5  SE"S  I/O  TIME  *  «2.6  SECS 
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Table  E8 


Exagple  2,  Predicted  Soli  Moisture  Content  In  Inches 


SURFACE  COMPOSITION  GROUP  0808 
LOCATION  -  RIFLE  RANGE.  VICKSBURG.  MISS. 

silt,  fine  grained,  organic  matter  0-7  percent 

INTERNAL  DRAINAGE-CLASS  1 
RECORD  -  13  APR  1951-31  DEC  1952 


SOIL 

MOISTURE 

IN/6 

IN 

6 

-  12  IN 

layer 

C  A  Y 

APRIL 

1951 

?-i« 

2.24 

2.22 

l'»“2i 

2.21 

2.19 

2.17 

2.15 

2.28 

2.27 

2.57 

2  -2d 

2.55 

2.53 

2.52 

2.50 

2.48 

2.46 

2.45 

U  vi  —  3.! 

2.43 

2.41 

CAY 

MAY 

1951 

1  -  / 

2.40 

2.38 

2.36 

2.34 

2.33 

2.31 

2.29 

a-i* 

2.27 

2.37 

2.41 

2.37 

2.33 

2.30 

2.26 

15-21 

2.22 

2.18 

2.15 

2.11 

2.07 

2.03 

2*00 

2^-2b 

1.96 

1.92 

1.88 

1.84 

1.80 

1.76 

1*73 

2«,-31 

1.69 

1.65 

1.61 

TAy 

JUNE 

1951 

1-  / 

1.57 

1.53 

1.49 

1.45 

1.41 

1.38 

1-34 

8-14 

1.30 

1.39 

1.77 

1-73 

1.69 

1.65 

1.78 

1'.  -21 

1.74 

1.84 

1.80 

1,88 

1.90 

1.86 

1-82 

22-2C 

1.78 

1.74 

1.70 

1.66 

1.62 

1.58- 

1-77 

29-3'.1 

1.84 

1 . 89 
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Table  E9 


Example  2,  Predicted  Soli  Moisture  Content  In  Percent 


IDENTIFICATION  NUMBER: 
PERCENT  MOISTURE  CONTENT 
6  -  12  IN  LAYER 


0808 


LAY 

JANUARY 

1952 

!  -  - 

32.33 

32.22 

32.11 

32.47 

32.36 

32.26 

32-15 

8-1  -i 

32.05 

32.53 

32.42 

32.32 

32.21 

32.11 

32.00 

15-21 

31.90 

31.79 

31.69 

31.58 

31.48 

32.21 

32.11 

22-28 

32.00 

31.90 

31.79 

31.69 

31.59 

32.26 

32*  15 

29-31 

32.85 

31.94 

31.84 

day 

FEBRUARY 

1952 

t-  / 

32.36 

32.25 

32.15 

32.52 

32.42 

32.31 

32.21 

6-1‘t 

32.10 

32.00 

31.89 

31.79 

32.34 

32.23 

32.13 

15-2i 

32.48 

32.83 

32.73 

32.63 

32.52 

32.88 

32.78 

22-28 

33.10 

33.00 

32.90 

33.10 

33,10 

33.00 

32.90 

CAY 

march 

1952 

,  -  / 

32.80 

32.69 

33.05 

32.95 

32.05 

32.74 

32.64 

8, -1-1 

32.54 

32.73 

33.08 

32.88 

32.68 

32.47 

32.27 

15-21 

32. C7 

31.87 

31.67 

32.29 

32.09 

31.89 

31>68 

22-28 

32.30 

32.09 

31.89 

31,69 

31.49 

31.29 

31.09 

29-31 

30.89 

30.69 

30.49 

C*Y 

APRIL 

1952 

i  -  / 

30.29 

30.09 

29.89 

31.07 

30.87 

30.67 

30*47 

t-14 

30.27 

30.07 

30.60 

30.40 

31.78 

31.58 

31.38 

1'-  *21 

31.18 

30.98 

30.78 

30.58 

30  >  36 

30.60 

30-40 

22-28 

30.20 

31.73 

31.50 

31.30 

31.10 

30.90 

30.70 

2^-38 

30.50 

30.30 

r  ay 

may 

1952 

i  -  / 

30.10 

31.66 

31.46 

31,26 

31.06 

30.86 

30.66 

e-i4 

30.46 

30.26 

31.13 

30,69 

30.25 

29.82 

29.38 

15-21 

28.95 

28.51 

29.08 

30,85 

31.96 

31.52 

31.08 

2?-28 

31.30 

32.14 

32.50 

32.05 

31.61 

31.17 

32-09 

29-31 

32.46 

32.01 

31.57 

cay 

JUNE 

1952 

3-  ' 

31.13 

30.69 

30.26 

29.62 

29.39 

28.95 

28-52 

e- 14 

28.08 

27.65 

27.21 

26.78 

26.34 

25.90 

25-46 

15-21 

25 . 02 

24.57 

24.12 

23.68 

23.22 

22.77 

22*32 

22-28 

21.86 

21.40 

20.94 

20.47 

20.01 

19.55 

19*08 

29-3C 

18,62 

18,16 
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Table  E10 


Predicted  Cone  Index 


IDENTIFICATION  NUMBER:  0808 

CONE  INDEX 
6  -  12  IN  LAY£R 


C*Y 

JANUARY 

1952 

1-  ' 

119. 

120. 

121. 

118. 

119. 

120. 

121. 

8-1-4 

122. 

117. 

118. 

119. 

l2Q . 

121. 

122. 

15-21 

123. 

124. 

125. 

126. 

128. 

120. 

121. 

22-28 

122. 

123. 

124. 

125. 

126. 

120. 

121. 

29-31 

122. 

123. 

124. 

L'Ay 

FEBRUARY 

1952 

t  -  / 

119. 

120. 

121. 

117. 

118. 

119. 

120. 

8-1A 

121. 

122. 

123. 

124. 

119. 

120. 

121. 

15-21 

Il8. 

115. 

116. 

116. 

117. 

114, 

115. 

22-28 

112. 

113. 

114. 

112. 

112. 

113. 

H4. 

T*Y 

MARCH 

1952 

i-  / 

115. 

116. 

113. 

114. 

115. 

115. 

H6. 

e  -i  -4 

117. 

116. 

112. 

114. 

116. 

118. 

120. 

*5-21 

122. 

124. 

126. 

120. 

122. 

123. 

125. 

22-28 

120. 

121. 

123. 

125. 

127. 

129. 

132. 

29-31 

134. 

136. 

138. 

day 

APRIL 

1952 

i  -  / 

141. 

143. 

145. 

132. 

134. 

136. 

139. 

8-1-) 

141. 

143. 

137. 

139. 

124. 

127. 

129. 

15-21 

131. 

133. 

135. 

137. 

140. 

137. 

139. 

22-28 

142. 

125. 

127. 

129. 

132. 

134. 

136. 

2  9  -  3  v1 

138. 

140. 

C*Y 

MAY 

1952 

1-  / 

143. 

126. 

128. 

130. 

132. 

134. 

136. 

8-14 

139. 

141. 

131. 

136. 

141. 

146. 

152. 

15-21 

158. 

164. 

170. 

134. 

123. 

127. 

132. 

2  2-28 

129. 

121. 

118. 

122. 

126. 

131. 

121, 

29-31 

Il8. 

122. 

127. 

D*Y 

JUNE 

1952 

;  -  / 

131. 

136. 

141. 

146. 

152. 

158. 

164. 

e-i-i 

17U. 

177 . 

184. 

192. 

200. 

209. 

218. 

15-21 

228. 

238. 

250. 

262, 

275. 

269. 

304, 

2r-28 

321. 

338. 

358. 

378. 

401 . 

425. 

452. 

29-30 

481. 

513. 
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Table  Ell 


Ex—le  3.  Pete  Pile  Printout 


0001 

1 

RIFLE 

TS0002  M  M 

1.35 

1.53 

0002 

2 

5  1  1 

0003 

61121 

1 .  5. 

78. 

17.  0.  0. 

0. 

0. 

0004 

1.71 

0.  4. 

0003 

61  122 

1 .  4. 

78. 

18.  0.  0. 

0. 

0. 

000ft 

1  .24 

0.  4. 

0014 

7 

TENTATIVE 

AVERAGE 

PREDICTION  RELATIONS 

0015 

7 

LOCATION  -  RIFLE  RANGE*  VICKSBURG*  MISS. 

0016 

7 

SOIL 

-  COLLINS  SILT  LOAM 

0017 

7 

VEGETATION 

-•  HERBACEOUS 

0018 

7 

RECORD  - 

1$  APR 

1951-31  DEC  1952 

3019 

"2 

13  4 

51 

0020 

9  8 

451 

77.77 

77.77 

77.77 

77.77 

77.77 

.02 

0021 

915 

451 

•  63 

2.69 

0022 

922 

451 

0023 

929 

451 

.02 

0024 

9  1 

551 

.03 

0025 

81 

10  5 

51 

3026 

9  8 

551 

.47 

.22 

3027 

915 

551 

.01 

0028 

922 

551 

0029 

929 

551 

3030 

9  1 

651 

0031 

9  8 

651 

.43 

1.79 

.02 

.64 

3032 

915 

651 

.47 

.40 

.13 

3033 

922 

651 

.86 

3034 

929 

651 

.40 

.24 

3035 

9  l 

751 

.81 

.72 

3036 

9  8 

751 

3037 

915 

751 

3038 

922 

751 

.28 

.05 

.06 

.01 

2.66 

3039 

929 

751 

3040 

9  1 

851 

3041 

9  8 

851 

3042 

915 

851 

.06 

1.61 

3043 

922 

851 

3044 

929 

851 

3045 

9  1 

9S1 

.58 

.03 

3046 

9  8 

951 

.71 

•  07 

2.  14 

3047 

82  ! 

21  9 

51 

3048 

915 

951 

.03 

3049 

922 

951 

.89 

3.82 

3050 

929 

951 
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Table  E12 


Example  3.  Teletype  Input 


FOR  INFORMATION  ON  RUNNING  THIS  MODEL  SEE  USER'S  MANUAL. 

WHEN  AN£"ERING  YES  OR  NO  QUESTIONS; 

USE  1  FOR  YES  AND  2  FOR  NO. 

IS  INPUT  RAINFALL  IN  CM? 

INPUT : 2 1 560 
?  2 

IF  YOU  WANT  OUTPUT  MOISTURE  CONTENT  IN: 

INCHES  ENTER  A  IS  CM  ENTER  A  2 
INPUT : 0 1 590 
?  1 

NAME  LAYERS  1  AND  2  SOIL  MOISTURE  CONTENT  FILES. 

EXAMPLE : XX 1 00 1 ; XX2 0 0 1 
INPUT: 01 620 
?  MC1001.MC2001 

DO  YOU  NEED  OUTPUT  MOISTURE  CONTENT  IN  PERCENT? 

INPUT: 01640 
?  1 

NAME  LAYERS  1  AND  2  PERCENT  MOISTURE  FILES. 

EXAM PLE : XX 1 0  0 1 ; XX2 0 0 1 
INPUT :0 1 680 
?  PM  1001 ;PM200 1 

IF  YOU  HAVE  AIR  TEMPERATURE;  ENTER  A  TJ  IF  NOT  A  9 
INPUT: 017 10 
?  9 

IF  YOU  HAVE  SNOW  DEPTH<M) ; ENTER  AN  SNI  IF  NOT  A  9 
INPUT: 01740 
?  0 

WHAT  IS  YOUR  INPUT  SOURCE? 

ENTER :  l  FOR  SURFACE  COMPOSITION  GROUP;  OTHERWISE  A  2 
INPUT : 0 1 780 
?  2 

ENTER  THE  NAME  OF  YOUR  INPUT  FILE 
INPUT :01810 
?  DTARIF 
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Table  E12  (concluded) 


ARE  MAX  AND  MIN  SOIL  MOISTURE  EQUATIONS  INCLUDED  IN  THE 
INPUT  DATA  FILE? 

INPUT  1 0 1 840 
?  2 

SITE  RIFLE  TS0002 
TEXTURE  M  M 

DENSITY  1.35  1.53 

SOIL  PROPERTIES  1.000 

0.000  0.000  0.000 

SOIL  PROPERTIES  1.000 

0.000  0.000  0.000 

MAX-MI N  MC  2.927  "  2.954  '  0.778  0.880  ' 

TENTATIVE  AVERAGE  PREDICTION  RELATIONS 
LOCATION  -  RIFLE  RANGE;  VICKSBURG;  MISS. 

SOIL  -  COLLINS  SILT  LOAM 
VEGETATION  -  HERBACEOUS 


RECORD 

-  13 

APR  1951-31  DEC  1952 

LAYER  =  1 

CONE  INDEX  EQ. 

A= 

1  . 3227  39227E+0 1 

B= 

-2.405180618E+00 

LAYER  =  1 

RATING  CONE  EQ. 

A= 

2.2378621 16E+01 

3= 

-5.421705970E+00 

LAYER  =  2 

CONE  INDEX  EQ. 

A= 

1  .362313798E+01 

B= 

-2.513599087E+00 

LAYER  =  2 

RATING  CONE  EQ. 

A= 

2.237869960E+01 

B= 

-5.400542561E+00 

DO  YOU  WANT  CONE  INDEX' FOR  FIRST  SOIL  LAYER?  SECOND  LAYER? 
RCI  FOR  FIRST  LAYER?  SECOND  LAYER?  SEPARATE  ANSWERS 
WITH  COMMAS.  EXAMPLE: FOR  RCI  ONLY  FOR  BOTH  LAYERS! 2; 2; 1; 1 
INPUT! 13 180" 

?  1  ;  1 ;  1  ;  1 

ENTER  NAMES  OF  SOIL  STRENGTH  FILES  TO  BE  USED; ENTER 
NODATA  IF  A  FILE  IS  NOT  TO  BE  USED. 

EXAMPLE:FOR  RCI  FOR  BOTH  LAYERSJNODATA;NODATA;XX1001;XX2091 
INPUT: 13220 

?  IC1001;IC2001;MR1001;MR2001 

STOP 

RUNNING  TIME:  88.6  SECS  I/O  TIME  :  39.7  SECS 


5.000 

1  .710 
4.000 

1  .240 


78  .000 

0.000 

78.000 

0.000 


17.000 

4.000 

18.000 

4.000 
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Table  E13 

Exaaple  3.  Predicted  Soil  Moisture  Content  In  Inches 


TENTATIVE  AVERAGE  prediction  relations 
LOCATION  -  RIFLE  RANGE.  VICKSBURG,  HISS. 
SOIL  ~  COLLINS  SILT  LOAM 
VEGETATION  >  herbaceous 
RECORD  -  13  APR  1951-31  DEC  1952 


son. 

MOISTURE 

IN/6 

IN 

6 

-  12  IN 

layer 

day 

APRIL 

1951 

8-14 

2.80 

2.68 

lb-21 

2.59 

2.50 

2.42 

2.35 

2.48 

2.40 

2.71 

2?~2» 

2.6l 

2.52 

2.44 

2.37 

2.29 

2.22 

2*15 

29-3U 

2.09 

2.02 

OAy 

may 

1951 

I-  / 

1.96 

1.90 

1.84 

1.79 

1.74 

1.69 

1*65 

e-i« 

1.61 

1.71 

1.74 

1.66 

1.57 

1.50 

1.43 

lb-21 

1.37 

1.31 

1.26 

1.21 

1.1? 

1.13 

1-10 

22-28 

1.9? 

1.05 

1.02 

1.00 

0,99 

0.97 

0*96 

29-31 

U.95 

0.94 

0.93 

CAy 

JUNE 

1951 

1-  ' 

P.92 

0.9? 

0.91 

0.91 

0.90 

0.90 

0.90 

8-14 

0.90 

0,98 

1.36 

1.31 

1.26 

1.21 

1.34 

11-21 

1.28 

1.38 

1.32 

1.40 

1.42 

1.35 

1.30 

22-28 

1.25 

1.20 

1.16 

1.13 

1.09 

1.07 

1.25 

2V-30 

1.33 

1.37 
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Example  3,  Predicted  Soli  Moisture  Content  in  Percent 


IDENTIFICATION  NUMBER:  RIFLE-  TSCTD02 
percent  mojsture  content 

6  12  IN  LAYER  - 


15-21 

22-28 

25-30 


2V76T  "2'9.5r — 29 

30.  f.  6  30.50  30 

29.52  29.39 — 29 

29.54  29.41  29 

30.31  '30.13  29 

31~'8  ~To  .83  3o" 

30.14  29.97  29 

31.21  T1.5T  -3l 

31.53  31.23  30 

30.91  30.68  31 


JANUARY 
TIT  30 , 84 
i.Jl  30.13 
.26  29* 14 

•  28  _  29.16 
.96 

FEBRUARY 
.61  30. #6 

.8?  ?9.67 

•27  31.00 

.97  31.46 

MARCH 
.36  31.08 


1952 
90.61 
29.96 
"29 . o3~ 
29.05 


30 .4"!  36.23 

_29 .81  29.66 

29 :6s~jr.w 

30.71  30.50 


195  2 

30.72  30.51  30> 32 

30.69  30.48  30*29 

30.76"  3T.39  "  SlTii” 
31.79  31.45  31.16 

4  aCo  “  —  —  -  — 


1952 

30.84 


30.62  30.41 


28.84 

30.49 

27.53 

28.43 

30.07" 

27.i_3_ 

28.03 

77.58 

26.74 

29.76 

29.15 

29.33  28.91 

28.74  "28.33" 

28.50 

27.93" 

26.34 

-TA-'t  1 — 

25.95 

APRIL 

25.57  26.66 

1952 

26.27  25.88 

25.49 

25.48  25.09  24.71  24, 3« 

23. 4«  25.ST~2T.AA  -  25.06 

23.57  23.22 

^  y~ 

22. b7  24.53  24.16  23.79 


18.48  17.83  17.20  20.16 

19.85  22.43  *23.88  "23fro 

24.52  23.74  22.97 

"  -  JUNE  " 

22,H  21.47  20.74  20.02 

rm  rsrr* — rra — ittst* 

13.76  1.3.38  13.03  J2.7i 

11.69  11.49"  11.31  11,14 

10.62  10.52 


23.97  24.17  23.81 

24.68- -24.31”  2379* 

1952”- . 

23.43  23.08  22*73 


21.08  20.36  19.65 

22.35"  21.-60"  237-36 

1952* - - 

19.32  18.64  17.98 

'*5.10  IT772 — ITT!* 

12.42  12.15  U.91 

10.99  10.86  "T&77T 
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Table  E15 


Example  3.  Predicted  Cone  Index 


IDENTIFICATION  NUMB'EHf  RIFLE  T5OO02 
CONE  INDEX 


6 

-12  IN 

I.  Iyer 

CAy 

JANUARY 

1952 

-  -  -■ 

— -  — 

-  7 ' 

l?5. 

M  %  If  iiiii  1  miii>  mu . . 

8-l« 

16' . 

154. 

156. 

159. 

161. 

163. 

165. 

15-21 

167. 

"1697 

”17  IT 

172‘. 

T74  V 

16T37 

1657  ~  — 

22-26 

167. 

168. 

170. 

172. 

174. 

151. 

154. 

29-31 

156. 

”  I59__ 

Tti: 

D*Y 

FEBRUARY 

1952 

1 "  7  -  ' 

14/. 

150. 

152; 

148V 

151. 

154. 

1>6. 

8-14 

158. 

161. 

163. 

165. 

151. 

154. 

156. 

1S-21 

146. 

141  . 

145. 

148. 

T5i. 

143. 

146. 

22-20 

141. 

145. 

148. 

142. 

139. 

142. 

146. 

day 

MARCH 

1952 

/ 

149. 

152. 

143. 

147. 

150. 

152. 

155. 

8  "l* 

161." 

159. 

146. 

152. 

159. 

1651 

l'l. 

1.5-21 

177. 

183. 

190. 

164. 

170. 

176. 

182. 

22-20 

‘  1541 

T60.  " 

■  W. — 

172. 

-179". 

res; 

1"92. 

29-31 

199. 

206. 

214. 

T 952" 

CAy 

APRIL 

l-  y 

222. 

230. 

239. 

216. 

224. 

232. 

241. 

a-iT" 

ii@nE£>fS 

?60. 

■KI  ONHH 

"2577" 

2247" 

”23H 

15-21 

241. 

251. 

261. 

271. 

281. 

275. 

286. 

22-20 

297". 

233. 

242, 

252. 

262 . 

272". 

2*2. 

29-30 

293. 

305. 

1952 “ 

BAy 

MAY 

1-  / 

3l7. 

266. 

276. 

287. 

298. 

310. 

322, 

8"-i4 

35c. 

38i. 

347. 

“379V 

413. 

“4517 

494 . 

15-21 

54  0. 

591. 

647. 

433. 

388. 

424. 

463. 

22-20 

452. 

"332. 

’  <T8’4  7 

309. 

336  • 

3657" 

300. 

29-31 

266. 

288. 

313. 

T952'  ‘ 

D*Y 

JUNE 

1  -  l 

341. 

37i. 

405. 

442. 

483. 

529. 

579. 

IF1* 

634. 

693 . 

760. 

VBTTT" 

7551 

"75151 

”7501 

15-21 

760. 

750. 

750. 

750. 

750. 

750. 

75o. 

22-28 

75". 

750. 

750V 

-750;  ‘ 

"  750 . 

750 . 

"750". 

29-3U 

75  . 

750. 
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6.1b 


5.74 


lb  3 

30  7  SENSOR  STUDY 

40  7  WEST  GERMANY  WURZBURG  AREA  18 

SO  7  SURFACE  COMPOSITION  COMPLEX  GROUP  0606 

60  7  SAND  WITH  FINES.  POORLY  DRAINED 

70  7 

80  83  1  1  58 


90  ' 

)  1 

1  58 

•  01 

•  1  3 

.31 

1.96 

100 

9  1 

I 58T  4.9 

1  .6 

-3.6  -4.3 

2.7 

6.7 

3.7 

no 

9  8 

158  .14 

•  30 

1.33  *01 

•  03 

.66 

.46 

IPO 

9  8 

1 58T  1.6 

2*6 

2.8  5.0 

2.3 

2.5 

2/0 

130 

915 

158  .02 

.02 

.10 

.64 

.09 

140 

915 

1  S8T  1 . 0 

0«8 

1.3  2.7 

2.0 

-2.3 

-1.9 

ISO 

915 

1  58S 

.  04 

.04 

160 

922 

158 

.01 

•  40 

•  01 

170 

922 

1 S8T  -0.6 

-3*4 

-3.8  -3.0 

-3.1 

-1.2 

-4.0 

180 

922 

158S  .04 

.02 

.02  .09 

.06 

.05 

.05 

190 

929 

158 

000 

929 

1  58T  -6.9 

-5-0 

-1  .4 

RIO 

929 

1C3S  .04 

•  04 

.04 

an  o 

9  1 

258 

.04 

.38 

1.18 

•  64 

a  so 

9  1 

2S8T  -3.2 

-1  .8 

1.2  -1.3 

3.0 

0.0 

-2.8 

:>4ii 

9  1 

2S8S  .04 

-04 

•  04 

.01 

.02 

.07 

2  SO 

9  8 

258  .40 

.60 

•  14 

1.08 

•  10 

.12 

260 

9  8 

258T  3.0 

7.2 

7.4  6.8 

10.6 

9*0 

14.7 

270 

9  8 

258S  .09 

280 

915 

258 

• 

CO 

.07 

•  53 

.97 

290 

915 

258T  13.3 

9.9 

CO 

• 

1 

• 

o 

-1.2 

1  .7 

3-3 

300 

915 

258S 

.01 

.07 

310 

922 

258  1.54 

.01 

1.74  1.00 

.03 

•  37 

320 

922 

258T  -0.6 

0.9 

4.0  8.2 

-0.6 

-3.4 

-1.7 

330 

922 

258S  .01 

.03 

340 

9  1 

358 

•  04 

•  59 

3S0 

9  1 

358T  -1.9 

1  *9 

2.0  2.6 

4.6 

2»  0 

-0.2 

.360 

9  8 

358  .06 

.17 

1.44  .09 

•  02 

.61 

370 

9  8 

358T  -0-4 

-2.2 

-4.2  -6.0 

-4.2 

-3*2 

-2.3 

380 

9  8 

358S  .01 

.02 

.12  .12 

.08 

.05 

.14 

390 

915 

358  .07 

400 

915 

358T  -3.6 

-2.7 

-2.3  -2.3 

-0.6 

-2.  1 

-3.7 

410 

915 

358S  .08 

.06 

.04 

420 

922 

358  .01 

.43 

430 

922 

358T  -4.6 

-3.3 

• 

n 

• 

o 

i 

5.4 

9.4 

9.6 

440 

929 

358  .55 

4bO 

929 

358T  10.6 

6*6 

1  .0 

460 

9  1 

458 

.45 

•  03 

.27 

470 

9  1 

458T  1.0 

2.  1 

8.4  10*7 

6.0 

6.5 

2.4 

480 

82 

10  4  58 
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Table  E17 


Example  4.  Teletype  Input 


FOR  INFORMATION  ON  RUNNING  THIS  MODEL  SEE  USER'S  MANUAL. 

'■WEN  ANSWERING  YES  OR  NO  QUESTIONS * 

USE  1  FOR  YES  AND  2  FOR  NO. 

IS  INPUT  RAINFALL  IN  CM? 

INPUT  SOI  560 
?  1 

IF  YOU  WANT  OUTPUT  MOISTURE  CONTENT  IN: 

INCHES  ENTER  A  1J  CM  ENTER  A  2 

INPUTS  0 1 590 
?  1 

NAME  LAYERS  1  AND  2  SOIL  MOISTURE  CONTENT  FILES. 

EXAMPLES XXI 00 1 *XX2001 

INPUTS  01 620 
?  CM1001 >CM2001 

DO  YOU  NEED  OUTPUT  MOISTURE  CONTENT  IN  PERCENT? 

INPUT  s  01 640 
?  1 

NAME  LAYERS  1  AND  2  PERCENT  MOISTURE  FILES. 

EXAMPLE : XXI 001  * XX20  0 1 

INPUTS01680 
?  MP1001 *MP2001 

IF  YOU  HAVE  AIR  TEMPERATURE*  ENTER  A  TI  IF  NOT  A  9 

INPUT  s  0 17 1 0 
?  T 

IF  YOU  HAVE  SNOW  DEPTH <M> *  ENTER  AN  SNI  IF  NOT  A  9 

INPUTS017A0 
?  SN 

WHAT  IS  YOUR  INPUT  SOURCE? 

ENTERS  1  FOR  SURFACE  COMPOSITION  GROUP*  OTHERWISE  A  2 

INPUTS01780 
?  1 
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Table  E17  (concluded) 
Example  4.  Teletype  Input 


ENTER  THE  NUMBER  OF  YOUR  SURFACE  GROUP 
INPUT : 0 1390 
?  0606 


SITE 
TEXTURE 
DENSITY 
MAX -MIN  MC 
ENTER  NAME  OF 
INPUT :0277e 
?  DT2RIF 


0606 
C  C 

1.55  1.65 

2 .420 
RAINFALL  FILE 


2.260 


0.500  0.540 


SENSOR  STUDY 

WEST  GERMANY  WURZBURG  AREA  18 
SURFACE  COMPOSITION  COMPLEX  GROUP  0606 
SAND  WITH  FINES*  POORLY  DRAINED 


ENTER  NAMES  OF  SNOW. FROZEN  SOIL  DEPTH  MODEL  INPUT  FILE 
EXAMPLE:  FRD001 
INPUT: 12190 
?  FR100DELETED 
?  FRD001 


DO  YOU  WANT  CONE  INDEX  FOR  FIRST  SOIL  LAYER?  SECOND  LAYER? 
RCI  FOR  FIRST  LAYER?  SECOND  LAYER?  SEPARATE  ANSWERS 
WITH  COMMAS.  EXAMPLE: FOR  RCI  ONLY  FOR  BOTH  LAYERS: 2. 2. 1.1 
INPUT : 1 3 1 80 
?  1.1. 1.1 

ENTER  NAMES  OF  SOIL  STRENGTH  FILES  TO  BE  USED. ENTER 
NODATA  IF  A  FILE  IS  NOT  TO  BE  USED. 

EXAMPLE: FOR  RCI  FOR  BOTH  LAYERS.'NODATA.NODATA.XX1001.XX2001 
INPUT: 13220 

?  IC1001.IC200 1 .CR1 001 . CR200 1 
STOP 

RUNNING  TIME:  175.1  SECS  I/O  TIME  :  51.8  SECS 
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Table  E18 


Example  4,  Input  Data  for  Freeze-Thcw  Model 


IDENT I r I  CAT  ION  NUMBER:  ,  0686 

FROZEN  DEPTH  OF  SOIL  (CRR£_L-1968  )_  DATA  INPUT  F I L.E*_ 

MEASURED  SNOW  DEPTH _ 


January  i95B 

density:  0-15  cm  LAYER,  1.55G/CMJ  lb-30  CM  LAYER,  1.656/CM 


oay 

1 

2 

3 

4 

b 

6 

7 

a 

9 

1" 

H 

1? 

13 

14 
lb 
1* 
1? 
l« 

2< 

?1 

2? 

"skew  pr  t  c  ip  -  ATfl 

DEPTH, M  ITATtC^CM  TEMP,  C  U 
• .  ..  u.0r‘  4,9n 

.  •  0.011  1.60 

u.01  -3.60 

0.0."  -4.30 

SOIL  TTOTSTURE 
-15  CM  LAYER.* 

25 . 79 

25.56 

25.33  — 

25.10 

■■■■simrHoisTURs —  ' 

lb-30  CM  LAYER, % 

22.69 

22.55 

22.41 

22.27 

>.  '  0T13 - 

0.31 

.....  1.96 

.  .  ,1  0.14 

.  '  0.30 

1.3’ 

:• . “  ■  o.oi 

..  •  0.0T 

0.66 
0.46 

•  0.02 

.  0.00 

"  rmr 

6.70 

3.70 

1.60 

2.60 

2.80 

24.87 

25.10 

25.43 

25.19 

“25.52 

25.84 

22713 

22.30 

- 22.61 

22.47 
“  '  22.77 

22.83 

Tnnr- 

2.30 

2.50 

2.00 

noo - 

0.80 

75  .'51 - 

25.38 

*25.70  - 

26.02 

“•25.79 

25.56 

22*69 

22.55 

22.83 

22.83 

“22.69  ’• 

22.55 

“8781 - 

.  .  0.0? 

. .  ■  0.1° 

.  4  0.64 

.  4  0.09 

.. .4  0.00 

"  1.30” 
2.70 

2.0" 
-2.30 
-1.90 
-  60 

2b  .33 

25.10 

24.87 

25.20 

24.-96 

24.73 

22741 

22.27 

'22.13 . 

22.44 

72.30 

22.16 

23 

• ,  ~7  ■  o . or 

'^7.  411 - 

— 24758 - 

22782 

54 

.2  0,00 

-3.80 

24.26 

21.88 

2b 

■ .  |9  0 , 4  n 

-3,00 

“24.88 

^2 .13 

26 

.  6  0.C(- 

-3.10 

24.65 

21.99 

2  7 

.5  0.00 

-1.20 

'  24.42 

71.85 

28 

..5  0. 01 

-4.00 

24.18 

21.71 

?9 

■  ,~t - 0.0  •  ~ 

'  -5T9TI - 

73795 

21757 

3'! 

.•4  0.0r 

-b.OO 

23.73 

21.43 

31 

..4  .j.m 

-1.40 

23.50 

21.29 
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APPENDIX  F:  GKAPHIC  DISPLAYS  OF  RESULTS 


1.  Appendix  F  contains  examples  of  graphic  displays  of  results 
for  the  first  three  examples  described  in  Appendix  E.  The  graphs  were 
run  on  the  CALCOMP  drum  plotter.  The  program  is  listed  in  table  FI. 

To  run  these  programs,  the  user  must  check  the  system,  program  compati¬ 
bility,  and  the  system  library  subroutines  called  by  the  program.  A 
run-time  input  by  teletype  is  shown  in  table  F2.  Figs.  F1-F6  are  graphs 
of  predicted  soil  moisture  content  and  Cl  for  the  first  three  examples 
described  in  Appendix  E. 

2.  Fig.  F7  illustrates  an  option  of  the  plotting  program.  Data 
for  each  of  two  (or  more)  years  are  plotted  on  one  axis.  Ten  or  more 
years'  output  could  be  plotted  with  this  option,  thereby  giving  a  long¬ 
term  visual  display  of  season  wetness  and/or  soil  strength. 


FI 
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MOISTURE  -  SECONO  LAYER 

I  soil  moisture  content  for  Example  1,  Appendix 


750 


o 


CD  in  CO  CM  *-» 

X3QNI  3N03 


CM 

ID 

CD 


u 
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TIME  IN  DAYS 

RIFLE  RANGE,  VICKSBURG,  MISS. 
CONE  INDEX  -  SECOND  LAYER 

Fig.  F2.  Predicted  Cl  for  Example  1,  Appendix 


co  csi 

3il fUSIOW  lN33iJ3d 
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o 

o 


Pd 

in 

O) 


CD  £ 
O  c 

cd  ce  S 
O  UJ5 


Q-  ccN* 

=  -*• 

o  w 

a:  cag¬ 
es  2« 

O  u 

•  On 
CD  UJO 

2OT- 

1  c 

•  Uii 
tS  a:  a 
oc 

=3  h-5 
CO  co» 

tn  wo 

ac  oB 

(_) 

►-«  O 

>  -J* 

Jr'S 

•  0“ 
UJ  (Ail 

O  s 

Z  »-“ 

a:  z°* 

ce  uj. 

U£ 
uj  oc 

-I  UJi 

u_  a.- 

M 

oc: 


750 


RIFLE  RANGE  VICKSBURG.  MISS.  G 
CONE  INDEX  -SECOND  LAY 

Fig.  F4.  Predicted  Cl  For  Example  2, 


Program  Listing  for  Graphic  Display  of  Results 

lUiiiUB>AX|Sj3. 

11  ■■  $L  !  B  •  L  AT  I  C » » 

:2'.slib»puott.  , 

13  DIMENSION  X<365)«  YC  <6f>).LABl<  i 6  ) , L AB2 ( 16 ) , L AB3 ( 6 > 

14 i  DIMENSION  mAxD<12),ICDM(2)»IDlS(2>#1DYR(2) 

15  COMMON  IDUMKY(40U),nUM(2U) 

16  DATA  JX/xH-/ 

1/  DATA  MaXD/31, ?8»3l, 30. 31.30,31 ,31. 3H.31, 30,31/ 
l»v  L-‘fV  =  3 

19-  P'xINT,»EnTfP  xmax  and  XMIN" 

20-  P" I  NT  * "  EXAMPLE:  400.-  0.” 

21 !  Rl-  ad.xmax.xmin 

22'.  S0X=(XMAX-XHlN)/8. 

23. :  SCaLEXi8.0/<XMAX-XMIN) 

24. ;  PKINT*"EnTER  YMAX  and  ymin" 

2t>-.  pR i nt * "  example:  ac..  o." 

26  ■  Br:AD,  YMAX#  YMIN 

2  I  •  SO'Y=  <  YmAX-YMIN)/6. 

28. 1  Si-ALEYi6.('/{  YMAX-YMIN) 

29, -  N :  1 

30.  CALL  PLOIS(LDEV) 

31'  PR  I  NT# "TYPE  IN  DATA  PILE  NAME." 

32',  RtAD  14  ,  F  IlNAM 
33n  !•*  F0HmAT<A6) 

34-  PkjNT«"EnTER  NUM8ER  OF  YEARS  TO  BE  PLOTTED." 

3  b  RL-aD.NPL 

36,  PR  I  NT  * "DO  YOU  WANT  ALL  THE  YEARS  PLOTTED  ON  ONE  GRAPH?" 

3/'  PRINT#"  TYPE  1  FOR  YES;  2  FOR  NO" 

3B(,  RlaD.NOP 

39-  CALL  OPEnF(I.PILNAM) 

40'.  PKJNT. "ENTER  FIRST  LINE  OP  HEADING  CENTERED  IN  48  SPACES." 

4  lii  READ  1  2  •  LaBI 

42  PRINT.  1,  "Enter  second  line  op  heading  IN  T:-*E  SAME  MANNER." 

43(1  PRINT. "SPACE  AT  THE  END  OP  FIRST  HEADING  LINE  IS  RESERVED" 

44-i  PRINT.  "BY  THE  PROGRAM  FOR  THE  YEAR." 

4b.i  RLaDi»2,lAB2 
46(i  i  -2  FORMAT  ( 16A3) 

47c  P^INT. "ENTER  VERTICAL  AXIS  TITLE  CENTERED  In  18  spaces." 

48c  RLaD  1  4.LAB3 
49t  1  A  FORM aT ( 6 a3 ) 

50(.)  DC  2  jO  K  =  1 .  NPl 
51-  L  s  1 
52(-  NI=1 
53:  Mb  s  0 

54.  PRINT. "ENTER  BEGINNING  DAY  AND  MONTH  AND  ENDING  DAY  AND  MONTH 
55o  PRINT. "FOR  THIS  YEAR.  EXAMPLE : BEG  I N  l3  APR,  AND  END  28  DEC." 
56n  print."  ENTER:  13.4,?8.12" 

5'/'.  RLAD,  IBEG.  IMB.  IEND.  IME 

58.  jDa-7 

59  IV)  lo5  I  J  =  1 . 5 

60o  I*  ( IBEG.LE. JDA) IBEG=< I J-l)*7*i 
610  PdBfcG.LE.JDA)  GO  TO  106 
62<>  1  -5  JDA  =  jDA«7 

63.1  l v. 6  PRINT. "ENTER  YEAR  AS  FOLLOW" 
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Table  FI  (continued) 

Program  Listing  for  Graphic  Display  of  Results 

64.1  pHlNT»"2i3LANKS»  4 NUMBER" 

6b  RLAD  lb, 1 D Y  R  < 1 ) » 1DYR<2  ) 

66,'  l  j  r0RMAT(2A3) 

67  D)  19b  H=1.2 

68  RlaD(1»1.  )  (DUM(NN)  .NN  =  1»9) 

6V,,  1  •  FORMAH  Ad) 

70,.-  DO  190  J  =  l»  6 

71"  Jl (M.E0.1.AND. J.LT. 1MB)ME=ME*MAXD(L) 

72,  1MM.EO.  l.AND.  J.LT.  1MB>GO  TO  1  88 
73  ■  (L.GT . IME1G0  TO  190 
74 )  Rl aD ( 1 > 1 . ) DUM ( 1 ) 

7b',  1F(l.EO.  IME)HE=IEND*ME 

76  It  ( l . bO . I  ME  >  GO  TO  107 

77  HbsME*MAxD(L ) 

78-  if (M.EO. l.ANO.J.EO. 1MB )N8EG=NI=N1*1 BEG-1 
7 V i  1'  7  RfcAD(l,2).<Y(I>.t=Nl.ME) 

80u  2  FORMAT  (1i,X#7F7.Q) 

81u  lf (M.EQ.l .AND. J.EQ. IMB)NI=NI*MAXD(L)-JBEG*1 
82<  H  (M.EO. l.AND. J.EQ. IMB)  GO  TO  189 
83'  1  >8  N I  =N  1  *MAXD  (  L  ) 

84  109  L=L*1 

8b  1 ,0  CONTINUE 
86.  RtAD(l,l  ) DUM ( 1 ) 

8  7;!  1  '5  CONTINUE 

88.  I- (NOP. 60. l.AND. N.GE. 2)  GO  TO  45 

89  >  If  (N-1)16»16.18 

900  10  CALL  PLOT (0. .-30, , -3? 

91  GO  To  8C 

92  1('  IF(MOD(N»2).EQ.O)GO  TO  70 

93  .  C*lL  PL0T<15. ,-12. ,-3) 

94  GO  TO  80 

9b  '.  7  CALL  PLOT(0.  .12.  ,-3) 

96i  GO  To  82 

97  8  CALL  PL0T(2. .3. .-3) 

9b;  6.  CONTINUE 

99-  CAlL  Ax  I Sl3 ( 0 • .0. .12HTIME  IN  DAYS»-12».14,B,,-1,0»XM1N*SCX»1,.1»0) 
10"0  CALL  AXlSl3(0. .0..LAB3. 18. 0.14,6. .-1.1, YMIN.SCY.1., 1,0) 
lOlO  Call  SYMBOL<0.71,-n.90,0.14»LABl,0.,48) 

1020  call  symBol(7. 10,-0.90.0.14. idyr.o, .6) 

1010  call  symbolic .71.-1 .18,0.14, lab2,o., 48) 
io-io  call  latIc(0.,o. .6.2,8. ,i.d 
lO‘-0  15  continue 
1060  call  ploTu.,o..3> 

1080  Do  40  i =nreg,mE 
10-0  x ( i ) = i 

11  ,0  X(I)S(x(I)-xm1N)*SCAlEX 

lllD  T(1)S(Y(I)-YMIN)*SCALEY 

1120  I PEN=3 

lD(J  CONTINUE 

nr)0  DO  50  I=Ne)EG,ME 

1160  iF(Yd).LE.O)  GO  TO  65 

11/0  f  F ( Y  < I J.LE.6.0)  GO  TO  60 
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Table  FI  (concluded) 

Program  Listing  for  Graphic  Display  of  Results 

UtiU  '.'AH.  PLOT  <X(  1 )  #Y(  1  )  >  |  PEN } 

11*U  I pENs3 
1 2(  0  -0  TO  5  • 

12U‘  1  PEN  =  3 

l2i?0  !,0  TO  5 

IP’-  *>0  CALL  PlOT(X(  I  >,Y(  I  ),  I  PEN  > 

1240  ;pEn»2 

l2t)U  0  CONTINUE 

1260  lFtN.E0.NPL>  GO  TO  9: 

1 2  0  N  =  N*1 

I2fl0  i)  0  CONTINUE 

12->U  >.)  CALL  PL0T(0.«0.  t  3) 

13mC  CALL  PLOT ,0. #999) 
l3j  (•  ‘jTOP 
132  0  t  nD 
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Table  F2 


Example  of  Teletype  Input  for  Graphic  Display  of  Results 


ENTER  XMAX  AND  XMIN 

EXAMPLE:  400.4  0. 

INPUT: 002 10 
?  400.40* 

ENTER  YMAX  AND  YMIN 

EXAMPLE:  48.4  0. 

INPUT: 00260 

?  60 .40. 

****FR0M  PLOTTER  -  PLEASE  INPUT  A  PLOT  FILE  NAME 
INPUT: 00 145 
?  AMSP01 


TYPE  IN  DATA  FILE  NAMt. 

INPUT: 00320 
?  PM2001 

ENTER  NUMBER  OF  YEARS  TO  BE  PLOTTED. 

INPUT :00350 
?  2 

DO  YOU  WANT  ALL  THE  YEARS  PLOTTED  ON  ONE  GRAPH? 

TYPE  1  FOR  YES 1  2  FOR  NO 
INPUT : 00380 
?  2 

ENTER  FIRST  LINE  OF  HEADING  CENTERED  IN  48  SPACES. 
INPUT: 004 10 

?  RIFLE  RANGE.  VICKSBURG.  MISS.  GROUP  0808 


ENTER  SECOND  LINE  OF  HEADING  IN  THE  SAME  MANNER. 
SPACE  AT  THE  END  OF  FIRST  HEADING  LINE  IS  RESERVED 
BY  THE  PROGRAM  FOR  THE  YEAR. 

INPUT: 00450 

?  PERCENT  SOIL  MOISTURE  -  SECOND  LAYER 

ENTER  VERTICAL  AXIS  TITLE  CENTERED  IN  18  SPACES. 
INPUT: 00480 
?  PERCENT  MOISTURE 
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Table  F2  (concluded) 

Exaaple  of  Teletype  Input  for  Graphic  Diaplay  of  Results 

ENTER  BEGINNING  BAY  AND  MONTH  AND  ENDING  DAY  AND  MONTH 
FOR  THIS  YEAR.  EXAMPLE  *  BEG I N  13  APR.  AND  END  28  DEC. 

ENTERS  13.4.28. 12 
INPUT  S00570 
7  13.4.31.12 

ENTER  YEAR  AS  FOLLOW 
2BLANKS. 4NUMBER 
INPUT: 00652 
?  1951 

ENTER  BEGINNING  DAY  AND  MONTH  AND  ENDING  DAY  AND  MONTH. 
FOR  THIS  YEAR.  EXAMPLE t BEG 1 N  13  APR.  AND  END  28  DEC. 

ENTERS  13.4.28. 1 2 
INPUT  S00570 
?  1.1.31.12 

ENTER  YEAR  AS  FOLLOW 
25LANKS . 4NUMBER 
INPUTS  00650 
?  1952 

STOP 

RUNNING  TIME:  21.7  SECS  I/O  TIME  <  12.6  SECS 
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APPENDIX  G:  PROCEDURES  FOR  CONVERTING  MOISTURE  CONTENT  AND 
CONE  INDEX  TO  TERMS  REQUIRED  FOR  THE  AIRFIELD  CONSTRUCTION  EFFORT  (ACE)  MODEL 


X.  The  Airfield  Construction  Effort  (ACE)  model  requires  ao  input 
the  moisture  condition  (wet  or  dry)  and  the  California  Bearing  Ratio  (CBR) 
of  the  upper  layer  of  surface  material.  Wet  or  dry  conditions  of  the  soil 
and  CBR  can  be  derived  from  percent  aoil  moisture  content  and  Cl  values, 
respectively,  that  are  output  from  the  SMSP  model.  As  previously  discussed, 
the  SMSP  model  outputs  data  for  the  0-  to  15-cm  and  15-  to  30-cm  layers, 
respectively.  Since  the  properties  of  the  15-  to  30-cm  layer  more  nearly 
approximate  the  average  properties  of  an  upper  layer  of  several  feet  of 
soil  thickness,  which  is  of  direct  pertinence  to  the  required  airfield  con¬ 
struction  effort,  this  second  layer  is  the  only  layer  considered  in  the 
conversion  procedures  that  follow. 

Conversion  of  Cl  to  CBR 

2 .  The  CBR  for  the  ACE  model  can  be  derived  from  the  following 
equation: 

CBR  =  0.02  C£  (i5-30  cm  layer) 

This  equation  was  derived  from  an  analysis  of  data  presented  in  plate 
43  of  reference  16,  and  closely  approximates  the  relation  shown  in  fig. 

4.11  of  reference  17.  The  following  table  shows  Cl  class  values  derived 
from  the  above  equation,  that  are  equivalent  to  CBR  classes  used  in  the 
ACE  model: 

Subgrade  CBR  Equivalent  Subgrade  Cl 

3-9  150-<500 

10-20  500-1000 
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It  should  be  noted  that  a  Cl  less  than  150  or  a  CBR  less  than  3  indi¬ 
cates  that  construction  effort  is  not  considered  because  the  soil  is 
too  soft  to  permit  efficient  operation  of  equipment. 

Conversion  of  Moisture  Content  to  Wet  or  Dry  States 

3.  The  soil  moisture  content  at  the  plastic  limit  (PL)*  is  con¬ 
sidered  to  be  the  boundary  line  between  the  dry  and  wet  states  of  the 

soil.  Soils  with  moisture  contents  above  PL  act  as  plastic  materials,  adhere  to 
equipment,  and,  in  general,  are  more  difficult  to  work  than  soils  with 
moisture  contents  below  PL. 

4.  The  wet  or  dry  conditions  of  the  soil  can  be  determined  from 
one  of  the  two  following  procedures,  depending  upon  the  availability  of 
PL  information  for  the  soil: 

«i.  PL  of  the  soil  is  known: 

%  MC  (15-  to  30-cm  layer) 

_ from  SMSP  Model _  Soil  Moisture  Condition 

<PL  Dry 

>PL  Wet 

t>.  PL  of  soil  is  not  known:  Use  data  from  the  following  tab¬ 
ulation  to  estimate  the  PL  of  the  soil,  and  then  compare  the  PL  with  the 

*Plastic  limit  units  are  in  percent  moisture  content. 

G2 

26 Z 


percent  moisture  content 

(15-  to  30-cm  layer) ,  as 

in  ji  above 

determine  the  dry  or  wet 

condition  of  the  soil: 

Region 

Soil  Type 

Average 

PL 

Tropical 

Sand 

16 

Silt 

26 

Clay 

33 

Laterite 

22 

Temperate 

Sand 

18 

Silt 

23 

Clay 

22 

Arid 

Sand 

16 

Silt 

21 

Clay 

21 

Data  in  the  above  table  were  obtained  from  table  18  of  reference  8. 
Soil  samples  from  1066  sites  were  used  to  derive  values  of  average  PL. 
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